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Abstract 

This article reviews the current theoretical and experimental status of the 
field of muon decay and its potential to search for new physics beyond the 
Standard Model. The importance of rare muon processes with lepton flavor 
violation is highly stressed, together with precision measurements of normal 
muon decay. Recent up-to-date motivations of lepton flavor violation based 
on supersymmetric models, in particular supersymmetric grand unified theo- 
ries, are described along with other theoretical models. Future prospects of 
experiments and muon sources of high intensity for further progress in this 
field are also discussed. 
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I. INTRODUCTION 



The muon was discovered in 1937 by Neddermeyer and Anderson (Neddermeyer and 
Anderson, 1937) in cosmic rays. The mass of the muon was found to be about 200 times the 
mass of the electron. The discovery was made just after Yukawa postulated the existence of 
the 7r meson, a force carrier of the strong interaction, in 1935 (Yukawa, 1935). But, it was 
demonstrated in 1947 that the muon did not react through the strong interaction, and thus 
it could not be the Yukawa n meson (Conversi, et al, 1947). The famous comment by Rabi 
"Who ordered that ?" might indicate how puzzling the existence of a new lepton was at that 
time. By then, it was known that the muon decays into an electron and a neutral particle. It 
was thought that if the muon were a heavy electron, it would also decay into an electron and 
a 7-ray. The first search for /i + — > e + 7 was made by Hincks and Pontecorvo in 1947 using 
cosmic-ray muons (Hincks and Pontecorvo, 1947). Its negative result set an upper limit of 
less than 10%. This was the beginning of the search for lepton flavor violation, i.e. violation 
of the lepton number conservation for each generation. In 1948, the continuous spectrum 
of electrons was established, suggesting a three-body decay with an electron accompanied 
by two neutral particles (Steinberger, 1948). Soon after, the search for a neutrinoless muon 
nuclear capture process (fJ>~N — > e~iV, where N is a nucleus capturing the muon) was 
also carried out, but with a negative result (Lagarrigue and Peyrou, 1952). Such searches 
were significantly improved when muons became artificially produced at accelerators. In 
1955, the upper limits of B(n — > ej) < 2 x 10~ 5 (Lokonathan and Steinberger, 1955) and 
B(/jl~Cu — > e~Cu) < 5 x 10~ 4 (Steinberger and Wolfe, 1955) were set at the Columbia 
University Nevis cyclotron. 

After the discovery of parity violation, it was suggested that the weak interaction takes 
place through the exchange of charged intermediate vector bosons (Feynman and Gell-Mann, 
1958). In 1958, Feinberg pointed out that the intermediate vector boson, if it exists, would 
lead to /i + — > e + 7 at a branching ratio of 10~ 4 (Feinberg, 1958). The absence of any 
experimental observation of the fi + — > e + 7 process with B(/j, — > cy) < 2 x 10~ 5 led directly 
to the two-neutrino hypothesis (Nishijima, 1957; Schwinger, 1957) in which the neutrino 
coupled to the muon is different from that coupled to the electron, and the /i + — > e + 7 process 
would be forbidden. The two-neutrino hypothesis was verified experimentally at Brookhaven 
National Laboratory (BNL) by confirming muon production, but no electron production 
from the scattering of neutrinos from pion decays (Danby, et al, 1962). This introduced 
the concept of the separate conservation of lepton flavors, electron number (L e ) and muon 
number (L^). 

Our understanding of modern elementary particle physics is based on the Standard Model 
(SM), which is a gauge theory of the strong and electroweak interactions. The formulation of 
the SM is based on many theoretical developments of gauge theory in the 1960s and 1970s. 
Since then, the SM has confronted various experimental tests and has had amazing success in 
explaining all measurements so far. In the minimal version of the SM, where only one Higgs 
doublet is included and massless neutrinos are assumed, lepton flavor conservation is an 
automatic consequence of gauge invariance and the renormalizability of the SM Lagrangian. 
It is the basis of a natural explanation for the smallness of lepton flavor violation (LFV) in 
the charged lepton process. 

In extensions of the minimal SM, however, LFV could occur from various sources. In 
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fact, in most cases, new physics or interactions beyond the SM would predict LFV at some 
level. LFV processes with muons of major interest are such as /i + — > e + 7 , fi~—e~ conversion 
in a muonic atom (fi~N — > e~N), fi + — > e + e + e~ , and muonium to anti-muonium conversion 
(Mu — Mu conversion). Historical progress in various LFV searches in muon and kaon decays 
is shown in Fig.|l], in which the experimental upper limits have been continuously improved 
at a rate of about two orders of magnitude per decade for about 50 years since the first LFV 
experiment by Hincks and Pontecorvo. The current LFV searches with muons have reached 
a sensitivity on the order of 10~ 12 — 10~ 13 in their branching ratios. In general, searches for 
rare processes could probe new interactions at high energy. For example, in the four fermion 
interaction, the LFV branching ratios could be scaled by (m^/m^) 4 , where mx is the mass 
of an exotic heavy particle responsible for the LFV interaction and mw is the mass of the W 
gauge boson. Thus, the present sensitivities for LFV searches in muon decays could explore 
a mass scale of several 100 TeV, which is not directly accessible by the present accelerators. 

Recently, considerable interest in LFV processes has arisen based on supersymmetric 
(SUSY) extensions to the SM, in particular supersymmetric grand unified theories (SUSY 
GUT). Since the three gauge coupling constants of the strong, weak, and electromagnetic 
interactions, which were measured at LEP and SLC, are shown to be consistent with the 
assumption that they are unified to a single 577(5) gauge coupling constant at the order of 
10 16 GeV scale in SUSY SM, the SUSY GUT model becomes a very attractive candidate 
for physics beyond the SM. In the SUSY models, in general, there is a new source of flavor 
mixing in the mass matrices of SUSY partners for leptons and quarks. This will induce 
LFV processes for charged leptons. The predictions of the branching ratios depend on flavor 
mixing in the mass matrix of sleptons, which are supersymmetric partners of leptons. In 
the SUSY-GUT scenario, the flavor mixing in the slepton sector is naturally induced at the 
GUT scale because leptons and quarks are in the same GUT multiplet (Hall, et al., 1986). 
It has been shown (Barbieri and Hall, 1994; Barbieri, et al., 1995a) that the surprisingly 
large top-quark mass determined recently has an impact on calculations of the branching 
ratios of fi + — > e + 7 and fi~—e~ conversion in SUSY GUT. The predictions are as large as 
just one or two orders of magnitude lower than the present experimental limits. 

There is considerable evidence for the existence of neutrino masses and their mixing based 
on the experimental results of the solar neutrino deficit (Cleveland, et al., 1998; Fukuda, et 
al, 1996; Hampel, et al, 1999; Abdurashitov, et al., 1996; Fukuda, et al., 1998a) and the 
atmospheric neutrino anomaly (Fukuda, et al, 1998b). Since neutrino oscillations imply that 
lepton flavor is certainly not conserved, LFV processes in muon decays are also expected 
to occur. In non-SUSY models, however, the neutrino mixing introduces only small effects 
on /i + — > e + 7 . For example, the branching ratio of /i + — > e + 7 is on the order of 10~ 50 
for a difference of the neutrino mass squared of AmJ ~ 10~ 3 eV 2 , since it is suppressed 
by (AmJ/m^) 2 (Petcov, 1977; Bilenky, et al, 1977). The situation changes drastically in 
SUSY models. In SUSY models with a neutrino-mass generation mechanism of the see-saw 
type, the Yukawa coupling constants among the Higgs doublet, lepton doublets and right- 
handed neutrinos could induce large flavor mixing effects in the slepton sector (Borzumati 
and Masiero, 1986; Hisano, et al., 1998a; Hisano and Nomura, 1999a). The resulting LFV 
rates can be as large as, or even larger than, the experimental upper bounds, depending on 
various parameters, especially on the Majorana mass of the right-handed neutrino. In such 
a case, the same Yukawa coupling constants for the right-handed neutrino are responsible 
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for both the neutrino oscillation and the LFV processes of charged leptons. 

Thus, there are many theoretical scenarios under which the predicted branching ratios 
for the muon LFV processes can be close to their present experimental upper bounds, and 
therefore they could be accessible and tested by future experiments. 

Experimentally, there has been much progress. First of all, several new experimental 
results have been obtained using the high-intensity muon beams now available, and on- 
going and proposed experiments are aiming at further improvements. Furthermore, in the 
long term, new attempts to create high-intensity muon sources have been initiated, based 
on the ideas arising from the collider project. The muon beam intensity aimed at such 

a new muon source would be about 10 12 — 10 13 / u ± /sec, several orders of magnitude higher 
than that presently available. With this increased muon flux, significant improvements in 
experimental searches can be anticipated. 

In this article, we review the current theoretical and experimental status of the field of 
muon decay and its potential for new physics beyond the Standard Model. We highly stress 
the importance of rare LFV processes of muons, especially within the framework of SUSY 
models. In addition, we cover precision measurements of normal muon decay. There have 
been many excellent review articles on muon decay and lepton flavor violation (Frankel, 
1975; Scheck, 1978; Vergados, 1986; Engfer and Walter, 1986; Depommier, 1987; Van der 
Schaaf, 1993; Depommier and Leroy, 1995). However, reflecting current renewed interest, 
this article has been written to bring up to date recent topics on muon decay and physics 
beyond the SM. The phenomenology and experimental status of some of the important muon 
processes are described in detail. 

This article is organized as follows. In Chapter 2, we give a short summary of the SM 
and the muon properties within the SM. In Chapter 3, LFV is discussed based on various 
theoretical models, including the SUSY models. Chapter 4 deals with the current status of 
precision measurements in normal muon decay, such as the muon lifetime, the Michel decay 
spectrum and its asymmetry, and e + polarization. In Chapter 5, the phenomenology and 
status of the most recent experiments on various lepton flavor violating muon decay modes, 
such as fi + — > e + 7 , ji + — > e + e + e~ , n~~—e~ and /i - — e + conversions in a muonic atom, and 
Mu — Mu conversion, are described. In Chapter 6, prospects on future experiments and 
high-intensity muon sources are briefly discussed. Some useful formulas are collected in the 
Appendices. 



6 



II. BASIC PROPERTIES OF THE MUON IN THE STANDARD MODEL 



A. Muon in the Standard Model 



1. The Standard Model 

The current view of elementary particle physics is based on a gauge theory of quarks 
and leptons. In the Standard Model (SM), three fundamental interactions, the strong, 
electromagnetic and weak interactions, are described as a SU(3)c x SU(2)l X U(l)y gauge 
theory. Quarks and leptons are classified in three generations. The three quarks of |e 
electric charge are the up quark (u), charm quark (c) and top quark (t), and those of — 
charge are the down quark (d), strange quark (s) and bottom quark (b). Correspondingly, 
there are three charged leptons of — e electric charge, electron (e), muon (/i) and tau (r), 
and three species of neutrinos of neutral charge are introduced, i.e. electron neutrino (u e ), 
muon neutrino (i/^) and tau neutrino (u T ). These six quarks and six leptons are given in 
Table | 

TABLE I. Quarks and leptons in the Standard Model 



Electric charge 


1st generation 


2nd generation 


3rd generation 


¥ 


u 


c 


t 




d 


s 


b 











-1 


e 




T 



In the SM, ferminonic fields, gauge fields, and a SU (2) L doublet Higgs field are introduced 
as an elementary field. They are listed along with their quantum numbers in Table [TI], where 
the SU(3)c, SU{2) L and U(l)y gauge fields are denoted as G M , and B^, respectively. 
The subscripts of L and R represent left-handed and right-handed chirality projections 
(Pl = (1 — 75)/2 and Pr = (1 + 7s)/2), respectively. H represents the Higgs doublet field. 
The suffix i (= 1 — 3) for the quark and lepton fields is the generation index. The SU(2)l 
doublet fields, such as qi L , k L , and H, have field components given by 



^=[Z)- H -{>)- « 

The SM Lagrangian, Csm, consists of three parts, which are for the gauge interaction, 
the Higgs potential, and the Yukawa interaction. It is given as 

£"SM £gauge ^Higgs ^Yukawa- (2) 

The Lagrangian for the gauge interaction, C gaugei is presented by 

£ 3 auge= £ F$ + £ lJ lL{R)1 »V ^ iL{R) + \V »H\\ (3) 

SU{3)c,SU(2) L ,U(l) Y quarks, leptons 

where F nu is the gauge-field strength, and X> M is a covariant derivative defined as 
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TABLE II. Quantum numbers of elementary fields in the minimal Standard Model. The 
SU(3)c, SU{2)l representation and U(l)y charge are given. 





G> 


4. 




QiL u iR di R 


hi 




H 


SU(3) C 
SU(2) L 
U(1) Y 


8 
1 




1 

3 




1 
1 




3 3 3 
2 11 

1 2 1 

6 3 3 


1 

2 

i 

2 


1 
1 

-1 


1 

2 

i 

2 






2>m 


= d^ + ig s 


y Gl + ig^Al + ig'QyB, 






(4) 



for the representations with SU (3)c-triplet and SU (2) L -doublet and Qy-U{1)y charge quan- 
tum numbers. g s , g, and g are the gauge coupling constants for SU(3)c, SU (2) L , and U(1)y, 
respectively. A a (a = 1 — 8) is the Gell-Mann matrix for a SU(3) group, and r a (a = 1 — 3) 
is the Pauli matrix for a SU(2) group. The terms for a singlet representation for either the 
SU(3)c or SU{2) L gauge groups are absent in the definition of V^. 
The Lagrangian for the Higgs potential, £m g gs, is given by 

C Hl99S = -(-fi 2 \H\ 2 + X\H\ 4 ). (5) 

For fi 2 > 0, the Higgs field develops the following vacuum expectation: 

where v = jj/\/X (= 246 GeV). The physical Higgs mass is given by mu = v^2At>. After 
electroweak symmetry breaking, the SU{2)l and U{1)y gauge fields form a massless photon, 
and massive W ± and Z° bosons. At the tree level, their masses are given as raw = \gv and 

m z = |V ' 9 2 + Q 2 v - The SU(3) gauge boson, gluon, remains massless. 
The Yukawa interaction part of the Lagrangian is given by 

^■Yukawa = {y e )ijH' f e'i R l jL + {yd)ijH ] d iR qj L + {y u )ijH^ui R qj L + H.c, (7) 

where (y e )ij, (Vdjij, and (y u )ij are the Yukawa coupling constants for the charged leptons, 
the down-type quarks, and the up-type quarks respectively, and, 

H = i T 2 H*=(^y (8) 

Substituting the vacuum expectation value for the Higgs field, the Yukawa interaction 
in Eq.(^) generates the mass terms for quarks and leptons, as follows: 

C mass = -[e- iR (m e ) i:j e jL + d lR {m d ) ij d jL + ul^m^^Uj^ + H.c, (9) 

where (m e )^ = -(y e ) i:j (v/<\/2), {m d ) i:j = -{y d ) i:j {v / ^/2) and (m M )y = -{y u ) i:j {v / ^/2). Each 
mass matrix is diagonalized by unitary transformations for the left-handed fermions and the 
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right-handed fermions with the same charge. Since the unitary matrices for the left-handed 
up- type quark and the left-handed down-type quark are generally different, flavor mixing is 
induced in the charged weak current interaction for quarks. It is given by 

C w - qg = -^(w iL ^(VcKM)ijd jL W+ + d lL ^{VcKMr 3l u 3L W-) , (10) 

where the {VcKM)ij represents the flavor mixing matrix for the quark sector, i.e. Cabibbo- 
Kobayashi-Maskawa (CKM) matrix (Kobayashi and Maskawa, 1973). In this equation and 
hereafter, the quark fields are presented in the mass-diagonalized basis. 

On the other hand, the charged lepton mass matrix in Eq.(|9]), equivalently the lepton 
Yukawa coupling constant, is fully diagonalized by unitary transformations on the lepton 
doublet fields (k L ) and the lepton singlet fields (e JR ). In the mass-diagonalized basis, the 
charged weak current interaction for leptons remain diagonal, as follows: 

Cwve = -^_{vi L ^e iL W; + ei Ll »v iL W-\ (11) 

In the above basis, the lepton flavors can be defined for each generation, and are thus 
conserved. They are the electron number (L e ), the muon number (L^), and the tau number 
(L T ), as defined in Table |TJ. 

TABLE III. Assignment of lepton flavors, the electron (L e ), the muon number (L^), and the 
tau number (L T ). 





e 






*V 


r 


v T 


e+ 


V~e 


M + 


*V 


r+ 






+1 


+ 1 














-1 


-1 














w 








+1 


+1 














-1 


-1 






















+1 


+1 














-1 


-1 



2. Interaction of the muon in the Standard Model 

At the tree level of the SM Lagrangian, the muon has three gauge interactions, namely 
those with the photon, the and Z° bosons, and also the Higgs interaction. They are 
given by 

C = eJlY^A, - A={tj^ l1 ^ l W+ + Wl^ L W-) 

- \jg 2 + 9' 2 {wr{-\ + sin 2 6 w )(jl l + MY sin 2 Q wm \zl - ^JlfiH, (12) 

— 1) 

where the Weinberg angle 9 W is defined by sm9 w = g / \j g 2 + g' 2 , and also e — g sin 6w 
the tree level. H denotes the physical Higgs boson field. In addition to the electromagnetic 
interaction, the second and third terms describe, respectively, the charged weak-current 
interaction mediated by the boson and the neutral weak-current interaction mediated by 
the Z° boson. The other charged leptons, electron and tau, have the same gauge interaction 
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as the above, and the coupling constant to the Higgs boson is proportional to their lepton 
masses. 

The muon decay in the SM is described by a charged weak-current interaction mediated 
by the W 1 * 1 gauge boson. The four fermion interaction is given by 

Cpermi = --^=[lyy /i (l - ~f 5 ) fxe-f^l - 75K + Iv/^l - ^efx^y^l - 75)^], (13) 

where Gf is the Fermi coupling constant. At the tree level of the SM, this is given by 

G ' = *& (14) 

where my/ is the W ± boson mass. This interaction describes the normal muon decays, 
fx + — > e + u f JU fl and fT — > e~77 e z/ M . 

Lepton universality is a fundamental property of the gauge interaction. The universality 
in charged weak currents has been tested from the combination of leptonic and semi-leptonic 
decays of r, and leptonic decays of fx, 7r and K mesons. The constraints on the ratio of the 
charged weak current coupling constants for electron (g e ) and muon (g^) are obtained from 
the ratios of B{t~ — > e~V^i> T ) / 'B(t~ — > fx~V^v T ) and B(ix~ — > e~V e ) / B(ix~ — > fx^V^). For 
the ratio of the charged weak current coupling of tau (g T ) and g^, T(r _ — > e~V e v T ) /Y[t~ — > 
fx~V^u T ), r(r _ — > z/ r 7r _ )/r(7r _ — > fTV^) and T(r~ — > v T K~)/Y{K~ — > fx~u^) were used. 
These ratios give a test of the equality of the relevant coupling constants at the level of 
0.1 - 0.5 %, except for those decays involving kaons, where the sensitivity is at the 2% 
level (Pich, 1997). A similar test can be performed for leptonic decays of the W bo- 
son at Tevatron and LEP. A recent measurement of leptonic branching ratios of W gives 
B(W -> ev) = 10.9 ± 0.4%, B(W -> fxv) = 10.2 ± 0.5% and B{W -> rz/) = 11.3 ± 0.8% 
(Particle Data Group, 1998). Thus, the sensitivity for the coupling ratio is a few %. The 
lepton universality for neutral weak currents has been tested at the Z° boson pole. At LEP 
experiments, the measurements of partial leptonic widths, leptonic forward-backward asym- 
metries, tau polarization and its angular distributions are sensitive to the vector and axial 
vector coupling constants for different lepton species. At SLC, the left-right asymmetry and 
forward-backward left-right asymmetry are also measured using a polarized beam. Lepton 
universality is now treated at the level of 5-10% for vector couplings and less than 0.2% for 
axial vector couplings (Pich, 1997). The difference in sensitivity is due to the fact that the 
lepton vector coupling is very small compared to the axial coupling. 



3. Neutrino mass and mixing 

Although the minimal SM does not allow massive neutrinos, there is increasing evidence 
for the masses of neutrinos from the solar-neutrino deficit (Cleveland, et ai, 1998; Fukuda, et 
al, 1996; Hampel, et ai, 1999; Abdurashitov, et ai, 1996; Fukuda, et ai, 1998a) and 
the atmospheric- neutrino anomaly (Fukuda, et ai, 1998b). If the solar- neutrino deficit 
is explained by neutrino oscillations, the mass-square difference is in the range of Ami ~ 
10~ 6 - 10~ 5 eV 2 for the MSW solution (Wolfenstein, 1978; Mikheyev and Smirnov, 1985) 
or Ami ~ 10~ n eV 2 for the vacuum oscillation (or "just-so" oscillation) (Glashow and 
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Krauss, 1987; Glashow, et al, 1999). Also, the atmospheric-neutrino anomaly suggests that 
the mass-square difference between the muon neutrino and the tau (or sterile) neutrino is 
on the order of Am 2 ~ 10_3 - 10 ~ 2 eV2 (Fukuda, et al, 1998b). In addition, the LSND 
(Liquid Scintillator Neutrino Detector) experiment has reported the V^ivy) — > 17 e (z/ e ) oscil- 
lation, suggesting |m 2 — ml e \ ~ 10 _1 — 10° eV 2 (Athanassopoulos, et al, 1998), although 
the KARMEN experiment has seen no evidence of such an oscillation (Eitel, et al, 1999). If 
neutrino mixing is the true interpretation of the anomalies, the SM must be extended. On 
the other hand, there are direct upper bounds on the neutrino mass for each species, which 
is 15 eV/c 2 for the electron neutrino, 170 keV/c 2 for the muon neutrino, and 18.2 MeV/c 2 
for the tau neutrino (Particle Data Group, 1998). Recently, improved measurements of the 
electron neutrino mass have been reported with a much better sensitivity of a few eV/c 2 
(Lobashev, 1998; Otten and Weinheimer, 1998). 

It is possible to accommodate the Dirac mass terms for the neutrinos if SU(2) singlet 
fields of the right-handed neutrinos v iR {% = 1 — 3) are included in the minimal SM field 
contents. Then, the following interaction can be added to Eq.(^): 

C VR = (y u ) ij H^V- iR l jL + H.c, (15) 

where (y u )ij is the Yukawa coupling for neutrinos. If it is very small, the small masses of 
neutrinos can be explained. For example, the Yukawa coupling constant should be O(10~ n ) 
for a neutrino mass of 1 eV/c 2 . Note that the total lepton number is conserved in this 
scenario, whereas the lepton flavor could be generally violated. 

A more natural explanation for the small neutrino masses is provided by the "see-saw 
mechanism" (Yanagida, 1979; Gell-Mann, et al., 1979). In this scenario, the Majorana mass 
term is also included for the right-handed neutrino, 

£» R = i.yu)ijH ] Vl R l jL - ■^i R {M^) i jVj c R + H.c, (16) 

where the charge-conjugation field is defined as ip c = —ip T C~ l and the charge-conjugation 
matrix (C) satisfies C~ 1 7 M C = —•y fJ,T . (Mg)ij is the right-handed Majorana neutrino matrix. 
Substituting the vacuum expectation value for the Higgs field, the neutrino mass terms 
become 

1 



C = - 2 ( (m)°,v- m ) ^ m ; m ° r j ^ ( „« . j + H,., (17) 

where the Dirac mass term is {mE>)ij = —(y v )ij(v/y/2). When the Majorana mass scale is 
much larger than the Dirac masses, the lighter neutrino masses are given by 



£v mass ~ --{ViL) C {m u )ijV jL + H.C. (18) 

and 

m v = m T D {M R )- l m D . (19) 

For example, if Mr is 10 15 GeV and the Dirac mass is in the range of 100 GeV, then the 
neutrino mass becomes naturally O(10~ 2 ) eV. 
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By diagonalization of the charged lepton and neutrino mass matrices, lepton flavor mixing 
is induced in the charged weak current interaction, as follows: 



r — y 

>~>Wve — ;= 



9_ 
V2 



ViLYiYuNS. 



w; + e iLl »{v MNS ) mL w; 



(20) 



where the (VMNs)ij is the flavor mixing matrix for the lepton sector, i.e. Maki-Nakagawa- 
Sakata (MNS) matrix (Maki, et ai, 1962). Note that as in the case of quarks, the lepton 
fields in Eq. (]20f) are written in the mass-diagonalized basis. The {Vmns)^ matrix element 
represents neutrino mixing which can be studied by neutrino oscillations. For a review on the 
neutrino mass and mixing, please see other references (Bilenky and Petcov, 1987; Fukugita 
and Yanagida, 1994; Mohapatra and Pal, 1998; Fisher, Kayser, and McFarland, 1999). 



B. Static properties of the muon 

1. Mass and lifetime 

The mass and lifetime of the muon are the fundamental inputs of the SM. The muon 
mass is given by (Particle Data Group, 1998) 

m M = 106.658389(34) MeV. (21) 

It is derived from the ratio of the muon mass to the electron mass, m^/rrie, which is measured 
in a muonium (/i + e - atom) with QED corrections (Cohen and Taylor, 1987). 
The experimental value of the muon lifetime is 

r M = 2.19703(4) x 10~ 6 sec. (22) 

In the framework of the SM, the muon lifetime (r M ) is related to the Fermi coupling constant 
(Gp), by including QED corrections, as follows (Kinoshita and Sirlin, 1959; Marciano and 
Sirlin, 1988): 



-i _ G 2 F m 5 m 2 3 ml 



aim,.) ,25 



^ 1927^ m^ 5 77% L 2tt 4 J 



where F(x) = 1 — 8x + 8x — x — Ylx lnx. m^ and m e are masses of the muon and the 
electron, respectively. The a value at the m^ scale, a(m^), is given by 



lnp) + -L « 136. (24) 

Air V m.- / hir 



aim^) 1 = a . 

Sir x -m P / bix 



From Eq . (|23|) , the Fermi coupling constant of Gp = 1.16639(1) x 10~ 5 GeV -2 is determined 
(Particle Data Group, 1998). The higher order two-loop corrections to the muon lifetime 
have been calculated (Van Ritbergen and Stuart, 1999). 
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2. Magnetic moment 



Since the muon is a Dirac particle, the g factor of its magnetic moment is 2, if radiative 
corrections are ignored. A deviation from 2, namely g — 2, is very important to investigate 
quantum corrections. The present experimental value of a u = (g n — 2)/2 is given by (Particle 
Data Group, 1998) 

a e * p = 11659230(84) x 1(T 10 (±7ppm). (25) 

A new experiment, E821, is on-going at Brookhaven National Laboratory (BNL) and the 
experimental error is aimed to be improved by a factor of 20. The result from the first run 
gives a™ p = 1165925(15) x 10" 9 (±13 ppm) (Carey, et al. } 1999a). 

Theoretically, this quantity is calculated very precise (Hughes and Kinoshita, 1999). 
The correction is divided into higher-order QED corrections, hadronic contributions and the 
electroweak (EW) contributions. A recent update of theoretical calculations for them gives 
(Czarnecki and Marciano, 1998) 

a% ED = 11658470.56(0.29) x 10~ 10 , (26) 

a hadron = 673 .C;( 6 . 7 ) x ( 27 ) 

a™ = 15.1(0.4) x 10~ 10 . (28) 
By adding them, the SM prediction is 



n SM _ n QED i hadron , EW 



= 11659159.6(6.7) x 10~ iU . (29) 

The theoretical prediction is in good agreement with the experimental value. 

To calculate the QED correction, the fine structure constant is needed as an input. 
Eq . (pBD was estimated by using the fine structure constant obtained from the Quantum Hall 
effect, which gives a~ 1 (qH) = 137.03600370(270). It is consistent with that determined 
from the electron anomalous magnetic moments (a e ) by assuming a theoretical evaluation 
based on the SM. It is a _1 (a e ) = 137.03599959(38), where the experimental values of the 
electron anomalous magnetic moments are a e ^- = 1159652188.4(4.3) x 10 -12 and a^+* = 
1159652187.9(4.3) x 10~ 12 . 

Although the electron g — 2 factor is better determined than the muon g — 2 factor, the 
latter is much more sensitive to short-distance physics. For example, the EW correction to a e 
is O(10 -14 ) level compared the a EW quoted in Eq. (p8|) . Therefore, the muon g — 2 experiment 
is much more important to search for the effects of new physics. For example, the on-going 
experiment E821 is expected to put very strong constraints on super symmetric (SUSY) 
models, since SUSY models could contribute to the muon g — 2 factor significantly when the 
slepton, charginos, and neutralinos exist in the mass range of a few hundred GeV (Lopez, et 
ai, 1994; Chattopadhyay and Nath, 1996; Moroi, 1996; Carena, et ai, 1997; Gabrielli and 
Sarid, 1997). 
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TABLE IV. Decay modes of muons. 



Decay mode 


Branching ratio 


References 


H~ -> e'VeUfj, 


~ 100% 




\i~ -> e"I7 e i/ M 7 


1.4 ± 0.4% ( for E y > 10 MeV) 


Crittenden, et al. (1961) 




(3.4 ±0.4) x 10~ 5 


Bertl, et al. (1985) 


/i~ -> e - !^ 


< 1.2% 


Freedman, et al. (1993) 


M~ -> e"7 


< 1.2 x lO" 11 


Brooks, et al. (1999) 


/i~ — > e~e~e + 


< 1.0 x 10~ 12 


Bellgardt, et al. (1988) 


/i~ — > e~77 


< 7.2 x 10~ n 


Bolton, et al. (1988) 



C. Decay modes of the muon 



The measured decay modes of muons are — > e~V e v^ (Michel decay), \i~ — * e~V e i>^y 
(radiative muon decay), and /i~ — > e~T7 e u ll e + e^ . The branching ratios for these modes 
and the upper bounds on the other exotic decay modes at 90% confidence level (C.L.) 
are summarized in Table |V|. Although these branching fractions have been measured in 
experiments with positive muon decays, they are listed for negative muons by assuming CP 
invariance. Since \i~ — *> e~V e v^ cannot be clearly separated from {i~ — ► e~z7 e f M with a 
soft photon, the branching ratio for the radiative decay is shown for E 1 > 10 MeV. There 
is no evidence of lepton flavor violating processes, such as /i - — ► e~7 , [i~ — > e~e~e + and 
pT — ► e _ 77. Also, an upper bound is set for those with |ALj| = 2, such as fi~ — > e'v^U^ 
decay, which is allowed if the lepton flavor is conserved multiplicatively instead of additively. 



1. Normal muon decay 

In the SM, the normal muon decay is described by the V— A interaction. In extensions 
to the SM, the energy spectrum of a decay electron (positron), its angular distribution if 
muons are polarized, and its spin polarization in pT — > e~V e v^ (or p + — > e + i> e U ll ) decay 
are sensitive to the type of interaction on muon decays, including new possible interactions 
besides the V— A interaction. If the general four fermion interactions with no derivatives 
are assumed, the muon differential decay rate is given with a few parameters by (Fetscher 
and Gerber, 1998) 

^dxaZ^ = i^^W* 2 - *§ (Fis(x)±P,cos9 e F AS (x))(l + P e (x,e e )-C), (30) 

where W efl = (m'f l +ml)/(2m fl ), x = E e /W efJi and x = m e /W etl {= 9.7 x 10~ 3 ) < x < 1. E e is 
the energy of the e^. m e and m M are the masses of the positron and the muon, respectively. 
The plus (minus) sign corresponds to decay. 9 e is the angle between the muon 

polarization (P^) and the electron (or positron) momentum, and ( is the directional vector 
of the measurement of the spin polarization. P e (x, 9 e ) is the polarization vector of the e^. 
The functions Fjs{x) and Fas(x) are the isotropic and anisotropic parts of the e energy 
spectrum, respectively. They are given by 
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2 

Fis(x) = z(l - x) + -p(4x 2 — 3x — Xq) + r]x (l - x), (31) 
9 



1 



3 



2 



^as(x) = -eV^ 2 - 4 1 - * + ~^x -3 + (Jl-x 2 - 1)) , (32) 



3 



where p, r], £, and 5 are called Michel parameters (Michel, 1950; Bouchiat and Michel, 1957). 

In the SM, these parameters are given asp=~, rj = Q, £=1 and S = |. When the 
electron (positron) polarization is not measured and Xq is ignored, the differential branching 
ratio in the SM in Eq. (|30|) leads to a simple form of 



(fT^ -> e ± z/I7) mlG F 2 



dxd cos 6 P 1927r 3 



x 2 [(3-2x) ±P M cos# e (2x- 1)]. (33) 



Fig.|2] shows the e + energy spectrum in p + — > e + u e V ll decay in the SM, for the cases of 
cos 9 = 0, cos# e = +1, and cos# e = —1 with 100% polarized positive muons. As can be 
seen, the spectrum is high at x fa 1, and the sign of asymmetry changes at x = 1/2. 

If the SM is not assumed, the muon lifetime in Eq. (p3|) should be replaced by (Scheck, 
1978; Fetscher and Gerber, 1995; Pich and Silva, 1995) 

G 2 P m!r„,m^ , m P „ / m? s 32, 3 x m 2 _ mt. 



r- 1 = — ^ 

" 1927T 3 



F(-f) + 4^G(-f)--(p- ,,,, 
3 _ aim,,) ,25 



x 



(1 + 54 )|1 + ^ ( T-^ (34) 

where G(x) = 1 + 9x — 9x 2 — x 3 + 6x(l + x)lnx. Regarding radiative corrections, since it can 
be assumed that the SM contribution dominates in the normal muon decay process, those 
based on the SM, as in Eq. (|23|) , are used. From Eq.([34]), the correction from the i] parameter 
is proportional to 0(m e /m^), whereas that from the p parameter is very small, in the order 
of 0(ml/mV). Since the present experimental accuracy of the rj parameter is about 1%, 
the uncertainty from the i] correction is about on the order of 10~ 4 to the estimation of the 
muon lifetime in the non-SM case. 

When the spin polarization of e + (e~) in the p + — > e + v e V^ (p~ — > e'v^pe) decay is 
detected, P e (x,9 e ) in Eq. (|30|) can be measured. It is given by 



Pe(x, 9 e ) = P Tl ■ ^f** + Pt 2 ■ -^4 +P L -^, (35) 

where z is the direction of the e momentum, and P^ is the muon spin polarization. P^, 
Pti, and Pt2 are, respectively, the e ± polarization component parallel to the momentum 
direction, that transverse to the e momentum within the decay plane, and that transverse 
to the e momentum and normal to the decay plane. A non-zero value of the triple T-odd 
correction, Pt2, would imply violation of time- reversal invariance. They are given by 

P (r 9\- P^sme e F T1 {x) 

PTl[X,V e ) — — — — — — , (dbj 

F IS (x) ± P M cos9 e F AS (x) 
P (r 9)- ^sinfl e F T2 (s) 

y T2 {x, v e — — — — — --, (6( 

F IS (x) ± Pf, cos9 e F AS {x) 
p / n n ±Fip(x) + P^ cos 9 e F AP (x) 

r L {x, V e - — — — -—— , (-38 J 

F IS (x) ± P M cos9 e F AS (x) 
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TABLE V. Experimental values of some of the Michel decay parameters. 



Michel parameter SM value Experimental value Sensitive observables 

p 3/1 0.7518 ±0.0026 Fis 

r) -0.007 ± 0.013 F IS and P T1 

5 3/4 0.7486 ± 0.0038 F AS and P L 

£ 1 1.0027 ± 0.0084 F\ s and P L 

f' 1 1.00 ±0.04 P L 

f" 1 0.65 ±0.36 P L 



t Only the product of £P M is measured. 



where the ± sign corresponds to fj, decays, and 



F T1 (x) = ^|-2[e" + 12(p-|)](l- x)x - 3 V (x 2 - xl) + V "(-3x 2 + Ax - sg) j, (39) 

F T2 (x) = i^ 2 -^§ {3^(1 " *) + 4^""^}' ( 4 °) 
F IP (x) = ^-Jx^-xl ^{-2x + 2 + yjl-xft + 4£(<f - |)(4ar - 4 + 0^8) }, (41) 
^p(^) = ^"(2x 2 - a - a$ + 4(p - ^)(4x 2 - 3x - x 2 ) + 2r/'(l - x)x ^. (42) 

where 77", {a /A), and [0 /A) are newly defined Michel parameters (Kinoshita and 

Sirlin, 1957b; Fetscher and Gerber, 1998). In the SM, ( = f" = 1 and 77" = (a'/A) = 

= 0. 

The muon-decay Lagrangian for the general four-fermion couplings with ten complex 
parameters is expressed as (Fetscher, et al, 1986) 



9RR(eRiy eL )(^ L ^R) + 9RL(e R iye L )(^ R ^L) 



+9 S LR(e L ^eR)(y; L m) + 9^1^)^^) 

+9RR(e R YveR)0^ R l^R) + g RL {e R ^v eR ) 

+9 V L R (zLYVeL)(^ R l^R) + 9LL&Ll' i VeL)fa L lp»L) 
T T 

+ 9 -^(e R <j^v eL )(v; R <j^ L ) + ^(e L a^v eR ){vr L <j^ R ) + H.c. ] , (43) 

where there is a normalization condition of 

\(\9rr\ 2 + \9L\ 2 + \9 S r L \ 2 + \9 S LR \ 2 ) + (\9 V rr\ 2 + \ 9 V LL \ 2 + \ 9 V RL \ 2 + \ 9 V LR \ 2 ) 

+ 3(\9 T r L \ 2 + \9 T lr\ 2 ) = 1- (44) 

Note that in the V— A interaction of the SM, g\ L = 1 and all the others are zero. 
The Michel decay parameters of p, 77, £ and 5 are given by 
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P = ^(\9 S nn\ 2 + \9 S ll\ 2 + \ 9 L ~ ^II? + \gL ~ 2g T LR f) + 3 -(\g v RR \ 2 + I^J 2 ) (45) 
V = Ir49rr9ll + 9 V ll9rr + gktfut + ^Ir) + gUsfa. + 6fl&)] (46) 
e = \(\9L\ 2 - \9 S rr\ 2 ) - \(\9 S L r\ 2 - \9 S r L ?) + (\9 V ll\ 2 - \9 V rr\ 2 ) + S(\g V LR \ 2 - \g V RL \ 2 ) 

+H\9 t lr\ 2 - \g T R L \ 2 ) + ±Re{g S L R g T R l - gLgm) (47) 
= ^(\g s LL \ 2 - \g S R R \ 2 + \g S R L - 2g T RL \ 2 - \g s LR - 2 5 Ll 2 )) + \{\g v LL \ 2 - \gU 2 ) W 

Table summarizes the present knowledge of the Michel decay parameters (Particle Data 
Group, 1998). Precise measurements of the Michel decay parameters would place constraints 
on various new physics, which would induce a small deviation from the V— A couplings. The 
current constraints on the general four-fermion couplings are summarized in Fetscher and 
Gerber (1998). 

Among many theoretical models which can be studied by precise measurements of normal 
muon decay, one example is the left-right symmetric models of SU(2)l x SU(2) R x U(l) 
symmetry (Herczeg, 1986; Langacker and Sankar, 1989). In this model, heavy right-handed 
gauge bosons (two charged W R and one neutral Z') are added. In general, they are mixed 
with each other, and form mass eigenstates W^ 2 - 



W l\_( cosC -wa.£\(Wt 
Wr)~ U^sinC e^cosCAwf 



(49) 



where ( is a mixing angle, and u is a CP-violating phase. One of the models is a manifest 
SU (2) l x SU (2) R x U (1) model, which has the same gauge couplings and CKM quark mixing 
matrix elements for both the SU{2)l and SU(2) R sectors. It, however, can receive a strong 
constraint from the K^—Ks mass difference, yielding a limit of the mass of W R , m R > 1.4 
TeV. The constraint can be much weaker in general cases (Langacker and Sankar, 1989). In 
the manifest SU{2) L x SU{2) R x U(l) model, the Michel parameters of p and £ would be 
related to £ and the masses of W^ 2 as follows: 

^(1-^P)=C 2 , (50) 

i(i-0 = C 2 + (^) 4 , (51) 

where m 12 are the masses of W^ 2 - The direct lower limit for the W 2 mass of m 2 > 720 
GeV/c 2 is obtained from DO (Abachi, et al, 1996). Together with the previous measurement 
of p and £ parameters (Jodidio, et al., 1986), the constraint on the ( — m R plane is given in 
Fig! 



2. Radiative muon decay 

The spectrum of the radiative muon decay, p^ — > e^vV'y , has been calculated by sev- 
eral authors (Kinoshita and Sirlin, 1957a; Fronsdal and Uberall, 1959; Eckstein and Pratt, 
1959). Within the framework of the V— A interaction, the differential branching ratio of 
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the radiative muon decay, where final electron (positron) and photon are emitted to energy 
intervals of dx and dy at solid angles of dfl e and dQ 1 respectively in the muon rest frame, is 
expressed by 



^^(3dx—dQ e dQ y [F(x, y, d) =f pP^ ■ p e G(x, y, d) =F P^ ■ p~ ( H(x, y, 



(52) 



P M is the muon polarization vector; p e and are the momenta of an electron (positron) 
and a photon in the muon rest frame respectively; p e and p 7 are their unit vectors defined 
by p e = p e /\p~e\, P-y = Pj/\p~i\ respectively; (3 is defined as (3 = \p e \/E e ; d is given by 
d=l- (3p e ■ p-f] x and y are normalized electron and photon energies, x — 2E e /m fl and 
y = 2Ey/m lt in the muon rest frame. From the four-body kinematics, the ranges of x and y 
are given by 



(1 



2y/r < x < 1 + r for < y < 1 — y/r, 
y) + r/(l — y) < x < 1 + r for 1 — y/r < y < 1 — r, 



(53) 



where r = (m e /m M ) 2 . -F(x, (i), G(x, y, d) and if (x, y, d) in the SM are given in Appendix A. 

The decay probability distribution in the x-y plane is shown in Fig.[|. The probability 
is high for an energetic e ± with a soft photon, namely x ~ 1 and y pa 0. In the soft-photon 
limit (y — > 0), the distribution has an infrared singularity which is canceled by the radiative 
correction of the Michel decay. 

The photon spectrum is obtained by integrating over the electron (positron) energy and 
angle variables. By neglecting the terms suppressed by m e /m M , it is given by (Kuno, et at, 
1997a) 



dB{^ -> e ± z/77 7 ) 1 



dyd cos 9 1 y 
where J+(y) and J-(y) are defined by 



J+(y)(l ± P» cos^ 7 ) + J_(|/)(l =F P^cos# 7 ) 



(54) 



-(i 

+ (21n 



I/) (31n- 



i-y 



r 
3 A 



a 
6n 



l-yf (31n 



1 



+(21n- 



^-^)(1 
r 12 A 
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-31n- 



+ 7)(l-y) 



+ (-41n^ + f)(l-y) 
r 3 



(55) 



(56) 



and # 7 is the angle between the muon spin polarization and the photon momentum. The 
photon spectrum for unpolarized muons is shown in Fig.|5|. Note that at the maximum 
photon energy (y ~ 1), the photon distribution is approximately given by (1 + P^cos^ 7 ) 
for the ji + — > e + v e V^f decay, because J+(y) has a first-order term in (1 — y), but J-(y) 
only contains the second and higher order terms. This fact is important for the suppression 
of accidental background in a fi + — > e + 7 search using polarized muons, as mentioned in 



Section VA5 
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In generalized interactions, the differential branching ratio of fi^ — > e ± vV r ) decay has 
been calculated (Lenard, 1953; Behrends, et al., 1956; Fronsdal and Uberall, 1959). Here, 
the spectra of e ± and a photon depend not only on the Michel parameters of p and 5 in 
the normal muon decay, but also on an additional parameter, rj, which should be zero in 
the V— A interaction in the SM. Also, the asymmetry of in /x ± — > e^vv^ from polarized 
muons is parameterized by another parameter, £ • k (Fetscher and Gerber, 1995). Measure- 
ments of these parameters would give additional constraints on the four fermion coupling 
constants (Eichenberger, et al, 1984). Time- reversal violation in radiative muon decay was 
also discussed (Pratt, 1958), but it was concluded that the T-odd effects have to include 
either the e ± polarization or those terms suppressed by the electron mass. 
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III. LEPTON FLAVOR VIOLATION AND PHYSICS BEYOND THE STANDARD 

MODEL 



In the minimal SM with vanishing neutrino masses, lepton flavor is conserved separately 
for each generation. This is not necessarily true if new particles or new interactions beyond 
the SM are introduced. In this section, we discuss various theoretical models with LFV in 
the charged-lepton processes. In particular, we mention those in which LFV effects could 
be large enough to be detected in on-going or future experiments on fi + — > e + 7 decay, 
fi + — > e + e + e~ decay, fi~—e~ conversion, and other LFV processes. Among the theoretical 
models which predict observable LFV effects, SUSY models have recently received much 
attention. In SUSY models, the origin of LFV could be interactions at a very high energy 
scale, such as the GUT scale or the mass scale of a heavy right-handed Majorana neutrino 
that appears in the sea-saw mechanism. Searches for rare muon decays, thereby, could 
provide a hint for physics at a very high energy scale. In the following, we first explain 
the effective Lagrangians for various muon LFV processes of |AL;| = 1, and then discuss 
LFV in SUSY models and other theoretical models. Finally, the muon LFV processes with 
\ALA = 2 are discussed. 



A. Effective Lagrangians for lepton flavor violating processes 



The effective Lagrangians for muon LFV processes of | ALj| = 1, such as ji + — > e + 7 decay, 
/i + — > e + e + e~ decay, and fi~—e~ conversion in a muonic atom, are discussed. Possible 
LFV contributions can be grouped into two types: photonic interaction and four fermion 
interaction. 

First, the effective Lagrangian for ji + — > e + 7 process is given by 



C 



4G F 
'75 



m fJ ,A R (x R (T f "'e L F ta , + m a A L ^ L (j t11 ' 'e R F uu + H.c. 



(57) 



e R l and //+ -> e^ 7 



where A R and A L are coupling constants that correspond to /i + 
processes, respectively. 

For /i + — > e + e + e~ decay and fi~—e~ conversion, off-shell photon emission also con- 
tributes. The general photonic fj, — e transition amplitude is, then, written as 



M, 



photonic 



q u q u 



eA;(q)u e (p e ) (f m (q 2 ) + l5 f M0 (q 2 )hAg" V ~ ~t) 



+(/Mi(g 2 ) + 75^i(g 2 ))-^ g 



m. 



(58) 



where p u and p e are the fj,~ and e~ four momenta, and q = p u — p e is the four-momentum 
transfer. The electromagnetic form factors (/eo, /ei, /mo an d /mi) are functions of q 2 . For 
ji + — > e + 7 decay, only /ei(0) and /mi(0) can contribute, whereas all of the four form factors 
could contribute to fi + — > e + e + e~ decay and fi~—e~ conversion. The coupling constants A R 
and Al are related to the dipole form factors as 
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A R = -^j(r E i(o) + r M1 m, m 

A L = ^j(r E1 (0)-fMM). (60) 

The direct four fermion interactions could introduce /i + — > e + e + e~ decay and fi~—e~ con- 
version, in addition to the photonic fi — e transition in Eg. (|58|) . For the fi + —>■ e + e + e~ decay, 
the general four fermion couplings are given by 



prion— photo *±\J E 

^ 3e ~^/2 



gi(^Be L )(e R e L ) + g2{^L,e R ){e L e R ) 



+9b{'pRl^e R ){ell^e L ) + g 6 (Wl fl e L )(e^'j^e R ) + H.c. 

where the Fierz rearrangement for the four fermion operators is used. For the /i" 
conversion process, the relevant interactions are written as 



(61) 



prion— photo _ l ~ r F 
^fi—e conv \p2, ' 



_ (9LS{q)ZLHR + 9RS(g)e R fi L )qq 

q=u,d,s.. 

+ {9LP{ q )el^R + 9RP( q )e^P>L)qi5<l 

+ (9LV(q)eI / y fl ^L + gRV^eRl^HRjll^q 

+ {9LA{ q )ell^^L + 5 , i?A(g)e^7 /i / i i?)97^75<? 

+^(9LT( q )eZ^ u ^R + gRT(q)e^^ u ^L)qcr U uq + H.c. 



(62) 



where gLXia) an d gRX(q) are the coupling constants for the left-handed and right-handed 
lepton currents, respectively, and X = S,P, V, A, T represent scalar, pseudoscalar, vector, 
axial vector, and tensor interactions, respectively. Here, the flavor- changing quark currents 
are not included. The four fermion coupling constants introduced in Eqs. flBTD and (|62|) are 
determined by specific contributions in some theoretical models beyond the SM. For exam- 
ples, they are box diagrams in supersymmetric models, tree diagrams of Z , super symmetric 
models with R— parity breaking, and others. 

The f E0 and /mo form factors contribute to off-shell photons and not to real photon 
emission. Therefore, they vanish in the q 2 — > limit. They could be rewritten by 



W) = ^f E0 (q 2 ), (63) 

m u 

fuo{q 2 ) = ^hio(q 2 ), (64) 

where f E o(q 2 ) and /mo(q' 2 ) are finite at q 2 — ► 0. If these transitions are induced by loop 
diagrams including heavy particles, f E o(q 2 ) and Jmo(<? 2 ) are regarded as slowly varying 
functions of q 2 . One example of such diagrams is shown in Figfl(a). In such a case, these 
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form factors could be translated into additional contributions which should be added into 
the corresponding four-fermion coupling constants in Eqs.(^) and (|62|) . Those additional 
contributions are 



y/2 e 



2 



Ag, = Ag 5 = -(/* (0) + /^ (0)), (65) 

4<_j p ill 



A^ 4 = Ag 6 = ^-^(/lo(O) - /mo(0)), (66) 



for g 3 , g A) g 5 , and g e , correspondingly, and 



2 



&9lv(u) = -2Ag LV{d) = -^^(f E0 (0) + / M0 (0)), (67) 



2 \[2 e 2 ~ 

A^v (u) = -2A^ W (d) = J (/bo(0) - /mo(0)), (68) 

for 9LV(d), 9rv{u) and g^y^) respectively. 

If these form factors are generated by penguin diagrams with a photon coupled to an 
internal line of a light fermion, as seen in Fig.|6](b), fEoil 2 ) and /mo(<? 2 ) have a logarithmic 
dependence of q 2 that is cut off by the light-fermion mass. Such examples are a model 
with a doubly charged Higgs boson and the SUSY models with i?-parity violation, which 
are discussed later in Section |III Oj. The logarithmic factor could enhance the rates of 



fi + — ► e + e + e~ decay and [i~—e~ conversion, but not that of /i + — ► e + 7 decay. 

If fEi{<l 2 ) and fh[i(q 2 ) dominate, the following simple relations among the branching 
ratios of [i + — > e + 7 , fi + —>■ e + e + e~ and fi~—e~ conversion can be derived: 

iW-eTi) ^ l. 5(/i+ _ e+7); (69) 



T(/iTj — > capture) 200 

5(/i + -> e+e+e") ~ y^5(/i + ^ e+7) " (70) 

These relations hold in some models of SUSY GUT, which are discussed in the next sub- 
section. Regarding fi~—e~ conversion, more detailed discussions on Eq. (|69|) , including the 
nuclear dependence, are given in Section |V C 1| . 



B. Supersymmetric models 

1. Introduction to supersymmetric models 

Phenomenological applications of SUSY theories have been considered since the late 70's 
in connection with the naturalness problem (or the hierarchy problem) in the SM. The SM 
model can be regarded as being a low-energy approximation of a more complete theory, and 
is replaced by this at a high energy scale. Supposing that this high energy scale is close to the 
Planck scale (~ 10 19 GeV), the quadratic divergence appearing in the radiative corrections 
to the Higgs scalar mass becomes problematic, because a precise fine tuning between the 
bare mass and the radiative corrections must be necessary to keep the electroweak scale well 
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below this high energy scale. This problem can be avoided if SUSY is introduced, since the 
quadratic divergence is canceled between the fermionic and bosonic loop contributions. For 
general reviews on SUSY models, see other references (Nilles, 1984; Haber and Kane, 1985). 

The minimal SUSY extension of the SM is called the minimal supersymmetric Standard 
Model (MSSM). In MSSM, SUSY partners (which have a different spin by 1/2 from the 
corresponding ordinary particles) are introduced for each particle in the SM. For quarks and 
leptons, complex scalar fields, squark (q) and slepton (I), are introduced. The superpartner 
of the gauge boson is a gauge fermion (or a gaugino), and that of the Higgs field is called 
a higgsino(iJ). The superpartners of gluon, SU(2) and U(l) gauge bosons are a gluino 
(G), a wino (W) and a bino (B), respectively. After electroweak symmetry breaking, the 
wino, bino and higgsino mix with each other and form two charged Dirac fermions, called 
charginos (xf) (i = 1,2), and four Majorana fermions, called neutralinos (xfl) (i = 1 — 4). 
As for the Higgs sector, the SUSY models contain at least two Higgs doublet fields. They are 
required separately for the mass terms for up-type quarks, and those for down-type quarks 
and charged leptons to eliminate any possible conflict with SUSY in the Yukawa couplings. 
The particle contents of the MSSM are listed in Table [VI]. 

TABLE VI. Particle Contents in the MSSM. 



Ordinary particles 


SUSY particles 


Particle Spin 


Particle Spin 


quark (q) \ 


squark (q) 


lepton (I) \ 


slepton (I) 


gluon (G) 1 


gluino (G) \ 


W ± , Z°, 7 1 
Higgs boson (h, H, A, H±) 


chargino (xf) (i = I — 2) \ 
neutralino (xf) (i = 1 — 4) \ 



The MSSM Lagrangian consists of two parts. They are the SUSY invariant Lagrangian 
and the soft SUSY-breaking terms, as follows: 

£ = C-SUSY inv + C-SUSY breaking j 

(71) 

The MSSM Lagrangian is described in more details in Appendix B. One of the important 
features of the SUSY invariant Lagrangian is that various bosonic and fermionic interactions 
are related to each other by the requirement of SUSY invariance. For example, the gauge 
coupling constants appear not only in the covariant derivative, but also in the gaugino- 
scalar-fermion interactions and the scalar self-couplings. 

The ordinary Yukawa coupling constants are included in a scalar function called the 
superpotential (W(4>i)). The Lagrangian specified by the superpotential (C SU per P otentiai) con- 
tains a set of fermionic interactions and scalar potentials, as follows: 

where the scalar field (<pf) and the left-handed Weyl field (V'il) form a chiral multiplet of 
SUSY. The superpotential of the MSSM is given by 
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w, 



MSSM 



(yehH^Lj + {y d ) ij H l D c i Q j + (y^.H.U^ - 



(73) 



where Ef and Li represent the supermultiplets of SU{2)l lepton singlets and doublets, 
respectively. Also, Qi Uf and are the supermultiplets for quark doublets, up-type quark 
singlets, and down-type quark singlets, respectively. Hi and H 2 are two Higgs doublet fields. 
From this superpotential in Eq . (|73|) , the following Yukawa couplings are induced: 



C 



Yukawa 



(ye) ij H 1 e iR l jL + (ycdijHidmqjL + (y u ) ij H 2 u iR q jL 



H.c. 



(74) 



In addition to Eq.(|74|), the superpotential in Eq . ([731) generates the higgsino mass term, 
various Yukawa-type couplings, and the three- and four-point scalar couplings, according to 



2. Flavor problems in super symmetric models 

In the MSSM, the masses of superparticles, i.e., squarks, sleptons, and gaugino, are 
generated by the soft SUSY-breaking mass terms, which are defined as SUSY-breaking 
terms that do not induce quadratic divergence. In general, the soft SUSY-breaking mass 
terms become a new source of flavor mixing in the MSSM, which is not necessarily related 
to the flavor mixing in the Yukawa coupling constant in Eq . (|73|) . For the slepton sector, the 
soft SUSY-breaking mass terms are given by 

£soft = -(m 2 E )ije* Ri e Rj - (m 2 L )ijlu hj - { m o{ A e)ijHie Ri l Lj + H.c.}, (75) 

where (ra^j and (m^)^ are the mass matrices for the right-handed sleptons (e R ) and 
left-handed sleptons (iz), respectively; m is a SUSY-breaking parameter and A e is a di- 
mensionless scalar trilinear coupling matrix. 

In the basis where the lepton mass matrix is diagonalized, if there exist non-zero off- 
diagonal matrix elements in the slepton mass matrix, LFV is introduced. From the diagram 
in Fig.[?|, the constraint on the off-diagonal elements on the slepton mass matrix, for instance 
the jl — e element (Am^), can be placed from the present upper limit of /i + — > e + 7 decay 
as follows: 



Me < 10 



m? ~ VlOOGeW 



(77^7^ (76) 



where mj is the mass of a slepton. Similar constraints on the squark mass matrix elements 
are obtained from the flavor-changing neutral current (FCNC) processes in the quark sec- 
tor. For example, the observed K° — TC° mixing places constraints that the possible SUSY 
contribution must be small, resulting in that the squarks in the first and second generations 
must be degenerate at the level of a few % in the case that the squark mass is a few 100 
GeV and that the squark mixing angle is of similar magnitude to the Cabibbo angle. These 
constraints from LFV and FCNC suggest that there should be a special suppression mech- 
anism on the flavor mixing of sfermions (squarks and sleptons) from the dynamics of SUSY 
breaking. It is called the SUSY flavor problem. 

There are several scenarios to solve the SUSY flavor problem: 
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• The soft SUSY-breaking mass terms have a universal structure at a very high en- 
ergy scale, such as the Planck scale (~ 10 19 GeV) (Gravity mediated SUSY breaking 
scenario) (Nilles, 1984). 

• The SUSY-breaking effects are mediated by the SM gauge interaction so that squarks 
and sleptons with the same quantum numbers receive the same amount of the soft 
SUSY-breaking mass terms (Gauge mediated SUSY breaking scenario) (Dine and Nel- 
son, 1993; Dine, et al, 1995; Dine, et al, 1996; Giudice and Rattazzi, 1998). 

• There is some approximate flavor symmetry which produces nearly degenerate masses 
for squarks and sleptons, at least for the first two generations (Flavor Symmetry Sce- 
nario) (Barbieri, et al, 1996). 

• Squarks and sleptons can be diagonalized in the same basis as the quarks and lepons 
(Alignment Scenario) (Nir and Seiberg, 1993). 

• The squark and slepton masses are heavy enough (10 -100 TeV) to avoid constraints 
from FCNC and LFV, at least for the sfermions in the first two generations (Effective 
SUSY scenario) (Cohen, et al, 1996). 

The minimal supergravity model (SUGRA) is a realization of the first scenario. There 
are many phenomenological analyses based on the supergravity model. In this model, all 
squarks and sleptons receive the same magnitudes of the soft SUSY-breaking mass terms 
through the coupling of supergravity at the Planck scale, so that the mass matrices are 
diagonal with the same diagonal elements. Therefore, there is neither FCNC nor LFV at 
that energy scale. This does not necessarily imply that LFV effects should not exist in this 
scenario. In fact, if there is some interaction which breaks lepton flavor conservation between 
the Planck (« 10 19 GeV) and the electroweak scales 10 2 GeV), the LFV effect could be 
induced in the slepton mass matrices through radiative corrections (Hall, et al, 1986). In 
recent years, it has been noticed that such an effect can induce muon LFV processes at large 
branching ratios in some models of SUSY GUT (Barbieri and Hall, 1994; Barbieri, et al, 
1995a). 

3. SUSY GUT and lepton flavor violation 

In SUSY GUT (Dimopoulos and Georgi, 1981; Sakai, 1981), the SM gauge groups of 
SU(3)c, SU(2) L and U(l)y are assumed to be unified by a larger group at a high energy 
scale. In recent years, SUSY GUT has attracted much attention because the three gauge 
coupling constants determined at LEP and SLC are consistent with the SU(5) GUT pre- 
diction if contributions from SUSY particles are taken into account in the evolution of the 
coupling constants. The three coupling constants are then unified at 2x 10 16 GeV (Amaldi, et 
al, 1991; Langacker and Luo, 1991; Ellis, et al, 1991). It suggests that SUSY GUTs with 
the SU (5) group or other gauge groups which include SU (5) could be a viable candidate of 
physics beyond the SM. 

Let us first discuss how LFV would be induced in SU(5) SUSY GUT. In this model, 
quarks and leptons are classified in the three generations of 5 and 10 representations of the 
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SU (5) group, where the 5 representation (Fi) contains d\ L and In and their superpartners, 
and the 10 representation (Tj) contains q^, u\ L , and e c iL and their superpartners. The 
Yukawa coupling constants at the GUT scale are determined by the superpotential, 



W SU (5) = Ayu)ijTi ■ Tj ■ H(5) + (y d )ijFi 



T, 



H(5), 



(77) 



where H(5) and S(5) are two Higgs fields associated with the 5 and 5 representations, 
respectively. By substituting the fermionic fields for Tj and T, and also the Higgs boson 
fields for H(5) and H{5), the Yukawa couplings responsible to the quark and lepton masses 
are obtained. Here, the matrix (y u )ij corresponds to the Yukawa coupling matrix for the 
up-type quarks, and (yd)ij to that for the down-type quarks, and the leptons. 

In the minimal supergravity model, all of the scalar partners of the quarks and leptons, 
namely squarks and sleptons, are assumed to have a common SUSY-breaking mass from the 
coupling of the gravity interaction. In addition to the mass terms, the scalar triple couplings 
also have a universal structure so that they are proportional to the corresponding Yukawa 
coupling constants. At the Planck scale, the soft SUSY-breaking mass terms are given by 



= t = 



soft 



-m 2 {T/T + T, - 



moMMijTi • Tj ■ H(5) + (y d ) l3 F t ■ Tj ■ H(5)} + H.c. 



(7f 



where m is the universal scalar mass, and A is the universal trilinear coupling. At this 
stage, there is no LFV in the slepton sector. When the lepton Yukawa coupling constants 
(in this case yd) are diagonalized by unitary transformations on the bosons and fermions of 
Fi and Tj for each generation indices, the soft SUSY-breaking mass terms for the sleptons 
can also become diagonal in the same basis. This is no longer true if we take into account 
the radiative corrections to the soft SUSY-breaking mass terms due to the Yukawa coupling 
constants. In particular, since the Yukawa coupling constant corresponding to the top-quark 
mass is surprisingly large, its effects on the soft SUSY-breaking mass terms are expected to 
be sizeable. In the basis where the matrix of up-type Yukawa coupling constants is diagonal, 
all members of T, including the right-handed slepton masses, are obtained by 



m 



m% — 



m 



and 



m 2 + Am 2 



Am 2 ^-A|( 2/u)33 | 2m 2 ( 3 + | Ao | 2)ln( ^ )) 

o7T IVlQ 



(79) 
(80) 



where Mp and Mq denote the reduced Planck mass (~ 2 x 10 18 GeV) and the GUT scale 
(~ 2 x 10 16 GeV). Am 2 arises from the evolution of the renormalization group equation 
between the Planck and the GUT scales through the diagrams in Fig.^|. 

Since the physical LFV effect is induced by a mismatch of the lepton and slepton diag- 
onalization, the off-diagonal terms of the slepton mass matrix can be examined in the basis 
where the Yukawa coupling constant for leptons is diagonalized. If it is diagonalized by 



VRVeVl = diagonal, 



(81) 



the off-diagonal elements of the right handed slepton mass matrix in this new basis are given 
by 
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2 -\ ^ ° l.33|2 



(VRh(V R )* 3 \yir<(3 + \A f) \n(-f). (82) 



It becomes a source of /i + — > e + 7 decay through the diagrams in Fig.|[ If the SU (5) GUT 
relation for the down-type quark and lepton Yukawa coupling constants given by 

Ve = yL (83) 
is used, Vr is given by the CKM matrix at the GUT scale as 

(V R )ij = {VS^uh, (84) 

where Vckm can be obtained from the CKM matrix at the electroweak scale by taking into 
account the effects of running coupling constants from the electroweak to the GUT scales. 



The prediction of the branching ratio of /z + — > e + 7 decay is presented in Fig.[T0] for typical 
SUSY parameters in SU(5) SUSY GUT. The branching ratio reaches the order of 10 -14 for 
a slepton mass of a few 100 GeV/c 2 . 

Some remarks on SU(5) SUSY GUT are presented in the following. 

• In the SU (5) SUSY GUT model, LFV appears only in the right-handed slepton sector 
for small or moderate values of tan f3, which is defined by the ratio of two Higgs vacuum 
expectation values as tan/? =< H® >/< >. This is because the renormalization 
effects contribute only to en, and not to II- As a result, the helicity of an electron 
(positron) in LFV processes becomes only right-handed (left-handed). For instance, 
only /i + — > e^7 decay occurs, not /i + — > e^7. These two processes can be distinguished 
when the angular distribution of the /i + — > e + 7 signal is measured by using polarized 
muons. 

• There is partial cancellation among the Feynman diagrams which contribute to the 
ji + — > e + 7 amplitudes in the SU (5) SUSY GUT (Hisano, et al, 1997). This cancelation 
can be seen in Fig.|l0"|. 



In Eqs.(|82D and (0), the off-diagonal elements of the right-handed slepton mass matrix 
are determined by the CKM matrix elements. When the favorable values of the CKM 
matrix elements are used, |A T | = |(Vr)i3(Vr)2 3 | of (3 — 5) x 10~ 4 is obtained. This 
results from the assumption that all of the Yukawa coupling constants are generated 
from the superpotential in Eq.([T7]). However, it is known that this assumption does not 
yield a realistic mass relation for the down-type quarks and charged leptons in the first 
and second generations. If higher dimensional terms or different SU(5) representations 
of Higgs fields are included to remedy this problem, the simple relationship between 
the slepton mixing matrix and the CKM matrix in Eq. (|84D would be lost, and A T would 
become essentially a free parameter. As a consequence, the predicted branching ratios 



could be different from those in Fig.|10|. For example, if |A r | ~ 10 , the branching ratio 



is enhanced by three orders of magnitude from FigJTO. In addition, for a large tan (3, 



a further enhancement can be expected (Arkani-Hamed, et al, 1996a). One example 
prediction of the fi + — > e + 7 branching ratio for a large tan/5 is shown in Fig.|TT| 
(Hisano, et al, 1998b). When the higher dimensional terms are included, the mass 
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matrix for left-handed sleptons also has off-diagonal elements owing to the large bottom 
Yukawa coupling constant, and therefore the branching ratio of /i + — > e + 7 decay is 
enhanced by (m T /m^) 2 m 10 2 , just as in the case of SO (10) SUSY GUT discussed 
below. This enhancement can be seen in the non-minimal case in Fig.|TT|. Furthermore, 
the destructive interference between the different diagrams in the minimal SU(5) SUSY 
GUT discussed before may disappear. 

The fi + — > e + e + e~ decay and the yT— e~ conversion process receive contributions from 
the off-shell photon, Z-penguin and box diagrams, in addition to the photonic penguin 
diagram contributing to the /i + — > e + 7 decay. In SU(5) SUSY GUT, the relative 
magnitudes between different contributions vary over the SUSY parameter space. This 
implies that the relations in Eqs.(^) and (|70|) may change for different SUSY input 
values. The predictions of the fi + —>■ e + e + e~ decay and the fi~—e~ conversion in 



SU(5) SUSY GUT are shown in Figs.[l2| and [L3|, respectively. 



Large LFV effects are also expected in SO (10) SUSY GUT. In the minimal 50(10) 
SUSY GUT, the superpotential is given by 

Wso(io) = ^(Vuh^i ■ • *i + \(Vd)i^i ■ $d • (85) 

where is the 16-dimensional representation of SO(10) and $ u and $^ are two 10- 
dimensional Higgs fields. In this model, both the left-handed and right-handed sleptons 
receive LFV effects. In particular, the diagrams shown in Fig.[l4] give a large contribution 
to the amplitude of /i + — > e + j decay because it is proportional to m T . Hence, the branching 
ratio is enhanced by (m T /m M ) 2 compared to the minimal SU(5) SUSY GUT. Owing to this 
enhancement, the branching ratios for muon LFV processes can become comparable to the 
present experimental upper bounds if the slepton mixing matrices are related to the ob- 
served CKM matrix elements. The predictions for the branching ratio of y + — > e + 7 decay in 
50(10) SUSY GUT are shown in Fig.|15|. In 5O(10) SUSY GUT, since the photon-penguin 
diagrams dominate in the amplitudes of fi + — > e + e + e~ decay and e~ conversion, their 
predicted branching ratios with respect to /i + — > e + 7 decay would follow Eqs.(|69|) and ( |70|) 
over a wide SUSY parameter space (Barbieri, et ai, 1995a). 

The rates of muon LFV processes depend on the structure on the Yukawa coupling 
constants at the GUT scale. The branching ratios for the muon LFV processes are calculated 
in various realistic 5O(10) SUSY GUT models (Arkani-Hamed, et ai, 1996a; Ciafaloni, et 
ai, 1996; Duong, et al. } 1996; Gomez and Goldberg, 1996). In the 5O(10) SUSY GUT 
models, the diagrams relevant to the fi + — > e + 7 amplitude would induce the electric dipole 
moments (EDM) of the electron and neutron (Dimopoulos and Hall, 1995). The branching 
ratios for muon LFV processes are compared with the prediction of EDMs, FCNC processes 
in the quark sector, and CP violations in B and K meson decays in both the SU(5) and 
5O(10) SUSY GUT models, and the leptonic signals are shown to be very sensitive to the 
interaction at the GUT scale (Barbieri, et ai, 1995b). The LFV process is also investigated in 
the breaking pattern of 50(10) -> SU(3)xSU L (2)xSU R (2)xU B _ L (l) -> SU(3)x SU L (2)x 
[/ y (l)(Deshpande, et ai, 1996) and in the SU(A) x SU L (2) x SU R (2) model without GUT 
unification (King and Oliveira, 1999). 
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In some theoretical scenarios in which the mass matrices for squarks and sleptons at very 
high energy scale are not universal, but have some correlations or alignments with corre- 
sponding fermion matrices, it is possible to avoid the SUSY flavor problem and at the same 
time still have muon LFV branching ratios large enough to be detected. One such realization 
has been investigated concerning the dynamical alignment mechanism (Rattazzi and Sarid, 
1996). An interesting possibility is a class of models based on U(2) flavor symmetry where 
both the Yukawa coupling constants and the soft SUSY-breaking mass terms are controlled 
by the same approximate symmetry. In this case, the branching ratio for fi + —>■ e + 7 decay 
is expected to receive large SUSY contributions (Barbieri, et ai, 1996). 



4- Super •symmetric models with a right-handed neutrino 

A large LFV effect can be expected if the supermultiplets of right-handed Majorana 
neutrino are included in the SUSY Standard Model (Borzumati and Masiero, 1986). As 



explained in Section |II A 3| , the smallness of the neutrino masses can be explained by the 
see-saw mechanism. To include this see-saw mechanism, part of the lepton sector in the 
Lagrangian in Eq.(|73D is replaced by 



W N = (y e )ijHiE^Lj + (y I/ ) i jH 2 N i Lj + ^(M^yJV^, (86) 

where Ni is the right-handed neutrino supermultiplets and (Mr)^- is the Majorana mass 
matrix, and a new Yukawa coupling constant matrix, y„, is introduced. Since there are two 
Yukawa coupling matrices (y e and y u ) in the lepton sector, flavor mixing would arise and 
lepton flavor would no longer be conserved separately for each generation, just as in the 
quark sector. In SUSY models with universal soft SUSY-breaking at the Planck scale, the 
flavor mixing in left-handed sleptons would induce sizeable LFV effects in muon and tau 
decays through the renormalization effects from the Planck to the Majorana mass scales. 

The expected magnitudes of the LFV effects depend on the Yukawa coupling constant 
and the flavor mixing in the lepton sector. In the basis where the Yukawa coupling constant 
matrix for charged leptons is diagonalized, the mass matrix of light neutrinos is given by 



UjtJ 



sin^ 3 

(yvhiM^hiyJu — ;rA (87) 



where v is the Higgs vacuum expectation value, and Q is the Higgs vacuum angle for the two 
Higgs doublet. On the other hand, the off-diagonal terms of the left-handed slepton mass 
matrix induced by the renormalization effect is 

K> - -^(y-TMkjm 2 (3 + \A \ 2 ) ln(^). (88) 



In general, there is no direct relationship between the neutrino mixing in Eq.(p7[) and the 
slepton mixing relevant to /i — > ej, r — > //7 and r — > in Eq.(^). If, however, we assume 
that the neutrino mixing mostly originates from the neutrino Yukawa coupling constants, 
(y u )ij, the information from the atmospheric and solar neutrinos can be related to the slepton 
mixing. Then, the branching ratios for /i + — > e + 7 and r — > /ry decays can be evaluated by 
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using the neutrino mixing parameters (Hisano, et ai, 1995; Hisano, et al., 1996; Hisano, et 
ai, 1998a; Hisano and Nomura, 1999a). Fig. 16 shows the predicted branching ratio for the 
yU + — > e + 7 decay for different solutions of the solar neutrino mixing. As can be seen in Fig.[16|, 
it can reach the present experimental bound for the Majorana mass larger than O(10 14 ) GeV 
if the large-angle MSW solution for the solar-neutrino problem is chosen. It is noted that 
from Eq . (jS7|) , given a fixed value of the light neutrino mass (m„), the Yukawa coupling 
constant (y u ) becomes larger for a larger value of the Majorana mass scale (Mr), resulting 
in that the LFV rate becomes larger in Eq.fl8~8|). The LFV rates increases approximately as 
the second power of Mr and, therefore, it could possibly probe the mass scale of the right- 
handed Majorana neutrino (Mr) in this scenario. Note that this prediction is in contrast 
to the see-saw mechanism without SUSY, in which the LFV rates for charged leptons are 
extremely suppressed, as discussed in Section |111 C\ 



5. Other supersymmetric models 

Observable effects of muon LFV processes may arise through the renormalization effects 
in the slepton mass matrix. Because these effects may come from anywhere between the 
Planck and the electroweak scale, it is possible to consider some other interaction at a high 
energy scale as a new source of LFV. 

Such an example is given in the context of the gauge-mediated SUSY breaking. If the 
mixing between the messenger fields of gauge mediation and the ordinary matter fields is 
allowed, a sizeable LFV effect can be generated through the renormalization of the slepton 
mass matrix (Dine, et ai, 1997; Dubovsky and Gorbunov, 1998). A similar effect appears in 
a supersymmetric E e type model, where the mixing between the ordinary leptons and the 
exotic leptons in the E e 27 representation could induce a large LFV effect via the slepton 
mixing (Kitano and Yamomoto, 1999). 



6. LFV in slepton production and decay 

If the charged sleptons are discovered in experiments at future colliders (like LHC and 
LC), LFV due to slepton mixing could be directly searched in their production and decay 
processes (Krasnikov, 1994; Krasnikov, 1996; Arkani-Hamed, et al., 1996b). For example, a 
process like e + e~ — > l + l~ — > e ± Xi/ i=F Xi breaks lepton flavor conservation, where the slepton 

is assumed to decay to a lepton and the neutralino Xv Direct searches for the e/i final 
states in e + e~ and e~e~ collisions can probe the slepton-mixing angle between the selectron 
(e) and the smuon (ft) of less than 0.1, which could be comparable to indirect searches such 
as the muon LFV processes (Arkani-Hamed, et ai, 1996b). 

The production cross section for the slepton LFV processes depends on the masses and 
mixing of the sleptons. If two sleptons, like e and ft, are almost degenerate in their masses, 
a possible oscillation between them would occur in the decay process. Another interesting 
possibility is a CP violating signal in the slepton oscillation, which may arise due to a new 
complex phase in the slepton mixing matrices (Arkani-Hamed et ai, 1997; Bowser-Chao and 
Keung, 1997). 
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In SU (5) SUSY GUT, LFV in the slepton pair production and decays, in particular with 
a tau in the final state, has been investigated, and it is found to be more sensitive than 
the searches for r — > e(/j)7 (Hirouchi and Tanaka, 1998). LFV signals in the production 
of left-handed sleptons at and e + e~ colliders are also considered in the framework of 

MSSM with right-handed neutrino motivated by the atmospheric-neutrino data. It is shown 
that the decay modes of either r + \x + 4 jets + missing energy or r + fi + I + 2 jets + 
missing energy could be useful to see the signals with suppressing any potential background 
(Hisano et al, 1999b). Other possible LFV searches at a e/x collider (Choi et ai, 1998) and 
also a e7 collider (Cao et al, 1999) are also considered. 

C. Other theoretical models 

In the late 70's, especially in 1977, there were many theoretical papers on models with 
heavy neutrinos to discuss LFV (Cheng and Li, 1977a; Cheng and Li, 1977b; Bjorken, et ai, 
1977; Lee, et ai, 1977a; Lee and Shrock, 1977b; Altarelli, et al, 1977), after responding to 
the false rumor on /i + — > e + 7 signals at SIN (Schweizerisches Institut fur Nuklearforschung) . 
Since then, various theoretical models have appeared. In this section, we discuss LFV effects 
in theoretical models other than SUSY models. Although many possibilities to induce LFV 
effects might exist, we discuss only a few examples of different types of LFV. For other 
references, see, for instance, Vergados (1986), and Depommier and Leroy (1995). 

1. Models with a massive neutrino 

The simplest way to violate lepton flavor conservation is to introduce neutrino masses 
and mixing. However, it has been known that the branching ratio of /i + — > e + 7 decay from 
ordinary neutrino mixing is very suppressed if the neutrino mass and mixing suggested by the 
atmospheric and solar-neutrino experiments are used. For example, the predicted branching 
ratio from the Dirac neutrino masses and mixing is given by (Petcov, 1977; Bilenky, et ai, 
1977) 

2 

(89) 

where (Vmns)™ is the lepton flavor mixing matrix (MNS matrix) defined in Eq. (|20"D . It is 
represented by 

VLa = Ys( V MNs)aiVlA, (90) 

i 

where v^ a is the neutrino field in the weak flavor basis, and Uu in the mass-diagonalized 
basis. Even if a 1 eV neutrino mass with maximal mixing is considered, Eq. (|89|) only gives 
a branching ratio on the order of 1CU 47 . 

For the Majorana neutrino model of see-saw type (Cheng and Li, 1980), the suppression 
factor of (m^/m^) in Eq . (|89"D is replaced by a factor of 0{m v J Mr), where Mr is the mass 
of a heavy Majorana neutrino. Then, the branching ratio is still O(10 -40 ) or less for m u = 1 
eV and Mr = 10 10 GeV. It is, therefore, difficult to expect observable LFV effects from the 
ordinary neutrino masses and mixing indicated by the atmospheric and solar neutrinos. 
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2. Models with a doubly charged Higgs 

There is a wide class of theoretical models which have an interesting enhancement for 
fi + e + e + e~ decay and fi~—e~ conversion. If LFV is induced by a penguin diagram with 
a heavy boson and a light charged fermion, the photonic form factors /£o(<? 2 ) and /mo(? 2 ) 
receive an enhancement factor of In {mb OS on/ m fermi<m) from the loop diagrams in which a 
photon is attached to the internal light fermion line (Marciano and Sanda, 1977a; Marciano 
and Sanda, 1977b; Wilczek and Zee, 1977; Raidal and Santamaria, 1998). Since this factor is 
absent for fEi(q 2 ) and fMi(q 2 ), the branching ratio of /i + — > e + j decay has no enhancement. 
One of such models is a double charged Higgs boson which has the following interaction of 

C = h ij ^e jR (j) ++ + H.c, (91) 

where the loop diagrams with internal charged leptons and doubly charged scalar <p ++ could 
induce a logarithmic enhancement (Raidal and Santamaria, 1998). 



3. Super symmetric models with R-parity violation 

Another important class of models which might induce logarithmic enhancement is the 
SUSY models with R parity violation. In the MSSM, if we only require gauge invariance to 
write all possible superpotentials, the following interactions would also be allowed: 

W = \ l0k UL 3 E c k + X' ijk LiQjDl + X-^UiDpi - ^UH 2 . (92) 

These interactions violate the baryon- or lepton-number conservations. To forbid proton 
decays that are too fast, a parity which distinguishes superparticles from ordinary particles 
is often required. This parity is called the R parity, and is defined as R = (_i) 3B - L + 2S ' ) where 
B, L, and S are the baryon number, lepton number, and spin, respectively. We can, however, 
consider models with .R-parity violation when combinations of the coupling constants which 
induce proton decays are highly suppressed, but the other coupling constants can be large. 

Some combinations of the coupling constants are severely constrained by the LFV pro- 
cesses. It is known that the allowed values of AA and AA' and A'A' still give large contri- 
butions at a tree level to the fi + — > e + e + e~ and fi~—e~ conversion processes (Kim, et ai, 
1982; Huitu, et ai, 1998; Faessler, et ai, 1999). Typical relevant tree diagrams are shown 
in Fig.0. 

In SUSY models with i?-parity violation, there also exist loop contributions for muon 
LFV processes, such as fi + — > e + 7 decay (Chaichian and Huitu, 1996), /i + — > e + e + e~ decay 
and fi~—e~ conversion. The latter two processes would receive a logarithmic enhancement 
(Huitu, et ai, 1998). Fig.|18| shows typical loop diagrams. From this one- loop diagram, the 



[i — e conversion process can also be induced from the AA coupling constants. For example, 
from the loop diagrams with internal leptons, the following four form factors are given as 




2 



/Ml(0 ) = T/Jn( 0) = -|^=f, (94) 
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where r = m?/(— q 2 ). F(r) takes a value of the order of 0(1), except for r > 1, F(r) = 
lnr + |. As a result, some of the AA combinations are more severely constrained by the 
fi~—e~ conversion than the tree-level process of n + — > e + e + e~ decay. 

The present constraints for the combinations of AA and A A couplings are summarized 
in Huitu, et al. (1998) and Faessler, et al. (1999). 

4- Models with Z' 

There are many models where non-photonic LFV transitions occur at tree levels. Typical 
examples are a model with a Z' that has flavor off-diagonal couplings, or models with 
extra fermion families which mix with the three families at the tree level so that the Z 
boson has a LFV couplings (Bernabeu, et al, 1993). In such cases, the fi + — ► e + e + e~ and 
e~ conversion processes are expected to be much more important than /x + — > e + 7 decay. 



5. Models with Lorentz non-invariance 

Recently, a possible violation of the Lorentz invariance has been suggested (Colemen 
and Glashow, 1999). In this context, the Lorentz transformation is not invariant, but only 
the translational and rotational symmetries are assumed to be exact in a preferred system. 
Thus, the maximum attainable velocity could be different for each species of particles, and 
this would cause many unique phenomena in particle physics and cosmic-ray physics. 

One of the good tests of the violation of Lorentz invariance is the muon LFV processes 
(Colemen and Glashow, 1999). If a small Lorentz-non-invariant interaction exists in the 
SM Lagrangian, flavor mixing couplings can be generally allowed in the photon-fermion 
interaction. The current limit on the fi + — ► e +/ y branching ratio puts a strong constraint on 
the relevant coupling constants. Another interesting effect is a change of the muon lifetime 
at a high energy. Since the \x —>■ decay width due to the Lorentz non-invariant interaction 
would increase with 7 4 , where 7 is the Lorentz factor, it would dominate over the ordinary 
muon decay that increases with 7. Therefore, the muon lifetime might start decreasing as 
7 -3 at a sufficiently-high energy. The current limit on the energy dependence of the muon 
lifetime has been given by the experiment of the muon anomalous magnetic moment. 



D. Lepton flavor violation with Polarized Muons 

In this subsection, we discuss the usefulness of polarized muons in searches for [i + — > 
e + 7 and n + — > e + e + e~ decays. As discussed later in Section fVj highly-polarized /i + s 
(surface muons) are available experimentally. Therefore, it would be useful to examine what 
kind of new information can be obtained by measuring the angular distribution of decay 
products with respect the muon polarization. 

When the initial muon is polarized in fi + —>■ e + 7 decay, the angular distribution of the 
positron is given by 

dB ^ + T +7) = 1927r ' A *! 2 t 1 " P m cos ^) + l^l 2 i 1 + ^* «***)) . ( 95 ) 

CI COS up 
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where 6 e is the angle between the muon polarization and the positron momentum in the 
muon rest frame. Ar and Al are given in Eqs.(^) and (|60|). P M is the magnitude of the muon 
polarization. A measurement of the e + polarization would give the relative amplitudes of Ar 
and Al, which correspond to the emission of right-handed e + (fi + — > e^7) and left-handed 
e + — > e^7), respectively. This is shown schematically in Fig.|T9|. 

Since Al and Ar are model-dependent, it would be useful to discriminate between dif- 
ferent LFV mechanisms. For instance, the minimal SU(5) SUSY GUT model predicts a 
vanishing Ar and a non-zero Al-, yielding a (1 + P M cos 9 e ) distribution. On the other hand, 
the 50(10) SUSY-GUT model predicts the helicity amplitudes for both right-handed and 
left-handed e + s. For non-unified supersymmetric models with right-handed neutrino, Ar is 
non-zero but Al is vanishing, giving a (1 — P^cos9 e ) distribution. Thus, a measurement 
of the angular distribution of e + with respect to the direction of muon polarization would 
provide a valuable means to clearly discriminate between these models. 

The /i + — > e + e + e~ decay with polarized muons would provide us with an interesting 
possibility of measuring T violation (Treiman, et al, 1977; Zee, 1985). A T-odd triple 
vector correlation, <t m • (pi x p 2 ), can be defined, where is the muon spin, and pi and 
f>2 are two independent momenta of the e + in the final state. If CPT invariance holds, 
the information on CP violation in the LFV interaction may be obtained from this T-odd 
correlations of the decay products. The T-odd asymmetry would arise as an interference 
between the photon-penguin terms and the four-fermion terms. Details are discussed in 
Section [VB 1| . In particular, the T-odd asymmetry is evaluated in the SU (5) SUSY GUT 
model based on supergravity (Okada, et al, 1998). It is shown that an asymmetry up to 20% 
is possible if CP-violating phases are introduced in the soft SUSY-breaking mass terms. It 
would give independent information from the electric dipole moment (EDM) of the electron 
and neutron. 

In the fi + — > e + e + e~ decay, parity-odd asymmetries can be also defined if the initial 
muons are polarized. These asymmetries are sensitive to the chiralities of the terms in the 
effective Lagrangian: both the photon-penguin terms (Al and Ar) and the four fermion 
coupling terms (gt). Measurements of the parity-odd asymmetries in fi + — > e + e + e~ decay, 
together with the branching ratios of /i + — > e + 7 and fi + — > e + e + e~ decays, are useful to 
distinguish different SUSY GUT models (Okada, et al, 1999). 



E. |ALj| = 2 processes 

The LFV processes with |AL»| = 2 are like the muonium to anti-muonium conversion 
(Mu — Mu conversion), the wrong-flavor muon decay (/i + — > e + u fl V e ), and others. The 
phenomenology of the Mu — Mu conversion can be described by an effective four fermion 
interaction. As an example, the interaction of (V — A)(V — A) type was considered by 
Feinberg and Weinberg (Feinberg and Weinberg, 1961). It is given by 

^Mulu = (^p>7A(l - 75)e/i7 A (l - 75)e + H.c, (96) 

in which G Mu j^ is a coupling constant characterizing the strength of the Mu— Mu conversion. 
As described in Section [VE3| , the present experimental limit on G Mu j^ < 3.0 x 10~ 3 Gp is 
given (Willmann, et al., 1999). 
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In general, there could be various combinations of different types of four fermion inter- 
actions. They are such as (V + A) (V + A) , (V - A) (V + A) , (S - P) (S - P) , (S + P) (S + P) , 
(S — P)(S + P), SS and PP, where V, A, S and P are vector, axial vector, scalar, psue- 
doscalar effective interactions, respectively. The types of interactions are determined by 
assumed theoretical models. In Fig.^0], example diagrams of speculative theoretical models 
responsible for the Mu — Mu conversion are shown. They are models of (a) the exchange of 
doubly charged Higgs bosons, (b) that of heavy Majorana neutrinos, (c) that of a neutral 
scalar particle (including a superparticle), and (d) that of a bilepton gauge boson. They are 
briefly described in the following. 

A simple example which induces the |ALj| = 2 process is provided by a model with 
a doubly- charged singlet scalar boson (Chang and Keung, 1989). In this model, only the 
(V A +A)(V+A) interaction is generated and the Mu— Mu conversion rate can be as large as the 
present experimental limit within the constraints from the measurements of the anomalous 
muon magnetic moment and the high-energy Bhabha scattering. A more general case with 
the doubly charged scalar boson is also considered (Swartz, 1989). 

In the left-right symmetric model with a triplet Higgs boson field, the Mu— Mu conversion 
could be induced by a doubly charged Higgs boson, A ++ (Halprin, 1982). In this model, if 
the mass of the muon neutrino is greater than 35 keV/c 2 and less than the present direct 
experimental bound of 170 keV/c 2 , there is an upper limit of the neutrino lifetime in order 
for the neutrino energy density in the universe not to exceed its present total energy density. 
From this requirement, a lower bound of G Mu mH is predicted as a function of the mass of 
the muon neutrino, m„ (Herczeg and Mohapatra, 1992). For the range 35 keV/c 2 < m v < 
170 keV/c 2 , a lower limit of G Mum > (1 - 40) x 10~ 4 G F is obtained. 

Supposing that neutrinos are of Majorana nature, the Mu — Mu conversion could take 
place by an intermediate pair of neutrinos. This coupling is related to neutrinoless double 
(3 decays, yielding G Mu j^ < lQ~ b G F (Halprin, 1982). 

There is a class of models with a neutral scalar boson, which has a flavor-changing 
coupling to introduce the | ALj| = 2 processes (Hou and Wong, 1996). An important example 
of this kind of models is the supersymmetric model with R— parity violation, where a tau 
sneutrino exchange induces the Mu— Mu conversion (Mohapatra, 1992; Halprin and Masiero, 
1993). In this case, the four fermion coupling is of the (S — P)(S + P) type. The present 
experimental limit of the Mu — Mu conversion gives a constraint on the relevant coupling 
constant, | A 132 A231 1 < 3 x 10" 4 , for a superpartner mass of order 100 GeV/c 2 . Also, the four 
fermion coupling constant for the /i + — > e + v ll i> e decay is predicted to be similar in magnitude 
as that for the Mu — Mu conversion. 

In some extensions of the SM gauge groups, there appear doubly charged gauge bosons 
(called bileptons), X , which couple only to leptons. They occur, for instance, in SU(15) 
GUT models (Frampton and Lee, 1990), or in a gauge model with SU(3)c x SU(3)l x 
U(1)y (311 model) (Frampton, 1992). In these models, singly-charged and doubly-charged 
bilepton gauge bosons appear from breaking of the SU(3)l gauge symmetry to SU{2) L of 
the SM gauge groups. The mass bound of the bilepton gauge bosons is obtained from a 
precise determination of the Michel parameters of the normal muon decay (Carlson and 
Frampton, 1992), the muoniun hyperfine splitting, and the ji + — > e + u f J7 e decay (Fujii, et 
ai, 1994). They give a lower bound of roughly 200 GeV/c 2 . In these models, the Mu — 
Mu conversion could occur by the exchange of a doubly charged bilepton (Horikawa and 
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Sasaki, 1996; Fujii, et al., 1993). The effective interaction is of the form (V — A)(V + A), 
which is in contrast to either the traditional (V — A){V — A) interaction by Feinberg and 
Weinberg. The bilepton interaction is given by 

!IM X--l-f lb Cl T + H.c, (97) 



2V2 

where / = e, /x, r, and C is the charge-conjugation matrix. The gauge coupling constant, g^, 
is on the order of 0(1), and is determined from theoretical models. The Mu — Mu conversion 
rate (P Mum ) is given by 

^MuMH = 4.5xl0 3 x(^) 4 , (98) 

where mx is the mass of the bilepton gauge boson. From the current experimental bound 
of the Mu — Mu conversion, a constraint of m x / gzi > 2.6 TeV/c 2 is obtained (Willmann, et 
a/., 1999). 

To study the |ALj| = 2 processes, the reaction e~e~ — > fi~fi~ has been known to be 
useful for a long time (Glashow, 1961). Such an experiments were carried out to place a 
constraint on G Mu ^ (Barber, et al, 1969). Recently, similar scattering processes at a high 
energy, such as at a ee linear collider or a /x/x collider, have been discussed, based on mostly 
bilepton models (Frampton, 1992b; Hou, 1996b; Raidal, 1998b). 
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IV. NORMAL MUON DECAY 



Normal muon decay remains the only purely leptonic process of weak interaction accessi- 
ble to precise measurements with high statistics. The studies are free from the complications 
of the strong interaction and hadronic structure. For this reason, precise studies of normal 
muon decay would provide information unambiguously interpreted. 

The experimental progress has benefited from the high-intensity muon beams available at 
the three meson factories, such as the Paul Scherrer Institute (PSI) in Switzerland, TRIUMF 
in Canada, LAMPF in U.S.A. (which was unfortunately shut down), and also lately the 
Rutherford Appleton Laboratory (RAL) in England. In particular, useful is a surface fi + 
beam, which arises from the decay of pions stopped at the surface of a pion production 
target. It has a monochromatic kinetic energy of 4.1 MeV (29.8 MeV/c in momentum) and 
100% muon polarization antiparallel to its momentum direction. This high polarization of 
the surface muons is useful for various measurements requiring muon polarization. 

In the past, studies of the normal muon decay have greatly contributed to develop the 
SM. It is now playing a role to probe for possible deviations from the SM. For example, see 
reviews by Fetscher and Gerber (1995) and Herczeg (1995). In the following, the current 
status of the studies of normal muon decay, in particular those which are being prepared or 
planned for the future, is discussed. They are such as measurements of the muon lifetime, the 
Michel spectrum, and the longitudinal polarization of e + s in polarized /i + — > e + u e U fl decay. 
The other important muon experiments, such as the muon anomalous magnetic moment 
and the muon electric dipole moment, are not discussed. 



A. Muon lifetime 

1. Phenomenology 

The Fermi coupling constant, Gf, is one of the three precisely measured inputs of the 
SM, along with the fine structure constant (a) and the Z bosom mass (mz)- Their updated 
values are given in Table [VTJ. Note that the mz value in the Particle Data Group (Particle 



Data Group, 1998) is mz = 91.187 ± 0.007 (77 ppm). However, after a recent improvement 
of mz (LEP and SLD Electroweak Working Group, 1999), the uncertainties of arid mz 
become comparable in the order, as can be seen in Table |V11| . Gp is determined from the 
muon lifetime (r M ) from Eq. (|23|) . The complete two-loop QED corrections to an estimation 
of the muon lifetime was calculated (Van Ritbergen and Stuart, 1999). The theoretical errors 
to derive Gp from the muon lifetime are now reduced to be negligible, compared with the 
experimental uncertainty in the measurement of the muon lifetime. It would be necessary 
to improve the accuracy of Gf, if experimentally possible. To test the SM, it is required 
to compare Gf from the muon lifetime with those determined from other measurements, 
such as tau leptonic decays, or the My/ and other observables at the Mz pole with similar 
accuracies (Marciano, 1999). 
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TABLE VII. Three fundamental values in the Standard Model. 



Parameter 


Experimental value 


Uncertainty 


1/a 


137.0359895 ± 0.0000061 


0.045 ppm 


G F 


(1.16639 ±0.00001) x 10~ 5 GeV~ 2 


9 ppm 


m z 


91.1867 ±0.0021 GeV/c 2 


23 ppm 



2. Experimental status 

Experimentally, measurements of the muon lifetime were carried out at TRIUMF (Gio- 
vanetti, et al., 1984) and at Saclay (Bardin, et al., 1984). Since then, for more than a 
decade, no experimental efforts to improve the situation have been made. The present value 
of r M is = 2.19703(4) x 10" 6 (±18 ppm) (Particle Data Group, 1998). Recently, however, 
an experimental proposal at Rutherford Appleton Laboratory (RAL) to reduce the uncer- 
tainty of Gp by an order of magnitude has been undertaken (Nakamura, et al., 1998), and 
two experimental proposals to PSI (Cavallo, et al, 1999; Carey, et al, 1999b) have been 
submitted. 

An on-going experiment at the RIKEN-RAL muon facility, R77, uses a pulsed muon 
beam. In previous experiments with a continuous muon beam, only one incident muon within 
the time window of measurement was allowed in order to avoid any possible confusion in 
the decay of other muons occuring in the same time window. It, otherwise, would introduce 
distortion of the time spectrum of muon decay. This requirement, however, would lead 
to a limitation on muon beam intensity, and the sensitivity would be statistically limited. 
To overcome this problem, R77 at RIKEN-RAL has adopted a pulsed beam, whose pulse 
interval of 20 msec is much longer than the measurement window. Since all muons come 
at the same time, multiple muon decays within the measurement time window are allowed. 
There are several sources of systematic errors. One of them is from counting losses due to 
pile-up e + events and the dead time of detection, since the instantaneous beam intensity is 
strong. To avoid this effect, segmentation of the detectors is required. In R77, multiwire 
proportional chambers are used. Another systematic error might come from a precession of 
muon spin under an earth field. A special magnetic material will be adopted for the muon 
stopping target to depolarize the muon spin polarization. It aims, by accumulating 10 11 
muon decays, to achieve statistical and systematic errors of about 3 ppm (a total of 4 ppm) 
in the r M measurement. Its initial phase, planning to accumulate 10 10 muon decays, have 
already started in 1999. 

There are two new planned experiments at PSI: R-99-06 (Cavallo, et al, 1999) and R- 
99-07 (Carey, et al., 1999b), where R-99-07 has just been approved. The both experiments 
are aiming at a factor 20 improvement over the current world average of r M . 

The approved experiment, R-99-07, is to use a chopped surface-muon beam at PSI. 
The muons are stopped in a sulfur target to reduce the residual muon polarization, and a 
transverse magnetic field of 75 G is applied to further de-phase it. The e + detector (fiLan 
detector) consists of 180 triangular scintillating tiles distributed within the 20 super-triangles 
of an icosahedral geometry centered on the target. Each scintillating tile is viewed by a 
photo-multiplier, and the signal is recorded by a waveform digitizer developed in the g — 2 
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experiment at BNL. This Air geometry of the pLan detector, especially a sum of the point- 
like symmetric tile pairs with respect to the center, would reduce further any asymmetries 
due to spin rotation. The experiment has just been approved and a physics run is expected 
in three years. 

The experiment, R-99-06, uses a fiber-active-scintillator target (FAST), which comprises 
plastic scintillating fibers of 0.5-mm diameter in an overall active volume of 20 x 20 x 20 cm 3 . 
A 7r + beam is stopped in the FAST target, and a 7r + — > p + — > e + decay chain is observed 
in the target. Owing to a high granularity and fast response of the FAST detector, many 
individual muon decays can be recorded in parallel in a high event rate such as 1 MHz. Since 
the p + s produced from tt + decays at rest have an isotropic distribution, possible errors from 
the muon polarization can be eliminated. In its first stage, a DC tc + beam is planned to use, 
and in the second stage, a pulsed n + beam will be used. 

B. Michel decay spectrum 

1. Phenomenology 

The spectrum of the normal muon decay is given in Eq.(|3(]) with the four Michel 
parameters of p, £, 5, and rj. It includes all possible Lorentz-invariant interactions. A 
precise determination of the Michel parameters would serve stringent testing of the (V — A) 
structure of electroweak interactions in the SM, and obtain a hint of new physics beyond 
the SM. 

2. Experimental status 

In the past, each of the Michel parameters in the normal muon decay was determined by 
several different experiments: for instance, the p parameter (Derenzo, 1969), the r\ parameter 
(Burkard, et al, 1985b), the 5 parameter (Balke, et al., 1988), P M £ (Beltrami, et al., 1987), 
and P^{6/p) (Jodidio, et al, 1986). 

A new experiment, E614 at TRIUMF, to measure the entire differential spectrum of 
positrons from the decay of polarized muons is being prepared (Abegg, et al., 1996). By 
accumulating 10 9 muon decays, the goal of the E614 experiment is to measure the Michel 
parameters (p, £, 5 and rj) to precisions of Ap < 1 x 10~ 4 , A6 < 3 x 10~ 4 and A(P M £) < 
2 x 10~ 4 . The aimed precisions are 3 to 10 times better than those previously achieved. Note 
that only the product P M £ can be experimentally determined. For example, the expected 
constraint on the left-right symmetric model from E614 is presented in Fig.[| 

A schematic view of the E614 detector is shown in Fig.^T[ The E614 spectrometer 
consists of a superconducting solenoid magnet with tracking chambers. In E614, a surface 
muon beam from the M13 beam channel at TRIUMF will be brought into the detector, and 
stopped in a muon-stopping target made of aluminum located at the center of the apparatus. 
An array of planar chambers, mounted symmetrically upstream and downstream from the 
target, will track the trajectories of e + s from muon decays under a magnetic field of 2 T with 
a homogeneity of better than 10~ 4 . The positron spectrum will be measured over a wide 
range of 0.4 < x < 1.0 for the positron energy, and 10° < 9 e < 70° and 110° < 9 e < 170° for 
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the angle between the muon spin direction and the positron momentum vector (9 e ), where 



x is defined in Section II C 1 



Experimentally, it is important to keep the muon polarization fully aligned with the 
magnetic field direction at the stopping target. Any reduction of P*, the projection of P^ on 
the field direction, must be minimized at a level of 10 -4 . Possible sources of reduction are 
a misalignment between the muon-beam axis and the magnetic field at the spectrometer, 
a fringing field of the spectrometer, a contamination of non-surface muons (such as cloud 
muons), multiple scattering of muons in the production target, and possible depolarization 
of the muon spin in the aluminum muon-stopping target. At the same time, as spectrometer 
requirements, low-mass tracking chambers to minimize multiple scattering of positrons and 
position accuracy of the detector assembly are crucial. Detector construction has been 
started, and physics data-taking is expected from 2001. 



C. Polarization of e + in fi + — > e + v e V^ L decay 
1. Phenomenology 

The longitudinal polarization of e + (Pl) in pA — > e + h , e j7 fJi decay is already given in 
Eq.(p8p. When the muon is not polarized (P M = 0) and the SM values ofp = 5 = 3/4 are 
taken, P L leads to 

P L (x,cos6 e )=C (99) 

It is independent of the values of x and cos# e . Therefore, a measurement of the longitudinal 
polarization of e + s emitted by unpolarized muons would provide good a direct determination 
of the parameter £ . 

If the muon is polarized with the SM values of the p and 5 parameters, Pl is given by 

* / £P„cos^ e (2s - 1) (£"-££') , N 

PAx, cosfl e = £ + . — — - ■ - — . 100 

V ' ' ^ (3-2x)+£P tl cos9 e (2x-l) £ V ; 

From this, the measurement of Pl as a function of energy (x) and angle (9 e ) would give other 
information about the combination of parameters, (£ —££)/£• m particular, for x ~ 1 and 
cos^g ~ —1, it leads to 

P L (X = 1, COs9 e = -1) « £' + — ^ • ^ (101) 

The combination (£ "—££ )/£ is multiplied by an enhancement factor of £-P M / (1 — £P M ), which 
could be large when P M is close to unity. 

The two transverse polarization components (Pti and P^) of the e + in pA — > e + i> e T> fl 
decay are given in Eqs. (|36D and (|37|) . If the time-reversal invariance holds, Pp2, which is 



the transverse e + polarization normal to the decay plane determined by P M and the e + 
momentum direction, should be vanishing. A non-zero P?2 would signal a violation of the 
time-reversal invariance. The electromagnetic final-state interaction, which mimics a T- 
odd effect, is known to be small. On the other hand, Pti, which is not forbidden by the 
fundamental symmetries, is sensitive to the Michel parameter, rj. It would be advantageous, 
since the determination of rj from the isotropic part of the e ± energy spectrum is difficult 
owing to the small Xq factor (~ 10~ 2 ) multiplied by r\. 
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2. Experimental status 



The latest measurement of Pl of e + in /z + — > e + v e V^ decay was carried out at SIN 
(Burkard, et ai, 1985a). A magnetized- iron foil was adopted as a polarimeter for the e + 
polarization, and was tilted by 45° with respect to the e + momentum direction. Either 
Bhabha scattering of e + off e~s (e + e~ — > e + e~) or annihilation in flight (e + e~ — > 77) in 
the magnetized iron foil was utilized, where the cross sections of those processes have a 
particular dependence on the relative angle between the e ± polarization directions. Since 
the magnitude and direction of e~ polarization in the magnetized-iron foil is known, the 
polarization of e + from fi + — > e + v e V ll decay can be determined. Either e + e~ or 77 pairs 
from the polarimeter were detected by four Nal(Tl) crystal detectors located behind. Both 
unpolarized and polarized muons were studied. From the case of unpolarized muons, Pl{= 
£ ) = 0.998 ± 0.042 was obtained, whereas the current average value in the Particle Data 
Group is Pl = 1.00 ± 0.04 (Particle Data Group, 1998), showing no strong evidence of 
the right-handed current. For the case of polarized muons, they obtained (£ — £ ■ £ )/£ = 
—0.35 ± 0.33 which, however, did not make better constraints of the coupling constants at 
this moment. 

The transverse polarization components, Pti and Pt2, of e + in /x + — > e + h , e T> fl decay were 
measured by the same group (Burkard, et ai, 1985b). A magnetized-iron foil was used again 
as a polarimeter, but was placed perpendicular to the e + momentum direction. The measure- 
ment was based on the fact that two photons from the annihilation of transversely-polarized 
e + with e"s in a magnetized-iron foil are preferentially emitted in the plane determined by 
the bisector of the e + polarization vector (Pt) and the e~ polarization vector in a magne- 
tized foil. Their results were < P Tl >= 0.016 ± 0.023 and < P T2 >= 0.007 ± 0.023. From 
the measured value of Pti, r] = —0.007 ± 0.013 was obtained. 

A new experiment, R-94-10 at PSI (Barnett, et ai, 1994), is in preparation to measure the 
transverse polarization, both P Tl and Pt 2 , of the positrons from polarized /i + — > e + v e V^ de- 
cay, with a precision of 3 x 10~ 3 . The experimental principle is the same as the previous one, 
mentioned before. Major improvements are expected to be from a higher muon- stopping 
rate due to a higher proton current at the PSI cyclotron, installation of two analyzing foils 
with additional wire chamber in between, and replacement of four Nal(Tl) crystals by 127 
hexagonal BGO crystal. An engineering run will start in 1999. 

Another new experiment, R-97-06 at PSI (Van Hove, et ai, 1997), is under development 
to measure the longitudinal polarization, Pl, of the positrons emitted antiparallel to the 
muon spin from polarized /i + — > e + u e V IJL decay. As explained in Eq. ( |101|) , Pl at x ~ 1 and 
cos# e ~ —1 is sensitive to the combination of (£ — )/£ with the enhancement factor. 
R-97-06 aims to measure this observable with an improvement of more than one order of 
magnitude over the previous experiments at SIN. It will use three solenoidal magnets to 
track e + s from — > e + u e V^ decay with double-sided Si strip detectors. Two magnetized 
iron foils with opposite sign of the analyzing power are used as a polarimeter, followed by 127 
BGO crystals to detect e ± s and also photons. The asymmetry for two different analyzing 
foils are compared for the two cases of polarized and unpolarized muons, giving a relative 
measurement to reduce systematic errors. The goal is to measure 
%. An engineer run with the complete setup is planed for late 1999. 
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V. LEPTON FLAVOR VIOLATING MUON DECAYS 



The muon system is one of the best places to search for LFV, compared with the others. 



In Table |VHIi the upper limits of various lepton flavor violating decays are listed. The 
sensitivity to LFV is superb in the muon system. It is mostly because of a large number of 
muons (of about 10 14 — 10 15 /year) available for experimental searches today. The theoretical 
frameworks for LFV were already presented in Section |TJ. In this section, phenomenology 
and experimental results, including the prospect for future improvements for each of the 
forbidden muon LFV processes, are reviewed. They are such as /u + — > e + 7 decay, /i + — > 
e + e + e~ decay, fi~—e~ conversion in a muonic atom, \x~ — e + conversion, and muonium to 
anti-muonium conversion. In the first three processes, lepton flavors change by one unit 
(|ALi| = l). 

TABLE VIII. Limits of the lepton-fiavor violating decays of muon, tau, pion, kaon and Z boson. 



Reaction 


Present limit 


Reference 


/i + — ► e + 7 


< 1.2 x 10" 11 


Brooks, et al. (1999) 


/i + — ► e + e + e~ 


< 1.0 x 1(T 12 


Bellgardt, et al. (1988) 


fj,~Ti — * e~Ti 


< 6.1 x i(r 13 


Wintz (1998) 


fi + e~ — ► /j,~e + 


< 8.3 x 10- 11 


Willmann, et al. (1999) 


t — ► B'y 


< 2.7 x 10~ 6 


Edwards, et al. (1997) 


T > M7 


< 3.0 x 10~ 6 


Edwards, et al. (1997) 


T — ► flflfl 


< 1.9 x 10~ 6 


Bliss, et al. (1998) 


t — ► eee 


< 2.9 x 10" 6 


Bliss, et al. (1998) 


7T° — > \xe 


< 8.6 x 10- 9 


Krolak, et al. (1994) 




< 4.7 x 10~ 12 


Ambrose, et al. (1998) 


K + — * 7r + fj, + e~ 


< 2.1 x 10~ 10 


Lee, et al. (1990) 


Kl -> 7rV + e" 


< 3.1 x 10" 9 


Arisaka, et al. (1998) 


Z° -» ne 


< 1.7 x 10" 6 


Akers, et al. (1995) 


Z°^re 


< 9.8 x 10" 6 


Akers, et al. (1995) 


Z°^th 


< 1.2 x 10" 5 


Abreu, et al. (1997) 



A. n + — * e + 7 decay 



1. Phenomenology of fi + — ► e + 7 decay 



The most popular process of lepton-flavor-violating muon decay is fi + — > e + 7. The 
Lagrangian for the /x + — > e + 7 amplitude is given by, as shown in Eq. 



C 



AG F 
V2 



m tt A R n R (T IJ,v e L F ia/ + m^A^a^ 1 'e R F^ v + H.c. 



(102) 



The differential angular distribution of /i + — > e + 7 decay is already given, as in Eq. fl95|) , by 
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dB (^ -> e + 7 ) 

^7 a~\ = 1927r 



\A R \ 2 (1 - P M cos# e ) + | A L | 2 (1 + P^ cos9 e ) 



(103) 



where # e is the angle between the muon polarization and the e + momentum vectors. P^ is 
the magnitude of the muon spin polarization. The branching ratio is presented by 



B{n+ - e + 7 ) 



r( M + - e+ 7 ) 



r(/i+ — > e + uv) 



l ^- = ZMn\\A R \ I + \A L \ 



(104) 



5. Event signature and backgrounds 

The event signature of fi + — > e + 7 decay at rest is an e + and a photon in coincidence, 
moving collinearly back-to-back with their energies equal to half of the muon mass (m M /2 = 
52.8 MeV). The searches in the past were carried out by using positive muon decay at 
rest to fully utilize its kinematics. A negative muon cannot be used, since it is captured 
by a nucleus when it is stopped in a material. There are two major backgrounds to a 
search for fi + —>■ e +r y. One is a physics (prompt) background from radiative muon decay, 
fi + — > e + i> e n ll '-y , when e + and photon are emitted back-to-back with the two neutrinos 
carrying off little energy. The other background is an accidental coincidence of an e + in a 
normal muon decay, fi + — > e + u e V fl , accompanied by a high energy photon. The sources 
of the latter might be either fi + — > e + u e V fl r y decay, or annihilation- in- flight or external 
bremsstrahlung of e + s from normal muon decay. These backgrounds are described in more 
detail in the following. 



3. Physics background 

One of the major physics backgrounds to the search for fi + — > e + 7 decay is radiative 
muon decay, /i + — > e + v e V^ (branching ratio = 1.4 % for E 1 > 10 MeV), when the e + 
and photon are emitted back-to-back with two neutrinos carrying off little energy. The 
differential decay width of this radiative muon decay was calculated as a function of the 
e + energy (E e ) and the photon energy (E 7 ) normalized to their maximum energies, namely 
x = 2E e /m fJ , and y = 2E 1 /m ll (Fronsdal and Uberall, 1959; Eckstein and Pratt, 1959). 
The ranges of x and y are already shown in Eq.(p3|). As a background to fi + — > e + 7, the 
kinematic case when x ~ 1 and y m 1 is important. In an approximation of the limit of 
x ~ 1 and y ~ 1 with an angle between e + and photon (0 e7 ) of almost 180°, the differential 
decay width of ji + — > e + u e V^y decay is given by (Kuno and Okada, 1996) 

1 -x) 2 (l -P^cos^) + (4(1 -x)(l -y) - 2 z2 )( x + P ^ cos 9 e )\dxdyzdzd(cos 6 e ), (105) 

where 9 e is the angle between the muon spin and the e + momentum direction. Gp is the 
Fermi coupling constant, a is the fine-structure constant, z = it — 9 n , and cos z is expanded 
in a polynomial of z, since z is small. In Eq. ( |105| ), only the terms of up to the second 
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order of a combination of (1 — x), (1 — y) and z are kept. At x ~ 1 and y « 1, the 
effect of the positron mass is found to be very small, on the order of (m e /m^) 2 , and is 
therefore neglected in Eq.( |105| ). The first term in Eq.( |105| ) represents the e + being emitted 
preferentially opposite to the muon spin direction, whereas in the second term the e + is 
emitted along the muon-spin direction. When x — 1 and y = 1 exactly, this differential 
decay width vanishes. However, in a real experiment, finite detector resolutions introduce 
background events which would ultimately limit the sensitivity of a search for fi + — ► e +r y. 

Given the detector resolution, the sensitivity limitation from this physics background 
can be estimated by integrating the differential decay width over the kinematic signal box. 
It is given by 



, , _ , 1 f 1 , f 1 , /•min(fc,2 v /(l- a: )(l-y)) dT(ll + -> e+UU-f) 

dB{ji — > e uuj) = _ ax ay dz 



r(//+ — > e + vu) Ji-Sx Ji-Sy Jo dxdydz 



a 



16n 



Ji(l-P M cos# e ) + J 2 (l + P M cos# e ) d(cos# e ), (106) 



where Sx, 5y and 5z are a half width of the /i + — > e + 7 signal region for x, y and z, 
respectively. — > e + vV) is the total muon decay width. Ji and J 2 are given as the sixth 
power of a combination of 5x and fo/. For the case of 5z > 2y / 5x5y, they are presented by 

J l = (5x)\5y) 2 and J 2 = ^(5xf(5yf . (107) 

When the angular resolution meets 5z < 2y/5x5y, they are given by 

J 2 = S(Sxf(Syf(^f - 8(fa)(%)(|)" + ?(|)«. (109) 

Experimentally, the resolution of the e + energy is better than that of the photon energy, i. e. 
5x < 5y. Also, the angular resolution, 5z, has been poor in past experiments. Thereby, J2 
is much larger than J\ for most cases. 

Fig.^ shows a fraction of the fi + — > e + z/ e 77 M 7 decay for the given 6x and 5y values with 
unpolarized muons in the case of 5z > 2y / 5x5y. From Fig.^, it can be seen that both Sx 
and Sy on the order of 0.01 are needed to achieve a sensitivity limit at the level of 10~ 15 . 

Radiative corrections to the radiative muon decay for the case of the physics background 
to ji + — > e + 7 decay have been calculated to be on the order of several %, depending on the 
detector resolution (Arbuzov, et ai, 1998). 



4- Accidental background 

With a very high rate of incident muons, the accidental background becomes more im- 
portant than the physics background. This is usually the case for the present and future 
experiments. The event rate of the accidental background normalized to the total decay 
rate (-B acc ) can be estimated by 
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5acc = R P ■ / e ° " / 7 ° " (At e7 ) • (^), (110) 

where i? M is the instantaneous muon intensity. /° and /° are the integrated fractions of 
the spectrum of e + in the normal muon decay and that of photon (such as from /i + — > 
e^VeVi/y decay) within the signal region, respectively. They include their corresponding 
branching ratios. At e7 and Au n are, respectively, the full widths of the signal regions for 
timing coincidence and angular constraint of the back-to-back kinematics. 

Given the sizes of the signal region, B acc can be evaluated. Let us take Sx, Sy, 59 ey , 
and 5t ei to be the half width of the signal region for e + , photon energies, angle 9 ei , and 
relative timing between e + and photon, respectively. /° can be estimated by integrating 
the Michel spectrum of the normal muon decay over 1 — 8x < x < 1, yielding /° « 2(&e). 
Given the angular resolution, 59 ey , the back-to-back resolution (Au; e7 /47r) is presented by 
(Ao; e7 /47r) = (59 n ) 2 /A. As for /?, if the radiative muon decay /x + — >■ e + u e V^ is considered 
as a source of the 52.8 MeV photon, it can be given by integrating Eq.(|54"D over 2ir for # 7 , 
and then over the photon energy within the width of the signal region (1 — 5y < y < 1). For 
unpolarized muons, it is given by 



111^ 



From Eq.( |lll| ), it is shown that /° for /i + — > decay is roughly proportional to 

The other sources of high-energy photons are annihilation in flight of e + s in the normal 
muon decay and external bremsstrahlung. The contribution from annihilation of e + in 
flight depends on the materials along the e + 's track path. Fig.E3| shows, for instance, the 



contribution of annihilation in flight for that case of e + s passing through a muon-stopping 
target of 50 mg in thickness. It indicates that its contribution from the target is smaller 
than the radiative muon decay, and only becomes important if the photon energy resolution 
becomes extremely good. However, it is dependent on the total amount of materials in an 
experimental setup. 

From the above, the effective branching ratio of accidental background is given by 



a 



B acc = R p ■ (25x) ■ —(5y) 2 (\n(5y) + 7.33) X — • (26t). (112) 



2tt 



59' 



V 4 



For instance, take some reference numbers such as the e + energy resolution of 1% (FWHM), 
the photon energy resolution of 6% (FWHM), Au; e7 = 3 x 10~ 4 steradian, At n = 1 nsec, 
and = 3 x 10 8 ;U + /sec, 5 acc is 3 x 10~ 13 . The accidental background becomes severe. 
Therefore, it is critical to make significant improvements in the detector resolution in order 
to reduce the accidental background. 



5. Muon polarization 

The use of polarized muons has been found to be useful to suppress backgrounds 
for fi + — > e + 7 search (Kuno and Okada, 1996; Kuno, et ai, 1997a). For the physical 
(prompt) background, as already discussed in Subsection |V A 3| , the coefficient of J 2 is 
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much larger than Ji, since the resolution of the photon energy is much worse than that 
of the e + energy detection. Therefore, the angular distribution of the physics background 
follows approximately (l+P^cos^) as long as Sy > Sx. Fig.|23| shows the angular distri- 



bution of — > e + i/ e 77 M 7 with, for instance, Sy/Sx = 4. If we selectively measure the 
e + s in fi + — ► e + 7 going opposite to the muon-polarization direction, the background from 
/i + — > e + z/ e z/ A1 7 would be significantly reduced for the n + —>■ e^7 search. Furthermore, by 
varying 5x and Sy, the angular distribution of the fi + — > e + u e V II '-f background can change 
according to Eq. (|106| ), thus providing another means to discriminate the signal from the 



backgrounds. 

Regarding the accidental background, the use of polarized muons has also provided a 
means for its suppression (Kuno, et al, 1997a). It happens that the sources of accidental 
backgrounds have a specific angular distribution when a muon is polarized. For instance, 
the e + s in the normal Michel /i + decay are emitted preferentially along the muon spin 
direction, following a (l + P^cos^e) distribution, as in Eq. (|33|) , whereas the inclusive angular 
distribution of a high-energy photon (e.g. > 50 MeV) from fi + — > e + i/ e V IJt r y decay follows 
a (1 + Pfj,cos9 7 ) distribution, as in Eq. (|54l) , where 9 1 is the angle of the photon direction 
with respect to the muon spin direction. It should be noted that this inclusive angular 
distribution was obtained after integrating the energy and direction of the e + s, which is in 
contrast to the case of the physics background, where only the extreme kinematics of the e + 
and photon being back-to-back in fi + — > e^vj?^ decay is relevant. It is further noted that 
the other sources of high-energy photons, such as external bremsstrahlung and annihilation 
in flight of e + s from the normal muon decay, also follow a (1 + P M cos# 7 ) distribution. 

This inclusive angular distribution of a high-energy photon in /i + — > e + ^ e 17 M 7 implies that 
the accidental background could be suppressed for ji + — > e^7, where high-energy photons 
must be detected at the opposite direction to the muon polarization. A similar suppression 
mechanism of accidental background can be seen for /z + — > e^7 when high-energy positrons 
are detected at the opposite direction to the muon polarization. As a result, the selective 
measurements of either e + s or photons antiparallel to the muon spin direction would give the 
same accidental background suppression for ji + — > e^7 and /x + — > e^7 decays, respectively. 
This favorable situation comes from the fact that the inclusive distributions of both high- 
energy e + s and photons, respectively, in the normal and radiative muon decays, follow a 
(1 + P^cos^) distribution, where 9 is either 9 e or 6* 7 . The suppression factor, r], is calculated 
for polarized muons by 

r]= f 1 d(cos#)(l + P M cos#)(l-P M cos#)/ f 1 d(cos9) 

J cos dp J cos 9d 

= (1 - Pi) + ip M 2 (l - cosfl D )(2 + cos0 D ), (H3) 

where 9d is a half opening angle of detection with respect to the muon polarization direction. 
i] is shown in Fig.^5] as a function of 9d- For instance, for 9d = 300 mrad, an accidental 
background can be suppressed down to about 1/20 (1/10) when P^ is 100 (97)%. 
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TABLE IX. Historical Progress of search for [i + — > e + 7 since the era of meson factories with 90 
% C.L. upper limits. The resolutions quoted are given as a full width at half maximum (FWHM). 



Place 


Year 


AE e 


/X-Evy At e ^ 




Upper limit 


References 


TRIUMF 


1977 


10% 


8.7% 6.7ns 




< 3.6 x 10" 9 


Depommier, et al. (1977) 


SIN 


1980 


8.7% 


9.3% 1.4ns 




< 1.0 x 10~ 9 


Van der Schaaf, et al. (1980) 


LANL 


1982 


8.8% 


8% 1.9ns 


37mrad 


< 1.7 x 10~ 10 


Kinnison, et al. (1982) 


LANL 


1988 


8% 


8% 1.8ns 


87mrad 


< 4.9 x lO -11 


Bolton, et al. (1988) 


LANL 


1999 


1.2%* 


4.5%* 1.6ns 


15mrad 


< 1.2 x lO" 11 


Brooks, et al. (1999) 



* shows an average of the numbers given in Brooks, et al. (1999). 



6. Experimental status of /i + — > e +r y decay 

Experimental searches for /x + — > e + 7 have a long history of more than 50 years. These 
searches have been actively promoted by intense muon beams available at the meson facto- 
ries. Experimental efforts have been devoted to improving the detection resolutions of four 
variables, namely the positron energy (E e ), the photon energy (E y ), the timing between the 
positron and photon (Ai e7 ), and the angle between the positron and photon (A6* e7 ). Various 
kinds of apparatus have been tried in the past. In Table |IX|, several experimental results of 
90% C.L. upper limit of /i + — > e + 7 decay in the past experiments are listed along with their 
achieved detection resolutions. 

The upper limit quoted in the Particle Data Group (Particle Data Group, 1998) is 
B(fi + — > e + 7) < 4.9 x 10~ n , which was obtained by an experiment with the "Crystal Box" 
detector (Bolton, et al., 1988) at Los Alamos National Laboratory (LANL). Its apparatus 
consisted of 396 Nal(Tl) crystals, cylindrical drift chambers surrounding a muon-stopping 
target in a zero magnetic field. 

Since then, a new experimental search for fi + — > e + 7 has been carried out by the MEGA 
collaboration at LANL. A schematic view of the MEGA spectrometer is shown in Fig.|2"lj|. The 
MEGA detector consisted of a magnetic spectrometer for the positron and three concentric 
pair-spectrometers for the photon. They were placed inside a superconducting solenoid 
magnet of a 1.5 T field. The positron spectrometer comprised eight cylindrical wire chambers 
and scintillators for timing. The positron energy resolution (FWHM) was from 0.5 MeV 
(0.95%) to 0.85 MeV (1.6%) for a 52.8-MeV e + , depending on the number of helical loops 
of e + tracks. For the pair-spectrometer, each layer had lead converters, MWPCs, drift 
chambers and scintillators. The photon energy resolutions (FWHM) were 1.7 MeV (3.3%) 
and 3.0 MeV (5.7%) for the outer and inner Pb conversion layers, respectively. A surface /i + 
beam of 29.8 MeV/c was introduced along the detector axis, and was stopped in the muon- 
stopping target made of a thin tilted Mylar foil. All of the charged particles from muon 
decays are confined within the positron spectrometer. The intensity of the muon beam was 
2.5 x 10 8 /sec with a macroscopic duty factor of 6%. The total number of muons stopped 
was 1.2 x 10 14 . By using the likelihood method, a new limit of 1.2 x 10~ n with 90% C.L. 
has been reported (Brooks, et al, 1999). 
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Recently, a new experimental proposal, R-99-05, aiming at a sensitivity of 10~ 14 in 
the yU + — > e + 7 branching ratio has been approved at PSI (Barkov, et al, 1999). The 
improvement will be expected by utilizing a continuous muon beam of 100% duty factor at 
PSI. With keeping the same instantaneous beam intensity as MEGA, the total number of 
muons available can be increased by a factor of 16. Further improvement is a novel liquid 
xenon scintillation detector of the "Mini-Kamiokande" type, which is a 0.8-m 3 volume of 
liquid xenon viewed by an array of a total of 800 photomultipliers from all the sides. The 
expected resolutions (FWHM) of the photon energy and position are about 1.4% and 4 mm, 
respectively. As the e + detection, a solenoidal magnetic spectrometer with a graded magnetic 
field is adopted, in which the magnetic field is arranged so that e + from the /i + — > e + 7 decay 
follows a trajectory with a constant radius, independently of its emission angle. It allows 
easier identification of the e + in the /i + — > e + 7 decay. Physics data taking is expected to 
start in year 2003. 

A search for fi + — > e +r y r ) was also undertaken simultaneously with the /i + — > e + 7 search. 
The present 90% C.L. upper limit of -B(/i + — > e +r yy) < 7.2 x 10~ n was obtained (Bolton, et 
al, 1988). 



B. /i + — > e + e + e~ decay 

1. Phenomenology of — ► e + e + e~ decay 

The decay width of /i + — > e + e~e~ is determined from the effective Lagrangian (at the 
scale) described in Eqs. fl57j) and ( |6~T1) in Section III. The relevant interactions are 
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fi— >eee 



4Gi _ _ 

J rgi{jTR~e L ){e^e L ) + g2{jiLe R ){ele R ) 

+g^{mYe R ){el-i il e L ) + g&{jnY^L){eEl^e R ) + H.c. 



(114) 



where the Jeo and /mo photonic contributions in Eq. (p8|) are included in the four fermion 
coupling constants. 

When muons are polarized, the kinematics of the /i + — > e + e + e~ decay is determined by 
two energy variables and two angle variables of the decay positrons (Okada, et al., 1998; 
Okada, et al., 1999). The energy variables are x\ = 2E\/m il and x<i = 2E 2 /m^, where E\ 
(E 2 ) is the higher (lower) energy of the decay positrons. The allowed regions of X\ and x 2 
are 5 < x\ < 1 and 1 — x 1 < x 2 < x±, if m e is neglected compared to m M . Let us take 
the coordinate as shown in Fig.^, where the z-axis is in the direction of the decay electron 
momentum (^3), and the z — x plane is the decay plane. The positive direction of the x-axis 
is chosen to be in the hemisphere of higher-energy positron. The two angles (6, tp) determine 
the direction of the muon polarization (P M ) with respect to the decay plane. 

In this coordinate, the differential branching ratio of the fi + — ► e + e + e _ decay is given by 
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dB(ff 



e T e T e 



Ciai(x 1 ,x 2 )(l + P^cos^) + C 2 ai(xi,x 2 )(l 



P M cos ( 



dxidx 2 d(cos9)d(p 2rr 

+ C 3 {a 2 (x 1 , x 2 ) + P^/^Xi, x 2 ) cos9 + P M 7i(^i, x 2 ) sin 6 1 cos y?} 

+ C 4 {a 2 (xi,x 2 ) - P M /?i(xi,x 2 ) cos# - P A1 7i(xi,x 2 ) sin#cosy?} 

+ C 5 {a 3 (x 1 ,x 2 ) + P^/M^i,^) cos# + P M 7 2 (xi,x 2 ) sin 6 1 cosy?} 

+ C 6 {a 3 (x 1 , x 2 ) - P ll l3 2 (x l ,x 2 ) cosO - P M 7 2 (xi, x 2 ) sin # cosy?} 

+ C 7 {a 4 (xi,x 2 )(l - P fM cos9) + P M 7 3 (:ri, x 2 ) sin 6> cosy?} 

+ C s {a 4 (xi, x 2 )(l + P fl cos9) - P M 7 3 (27,x 2 ) sin 9 cosy?} 

+ C 9 {a 5 (:Ei, x 2 )(l + P M cos#) - P M 7 4 (xi, x 2 ) sin # cosy?} 

+ C 10 {a 5 (x l ,x 2 )(l - P/i cosO) + P At 7 4 (xi, x 2 ) sin 6 1 cosy?} 

+ C u P^ 3 (x l} x 2 )sm6smip - C 12 P^4(xi,x 2 )sm6smip + H.c. , (115) 

where P M is the magnitude of the polarization vector. Functions of a>i, $ and 7^ are presented 
in the Appendix C. The coefficients of C« are expressed by gi (i = 1 — 6), At, and A R . They 
are given by 
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= N 2 , c,= 


£ 6 | 2 ,C 5 = leAjjJ 2 , C 6 


= \eA L \ 2 , 
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= Re(eA R gl), 


C 8 = Re(eA L g* 3 ), C 9 


= Re(eA R g*), C w = Re(eA L g* 5 ), 
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lm{eA R gl + eA L g* 3 ), C 12 = lm(eA R g* + eA L g* 5 ). (116) 

In Eq. ( |115| ), there are four types of contributions which have different angular depen- 
dences with respect to the muon polarization. They are the contributions with an isotropic 
angular distribution, which are the even parity (P) and time-reversal (T), those proportional 
to either cos# or sin 9 cosy?, which are P-odd and T-even, and finally, those proportional to 
sin 9 sin y?, which are P-even and T-odd. 

The integrated branching ratio and the T-odd asymmetry for fi + — > e + e + e~ decay are 
given by (Okada, et al, 1999) 
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The T-odd asymmetry turns out to be proportional to Im(eA R gl + eA^gl) (= Cn) and 
Im(eA^5fg + eAx,(yf5) (= Ci 2 ). It arises from interference between the on-shell photon-penguin 
terms and the four-fermion terms. 
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If only the photon-penguin diagrams contribute to fi + — > e + e + e~ decay (namely in the 
case of C5 ^ 0, Cq ^ 0, and the others =0), a model-independent relation between the two 
branching ratios can be derived, as follows: 

B(a + ^e + e + e-) a „ ,m 2 u . 11. 

„ , , -— - ~ — In — £) = 0.006. 119 



2. Event signature and backgrounds 

The event signature of the fi + — > e + e + e~ decay is kinematically well constrained, since 
all particles in the final state are detectable. Muon decay at rest has been used in all 
past experiments. In this case, the conservation of momentum = 0) and energy 

(J2i Ei = rri/j) could be effectively used together with the timing coincidence between two 
e + s and one e~, where pi and Ei (i — 1,3) are the momentum and energy of each of the 
three e's, respectively. 

One of the physics background processes is an allowed muon decay, /i + — > e + e + e~u e 17 lx , 
which becomes a serious background when v e and have very small energies. Its branching 
ratio is B(fi + — > e + e + e~z/ e z/ M ) = (3.4 ± 0.4) x 10~ 5 . The other background is an accidental 
coincidence of an e + from normal muon decay with an uncorrelated e + e~ pair, where a e + e~ 
pair could be produced either from Bhabha scattering of e + , or from the external conversion 
of the photon in fi + — > e + h> e P fl '-f decay. Since the e + e~ pair from photon conversion has a 
small invariant mass, it could be removed by eliminating events with a small opening angle 
between e + and e~. This, however, causes a loss of the signal sensitivity, in particular for 
theoretical models in which /i + — > e + e + e~ decay occurs mostly through photonic diagrams. 

The other background, which comes especially at the trigger level, comprises fake events 
with an e + curling back to the target, which mimics an e + e~ pair. For this background, an 
e + e~ pair forms its relative angle of 180°, and can therefore be rejected. 



3. Experimental status of [i + — > e + e + e decay 



TABLE X. Historical progress and summary of searches for fi + — ► e + e + e decay. 



Place Year 90%C.L. upper limit Reference 

JINR 1976 < 1.9 x 1(T 9 Korenchenko, et al. (1976) 

LANL 1984 < 1.3 x 1(T 10 Bolton, et al. (1984) 

PSI 1984 < 1.6 x 1(T 10 Bertl, et al. (1984) 

PSI 1985 < 2.4 x 1(T 12 Bertl, et al. (1985) 

LANL 1988 < 3.5 x KT 11 Bolton, et al. (1988) 

PSI 1988 < 1.0 x KT 12 Bellgardt, et al. (1988) 

JINR 1991 < 3.6 x KT 11 Baranov, et al. (1991) 
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After the pioneering measurement in 1976 using a cylindrical spectrometer, which gave 
an upper limit of B(n + — > e + e + e~) < 1.9 x 10 -9 (Korenchenko, et ai, 1976), various ex- 
periments to search for /i + — > e + e + e~ decay have been carried out, as shown in Table [X|. 
In particular, a series of experimental measurements with the SINDRUM magnetic spec- 
trometer at SIN (Bertl, et ai, 1984; Bertl, et ai, 1985; Bellgardt, et ai, 1988) were carried 
out. A surface fi + beam with 5 x 10 6 /i + /sec was used, and the muons were stopped in a 
hollow double-cone target. The e + s and e~s were tracked by the SINDRUM spectrometer, 
which consisted of five concentric multiwire proportional chambers (MWPC) and a cylin- 
drical array of 64 plastic scintillation counters under a solenoid magnetic field of 0.33 T. 
The momentum resolution was Ap/p = (12.0 ± 0.3) % (FWHM) at p =50 MeV/c. This 
experiment gave a 90% C.L. upper limit of B(/j + — > e + e + e~) < 1.0 x 10~ 12 , assuming a 
constant matrix element for the /i + — > e + e + e~ decay (Bellgardt, et ai, 1988). They also 
observed 9070 ± 10 events of ji + — > e + e + e~u e V^ decay. A detailed analysis of the differential 
decay rate of /i + — > e + e + e~h , e V fl decay was studied, and was found to be consistent with the 
V — A interaction (Kersch, et ai, 1988). 

Another recent experiment to search for fi + — > e + e + e~ was performed at Joint Institute 
for Nuclear Research (JINR), Dubna, Russia (Baranov, et ai, 1991). A magnetic 47r spec- 
trometer with cylindrical proportional chambers was used. They obtained an upper limit of 
90% CL of -E>(/! + — > e + e + e~) < 3.6 x 10~ n , where the matrix element of /i + — > e + e + e~ was 
assumed to be constant. 



C. fi~ — e~ coherent conversion in a muonic atom 

1. Phenomenology of [i~—e~ conversion 

Another prominent process concerning lepton flavor violation is /i~ — e~ conversion in a 
muonic atom. When a negative muon is stopped in some material, it is trapped by an atom, 
and forms a muonic atom. After it cascades down in energy levels in the muonic atom, a 
muon is bound in its Is ground state. The fate of the muon is then either decay in an orbit 
{^T — > e^UfjVf,) or capture by a nucleus of mass number A and atomic number Z, namely 

fi- + (A,Z) ^^ + (A,Z-l). (120) 

However, in the context of physics beyond the Standard Model, the exotic process of neu- 
trinoless muon capture, such as 

tx- + (A,Z)^e- + (A,Z), (121) 

is also expected. This process is called e~ conversion in a muonic atom. It violates the 
conservation of the lepton flavor numbers, L e and L M , by one unit, but conserves the total 
lepton number, L. 

The branching ratio of — e~ conversion can be given by 

B(„- + (A, Z)^e- + (A, Z» S F '; f - + '*f > - e ' + *» , (122, 

T(fi- + (A, Z) — > capture) 

where F is the corresponding decay width. 
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The final state of the nucleus (A, Z) could be either the ground state or excited states. 
In general, the transition process to the ground state, which is called coherent capture, is 
dominant. The rate of the coherent capture process over non-coherent ones is enhanced 
by a factor approximately equal to the number of nucleons in the nucleus, since all of the 
nucleons participate in the process. 

The possible contributions to pr —e~ conversion in a muonic atom can be grouped into 
two parts, which are the photonic contribution and the non-photonic contribution. There- 
fore, in principle, this process is theoretically interesting, since it does occur by mechanisms 
which do not contribute to the p + — ► e + 7 process. The study of the photonic contribution 
was initiated by Weinberg and Feinberg (Weinberg and Feinberg, 1959). The non-photonic 
contribution was studied later, for instance by Marciano and Sanda (1977b). 

Let us first discuss the photonic transition. The effective Lagrangian for the photonic 
transition is written as 

£photo — ~ e JphotoA/i- (123) 

The matrix element of the /i~(p M ) — > e~ {p e )l* (?) transition, where p M , p e and q = p^ — p e 
are the muon, electron and virtual photon four-momenta respectively, is given by 



M photonic = -eA*(q) < e (pe)\Jp ho to(°)\v (P/J >> 



-eA;(q)u e (p e )[(f E0 (q 2 ) +1 J M0 (q 2 )) lu (g^ - O-, 
Hf^n{q 2 )+l,fEl{q 2 )) l -^]uM. (124) 



Based on Eq. (|124|) , the branching ratio of the coherent fi —e conversion through the 



photonic contribution is given by (Weinberg and Feinberg, 1959) 

B{prN - e-N) = (8a 5 m^ eff ZF*e) ■ (125) 

J- capt 

where T cap is the total muon capture rate. £ 2 is given by 

e = 1/BoC-mJ) + fMi{-ml)\ 2 + \fEi(-ml) + /Mo(-mJ)| 2 . (126) 

It is noted that in photonic diagrams, in contrast to /i + — > e + 7 , not only /ei and /mi, 
but also Jeo and /mo, can contribute to the fi~—e~ conversion. Z e ff is an effective atomic 
charge obtained by averaging the muon wave function over the nuclear density (Chiang, et 
ai, 1993). This is defined as 

^^1L = I d 3 x\%(x)\ 2 Pp (x) =< * u > 2 , (127) 

71 Z J 

where $^(x) is the non-relativistic muon wave function for the Is state of the muonic atom 
and Pp(x) is the proton density in the nucleus normalized as 

d 3 xp p (x) = 1. (128) 
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Fp is the nuclear matrix element squared, given by 



F p = [ d 3 xe- ipeX p p (x) = 4tt / p p (r f mmiir r 2 dr. (129) 
J J m,,r 



In Eq.( |125 ), the /x — e conversion process is roughly proportional to (Z e ff) 4 Z, whereas the 



normal muon capture, T capt , is proportional to (Z e //) 4 . The enhancement by a factor of Z 
in the e~ coherent conversion is evident. 

Let us next discuss the non-photonic contribution. The general four fermion interaction 
of the fi~—e~ conversion at the quark level is given by 



"non— photo "\/2 



q=u,d,s... L 

+ (9LP( q )eI^R + 9RP{ q )eRH-L)qibq 

+ {9LV{qjeZl ll HL + gRV^ZRH^HR)!!^ 

+ {9LA{q)eZl^PL + 9RA{q)e^YpR)qi^q 

+ o {9LT{q)eE^ v l^R + 9RT(q)eH^yL)q^^q + H.c. 



(130) 



At first, this effective Lagrangian at the quark level is converted into that at the nucleon 
level by using the nucleon form factors (Vergados, 1986; Bernabeu, et al, 1993; Faessler, et 
ai, 1999). Since the momentum transfer in the fi~—e~ conversion process is smaller than 
the size of the nucleon structure, the momentum dependence of the nucleon form factors 
can be neglected. In such a case, the matrix element of the quark currents can be replaced 
by the corresponding nucleon current by using 

<p\qT K q\p> = G^ p) pT K p, 

< n\qT K q\n > = G^'^nT K n, (131) 



with T K = (1,75,7^)7^75,^) for K = ( 5 ' P i V -> A -> T )- For the vector current, Gy' p) = 
Gy^ = 2 and Gy' p ^ = Gy' n ^ = 1. In general, the isospin invariance requires the relation 1 

G (u,p) = G (d,n)^ G (u, P ) = G (d,n) G (s,p) = G (s,n)_ ^ eflfective Lagrangian at the nuc l e oi 

level is given by 

£non-photo = —J= eZpR^{(9LS + 9lI T ^) + ($LP + 9LP T ^)l^ 

+^M(9rI + 9 { ilrs) + (9 { rp + 9^3)75}* 
+eI 7 >L^{(^ + g^rs) + (9 { la + ^3)75}* 
+eHl^R^l»{(gfy + g { ilr 3 ) + («$> + sgkhs}* 

+^a^tM R *a^{(9LT + g ( L±T 3 )}V + U^a^^a^iigW + ^r 3 )}* + He. 

where \1/ = {p,n) T and the isoscalar and isovector coupling constants g^n and g^ K (X = 
L, R, and K = S,P, V, A, T) are respectively given by 
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& = ~ E (9xK {q) G ( r ] + 9x K{q) G%>% (133) 

q=u,d,s 

9% = \ E (9xK( q )G% ;p) - gxK^G^l (134) 

For coherent p~—e~ conversion, only the scalar and vector coupling constants can be 
kept. By using a non-relativistic approximation for the muon wave function, the transition 
rate is given by 

B(jTN - e-N) = ^^(\X L (p e )\ 2 + |X^ e )| 2 )-^, (135) 

o7T 1 capt 

where 

X L (Pe) = (g^s + M + g { & + g { Ll)ZM p ( Pe ) + (g® - g<£> + - gH)NM n (p e ), (136) 

X R ( Pe ) = (<$] + g { R ] s + gtl + g { nl)ZM p (p e ) + (g§ - «jg£ + - g^)NM n (p e ) 7 (137) 

and N = A — Z is the number of neutrons in the nuclei. M p (p) andM„(p) are given by 

M p (p) = J d'xe'^p^x^^x), M n (p) = J ^xe-^pnix^^x), (138) 

with the proton and neutron densities, p p (x) and p n (x), normalized to unity. If it is assumed 
that the proton and neutron densities are equal and that the muon wave function does not 
change very much in the nucleus, by using Eq. ( |127| ) , Eq. (|135|) can be deduced into 



+| (Z + Af)( 9 ™ + <?S) + (Z - JV)( 9 « + 9 « ) I 2 } ' tt—- (139) 

J J- capt 



This equation corresponds to Eq. ( |125| ) for the photonic case, which was derived in a similar 



approximation. In Eq.( |139| ), the coherent process is enhanced by a factor of the number of 



nucleons, as can be seen in Eq. (|125|) . 

In general, both photonic and non-photonic contributions might exist. If the non- 
relativistic approximation for the muon wave function is used and the momentum transfer 
of q 2 is replaced by — m 2 , the photonic contribution can be regarded as additional terms to 
the vector- coupling constants. They are given by 

= Asg> = ^^(ho(-ml) + fMi(-ml) + /^(-mj) + M-m£)), (140) 



A^°i = Ag£l = 2 -^^(fEo(-ml) + /^(-mj) - /^(-mj) - /*i(-mj)). (141) 



These contributions are added to the corresponding vector coupling constants in Eq.flTB! 
In such the interference terms should be taken into account as well. 
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So far, the non-relativistic approximation for the muon wave function and a plane wave 
for the emitted electron have been used to derive the e~ conversion rates. Possible 
corrections for this approximation turn out to be important for heavy nuclei. A relativistic 
treatment based on the Dirac equation was considered and the corrections to the Weinberg- 
Feinberg formulas were calculated (Shanker, 1979). Recently, the photonic transitions due 
to /mi and J'ei were further examined by properly treating the electric potential in the 
muonic atom (Czarnecki, et ai, 1997). 

In the case that the photonic contributions of f'Ei{<l 2 ) and /mi(? 2 ) dominate over the 
other contributions, the rate of [i~—e~ conversion can be parameterized by (Czarnecki, et 
al, 1997) 

B(fi-N -> e~N) = 3 ■ 10 12 (|/ £1 | 2 + \f M1 \ 2 )B(A } Z), (142) 

where B(A, Z) represents the rate dependence on the mass number (A) and the atomic 
number (Z) of the target nucleus. This particular case becomes important, for instance, in 
5*0(10) SUSY GUT. The values of B(A, Z), based on different approximations, are tabulated 
in Table K|, where Bwf{A, Z) is from the Weinberg- Feinberg approximation (Weinberg and 
Feinberg, 1959), Bs(A,Z) is from Shanker (1979), and Bcmk(A, Z) is from Czarnecki, et 
al. (1997). From Eq.flTjTjp, the ratio of B(fj,+ e + -f)/B(fi-N -> e~N) is given by 

5Gu+^e+7) _ 96tt 3 « 1 428 

BlJFN -»• e~N) ~ Gfm* ' 3 • 10 12 B{A, Z) ~ B{A, Z) ' ^ ' 

By using the values in Table [XI], the ratio B(fi + — > e +, y)/B(n~N — > e^A^) for different target 
nuclei can be calculated. It varies from 389 for 27 A1, 238 for 48 Ti, and 342 for 208 Pb. This 
result indicates that the rate of — e~ conversion has a maximum around the medium nuclei 
(A ~ 60), and flattens out or slightly decreases for heavy nuclei. However, the calculations, 
which took into account the nuclear effect, show a different Z dependence (Chiang, et ai, 
1993; Kosmas and Vergados, 1996; Kosmas, et ai, 1998). 

TABLE XI. Z dependence of the photonic contribution in the fi~—e~ conversion estimated by 
various theoretical models (after Czarnecki, et al, (1997)). 



Models 


Al 


Ti 


Pb 


Reference 


B\vf(A, Z) 


1.2 


2.0 


1.6 


Weinberg and Feinberg (1959) 


B S (A, Z) 


1.3 


2.2 


2.2 


Shanker (1979) 


Bcmk(A, Z) 


1.1 


1.8 


1.25 


Czarnecki, et al. (1997) 



The /j,~—e~ conversion rates to the ground state and all excited states have been calcu- 
lated by either the shell-model closure approximation (Kosmas, et ai, 1990) or the quasi- 
particle random-phase-approximation (QRPA) (Kosmas, et ai, 1994). The fraction of the 
coherent transition to the ground state is dominant. It was calculated specifically for 48 Ti to 
be (95-99)% in the QRPA approximation, which is even larger than in the shell-model closure 
approximation. It is also experimentally advantageous, since the background from excited 
states induced by the reaction can be minimized. It was also found that among the transi- 
tions to excited states, the dipole 1~ state is large both in the photonic and non-photonic 
contributions (Kosmas, et ai, 1994). 
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2. Event signature and backgrounds 



The event signature of the coherent \x — e conversion in a muonic atom is a mono- 
energetic single electron emitted from muon capture with an energy of 



E, 



fie 



(144) 



where m M is the muon mass, and B^ and E® ec are the binding energy of the Is muonic 
atom and the nuclear-recoil energy, respectively. The nuclear-recoil energy is approximately 



El 



(m, 



- B fM ) 2 /(2M A ), where M A is the mass of the recoiling nucleus, which is small. 



Since B^ is different for various nuclei, the peak energy of the pT— conversion signal 
changes. For instance, it varies from E^ = 104.3 MeV for titanium to E^ = 94.9 MeV for 
lead. 

From an experimental point of view, e~ conversion is very attractive. First, the e~ 
energy of about 105 MeV is far above the end-point energy of the muon decay spectrum (~ 
52.8 MeV). Second, since the event signature is a mono-energetic electron, no coincidence 
measurement is required. The search for this process has a potential to improve the sensi- 
tivity by using a high muon rate without suffering from accidental background, which would 
be serious background for other processes, such as /i + — > e + 7 and fi + — > e + e + e~ decays. 

One of the major backgrounds is muon decay in orbit from a muonic atom (also called 
a bound muon decay), in which the e~ endpoint energy is the same as the energy of the 
signal. It is discussed in more detail below. The other background sources are (i) radiative 
pion capture {n~ + (A, Z) — > (A, Z — 1) + 7) or radiative muon captures (/i~ + (A, Z) — > 

+ (A, Z — 1) + 7) followed by internal and external asymmetric e + e~ conversion of the 
photon (7 — > e + e~) with e + undetected, (ii) electrons in the beam scattering off the target, 
(hi) muon decay in flight, and (iv) cosmic rays. Note that the maximum e~ energy (EJ^ ax ) 
from the background of radiative muon capture is given by 



rpmax 



"V - Bp — E rec - A z _i 



m,. 



B„ 



(145) 



where Az-i is the difference in the nuclear binding energy of the initial and final nucleus 
involved in radiative muon capture. Therefore, an appropriate target with a large A Z -i 
can be selected so as to keep a wide background-free region for the coherent signal. The 
typical values of EJ^ ax are 89.7 MeV and 91.4 MeV for 48 Ti and 46 Ti, respectively, whereas 
E^ e is 104.3 MeV. In general, to eliminate these backgrounds, the purity of the beam (with 
no contamination of pions and electrons) is crucial, together with a highly efficient veto for 
cosmic rays. 

When the muon is polarized, the angular distribution of e~ in the coherent \x~ — e~ con- 
version process is given by 



dB{n~N -> e~N) p e E e G 



d(cos 9 e ) 



16tt 



\X L {p e )\\l -P M cosfl e ) + \X R ( Pe )\\l + P^cos9 e 



capt 



(146) 
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where 6 e is the angle between the e~ direction and the muon spin direction. Since the nucleus 
does not change for the coherent process, the conversion electron carries the original muon 
spin. X L and X R are given in Eqs.( |136|) and ( |137| ), and correspond to the emission of left- 
handed electrons (eTj and right-handed electrons (e^), respectively. As in polarized /x + — > 
e + 7 decay, in principle, the angular distribution would be useful to discriminate between 
theoretical models. However, even if negative muons in the beam are 100% spin polarized, 
they are depolarized during their atomic cascades down to the Is ground state. For a nucleus 
with zero nuclear spin, the residual polarization is about 16% (Evseev, 1975). For a nucleus 
with non-zero nuclear spin, it becomes much smaller. It would make a measurement of the 
angular distribution difficult unless high statistics is accumulated. If the fi~ polarization 
is restored, however, it might provide useful information. One possible way to repolarize 
a negative muon in a muonic atom is to use a polarized nuclear target (Nagamine and 
Yamazaki, 1974; Kuno, et al, 1986). 



3. Muon decay in orbit 

Muon decay in orbit (Porter and Primakoff, 1951) is one of the important background 
sources in the search for —e~ conversion in a muonic atom, since the end point of the 
electron spectrum comes close to the signal region of fi~—e~ conversion. Only the high- 
energy end of the electron energy spectrum is of interest for [i~—e~ conversion experiments. 
At the high-energy end, the effect of the nuclear-recoil energy plays an important role (on 
its phase space). There have been several studies on its electron energy spectrum with 
nuclear- recoil energy taken into account (Hanggi, et al, 1974; Herzog, et al, 1980; Shanker, 
1982). With the approximation of a constant nuclear-recoil energy, the electron spectrum 
with an expansion in powers of the electron energy (E e ) at the end-point energy is given by 
(Shanker, 1982) 

N{E e )dE e = (^) 2 (^) 5 \D + E .(^)+F. {—)}dE e) (147) 



where 5 = E^ e —E e and Si = {m il —B il )—E rec —E e . E^ e is the e energy of the fi — e conver- 
sion signal defined in Eq. (|144|) . E rec is the nuclear-recoil energy given by E rec « E^/(2Ma). 



It should be stressed that the spectrum falls off sharply as the fifth power of Si towards 
its end point. The coefficients D, E and F as well as the end-point energy are given in a 
numerical table (Shanker, 1982). The contributions of the E and F terms to the total rate 
are about 4% and 8% respectively for Z = 29 and E e = 100 MeV. Eq.( |147| ) agrees with 



those in Hanggi, et al. (1974) and Herzog, et al. (1980). In the evaluation of the leading 
term D, important are (1) the use of a correct electron wave function incorporating the 
finite nuclear charge distribution, (2) the use of the Dirac muon wave function, and (3) the 
use of the small component of the muon relativistic wave function. In particular, the effect 
of (1) is large (Shanker and Roy, 1997). 

Experimentally, to avoid any background from muon decay in orbit, the momentum 
resolution of e~ detection must be improved. Fig.^8| shows the effective branching ratio 
of the muon decay in orbit as a function of E e for the case of a titanium target, where 
E^ e = 104.3 MeV. It was calculated using Eq.( |147| ). For a resolution better than 2%, the 



contribution from muon decay in orbit occurs at a level below 10 



-14 
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What is the asymmetric angular distribution of electrons in muon decay in orbit, if muons 
are polarized ? Numerical calculations can be made by taking into account the angular 
distribution of electrons from polarized muon decay in orbit (Watanabe, et al, 1987). It is 
given by 

N(E e ,d e )dE e (^)=N (E e )(l + a(E e )P,cosd e )(^), (148) 

where a(E e ) is the asymmetry parameter, which becomes a(E e ) = —1 at the end-point, 
giving a distinct (1 — P^cosOe) distribution; namely, electrons are likely to be emitted 
opposite to the muon polarization direction. At low energy, a(E e ) becomes positive and 
electrons are emitted along the muon polarization. The Coulomb effect is significant for 
heavy nuclei, like 208 Pb, but very small for light nuclei, like 16 0. The calculated results 
of the decay rate, emitted electron energy spectrum, and asymmetry parameters for muon 
decay in orbit are given in numerical tables for some nuclei (Watanabe, et al, 1993). 



4- Experimental status of [i — e conversion 

The SINDRUM II collaboration at PSI is carrying out experiments to search for 
e~ conversion in various nuclei. A schematic view of the SINDRUM II spectrome- 



ter is shown in Fig.gTJ. It consists of a set of concentric cylindrical drift chambers inside a 
superconducting solenoid magnet of 1.2 T. Negative muons with a momentum of about 90 
MeV/c were stopped in a target located at the center of the apparatus, after passing a CH 2 
moderator and a beam counter made of plastic scintillator. Charged particles with transverse 
momentum (with respect to the magnetic field direction) above 100 MeV/c, originating from 
the target, hit two layers of plastic scintillation arrays and then two layers of drift chambers, 
and eventually hit plexiglass Cherenkov hodoscopes placed at both ends. Charged particles 
having transverse momentum below about 100 MeV/c were contained inside, and could not 
reach the tracking region under a magnetic field of 1.2 T. A momentum resolution of about 
2.8% (FWHM) for the energy region of conversion electrons was achieved. For the back- 
ground rejection, the e~ energy (E e ), a time delay between the times of charged particle 
tracks in the spectrometer and the beam-counter signal (At), the position of the origin of 
the reconstructed trajectory (Az), the polar track angle, are used in an off-line analysis. 
Events with small At were removed so as to reject prompt backgrounds, such as electron 
scattering and radiative pion capture. 

In a 1993 run with a titanium target, a total of 3 x 10 13 stopped [i~s were accumulated 
at a rate of 1.2 x 10 7 yU^/sec from the fiEl beam line at PSI. The overall efficiency was 
about 13 %. The e~ momentum spectrum for the Ti target in the 1993 data is shown 



in Fig.^TJ, where the successive background rejections by prompt veto (i.e. At cut) and 
cosmic-ray suppression are shown. Since no events were found in the signal region, a 90% 
C.L. upper limit of 6.1 x 10~ 13 was obtained (Wintz, 1998). Also, for a lead target, it gave 
B(n~Pb -> e~Pb) < 4.6 x lO" 11 (Honecker, et al, 1996). 

A next round of the SINDRUM II experiment is under preparation at the 7rE5 beam line 
at PSI. The key element of the next stage is a pion-muon converter (PMC) to eliminate con- 
tamination of pions and electrons in the muon beam. It is needed since a veto of secondary 
pions and electrons by a beam counter will no longer be working with a high rate, such as 
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TABLE XII. History and summary of [x —e conversion in various nuclei. 



Process 






90% C.L. upper limit 


place 


year 


reference 




+ Cu- 


-> e" 


- + Cu 


< 1.6 x 1(T 8 


SREL 


1972 


Bryman, et al. (1972) 




+ 32 S- 


* 6 


,32 c 
1 l-> 


<r 7 y ID -11 


Oil > 




uajLiei locnei , et tit. i lyoi i 




+ Ti- 


> e~ 


+ Ti 


< i.6 x kt 11 


TRIUMF 


1985 


Bryman, et al. (1985) 




+ Ti- 


> e~ 


+ Ti 


< 4.6 x 1(T 12 


TRIUMF 


1988 


Ahmad, et al. (1988) 




+ Pb- 


-> e" 


- +Pb 


< 4.9 x 1(T 10 


TRIUMF 


1988 


Ahmad, et al. (1988) 




+ Ti- 


-> e~ 


+ Ti 


< 4.3 x 1(T 12 


PSI 


1993 


Dohmen, et al. (1993) 




+ Pb- 


-> e" 


' +Pb 


< 4.6 x 10" 11 


PSI 


1996 


Honecker, et al. (1996) 




+ Ti- 


-> e~ 


+ Ti 


< 6.1 x 1(T 13 


PSI 


1998 


Wintz (1998) 



10 s /x~s/sec, at the 7rE5 beam line. The PMC consists of a long-straight superconducting 
solenoid magnet with length of 8.5 m and an inner diameter of 0.4 m. It is located between 
the pion target and the SINDRUM II spectrometer, and produces the same magnetic field 
as in the SINDRUM II spectrometer, 2 T. Low-energy negative muons (called cloud muons) 
from the production target are injected into the PMC. After the 8.5-m flight length, most of 
the pions in a beam would decay out, resulting in a negligible pion contamination. Unfor- 
tunately, the original PMC magnet did not fulfill the specification at the initial installation 
stage, and thus caused a severe delay. After its new assembly, the magnet finally met the 
goal. With a lower magnetic field of the PMC magnet, data with a gold target were taken in 
1998. A new run for a Ti target is expected to start in spring, 1999, with aimed sensitivity 
of B{iTTi -> e'Ti) < a few x 1(T 14 . 

A new experiment, E940, at Brookhaven National Laboratory (BNL) AGS, the MECO 
(Muon Electron Conversion) experiment, was prepared (Bachman, et al., 1997). MECO 
aims to search for fi~ + Al — >• e~ + Al at a sensitivity below 1CU 16 . It will use a new high- 
intensity pulsed muon beam, which could yield about 10 11 /i~s/sec stopped in a target. A 
schematic layout of the MECO detector is shown in Fig|?T|. The MECO apparatus consists 
of a superconducting (SC) solenoid magnet to capture pions from the production target 
(production solenoid), a curved transport SC solenoid magnet system (transport solenoid), 
and a SC solenoid spectrometer, which observes only the 105-MeV signal electrons (detector 
solenoid). Based on the solenoid capture scheme originally proposed by MELC (Dzhilkibaev 
and Lobashev, 1989), it has an axially graded magnetic field (from 3.5 T to 2.0 T) to effi- 
ciently capture pions from a tungsten target located on the axis of the solenoid magnet. The 
curved transport solenoid will capture muons from pion decays, and select the momentum 
and sign of charged particles by using collimators at three positions. Layers of thin aluminum 
targets where /i~s are stopped are placed in the detector solenoid with an axially graded 
magnetic field. The conversion electron of 105 MeV is momentum analyzed with a resolution 
of 300 keV (RMS) and an acceptance of 25% in a straw tracking chamber. A pulsed proton 
beam of about 1 MHz repetition with a pulse length of 30 nsec can be extracted at the 
AGS. A high extinction between the beam pulses (the ratio of a number of protons between 
pulses to that in the beam pulse) of 10~ 9 is needed to eliminate severe beam background at 
a high rate. They expect to observe 6 signal events for B(fi~Al — > e~Al) w 10 -16 during a 
one-year run, with an expected background of 0.4 events. 
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D. y~ — e + conversion in a muonic atom 

1. Phenomenology of [i~ — e + conversion 

The other neutrinoless muon-capture process is a charge-changing reaction, such as 

y7 + (A, Z)^e + + (A, Z - 2)*, (149) 

which violates the conservation of the total lepton number as well as the lepton flavor 
numbers, L e and L^. This process is closely related to neutrinoless double (3— decay (/3/3 „), 
since both processes require a mechanism involving two nucleons. The final state of the 
nucleus (A, Z — 2)* could be either the ground state (gs) or excited states (ex). Since the 
final nucleus is not the same as the initial nucleus, no coherent enhancement, even for the 
transition to the ground state, is expected. The branching ratio of the y~—e + conversion is 
defined by 

w - + <A 7\ + + <A7 nW _ r(^ + (A,Z)^e + + (A,Z-2y) 

B(y +(AZ)^e + (A,Z-2)) = r(^~ + (A, Z) -> capture) ' (150) 

Various theoretical models predict the rates accessible experimentally. One is the mini- 
mum super symmetric models (MSSM) with R-parity breaking, which allows the predicted 
branching ratio of the y~— e + conversion on the level of 10 -12 , since the A and A' parameters 
involved are not constrained (Babu and Mohapatra, 1995). The left-right symmetric models 
with a low-mass Wr also predict the y~—e + conversion branching ratio of 10 -14 , estimated 
by the same authors. 



2. Event signature and backgrounds 
The energy of the positron from the y~— e + conversion is given by 

Efj, e + — m p. — Bp — E rec — A Z -2 

^m^-B^- A z _ 2 , (151) 

where Az-2 is the difference in the nuclear binding energy between the (A, Z) and (A, Z — 2) 
nuclei, with the excitation energy in the final nucleus taken into account. In y~—e + con- 
version, because of the absence of coherent enhancement, the final nucleus could be either 
in the ground state or in excited states. Usually, it is assumed that a large fraction of the 
final nucleus could be in the giant dipole resonance state, which has a mean energy of 20 
MeV and a width of 20 MeV. Therefore, the e + from y~—e + conversion will have a broad 
momentum distribution corresponding to the width of giant dipole resonance excitation. 

The major background is radiative muon capture (RMC) or radiative pion capture 
(RPC), followed by asymmetric e + e~ conversion of the photon. For some nuclei, the end- 
point of the RMC background in Eq.( |145| ) can be selected to be much lower than the signal. 
For instance, for a titanium target, the maximum endpoint of RMC (of about 90 MeV) is 
about 10 MeV lower than the signal energy of about E^ e + w 100 MeV. The background 
from RPC must be reduced from the rejection of pions in the beam. 
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3. Experimental status of \i —e + conversion 



TABLE XIII. Historical progress and summary of e + conversion in various nuclei, gs and 
ex denote the transitions to the ground state and excited states (mostly giant dipole-resonance 
states) respectively. 



Process 






90% C.L. upper limit 


place 


year 


reference 




+ Cu 


— > e 


+ + Co 


2.6 x 10" 8 


SREL 


1972 


Bryman, et al. (1972) 




+ s- 


e+ 


+ Si 


9 x 1(T 10 


SIN 


1982 


Badertscher, et al. (1982) 




+ Ti- 


-e* 


' + Ca(gs) 


9 x KT 12 


TRIUMF 


1988 


Ahmad, et al. (1988) 




+ Ti- 


-e^ 


' + Ca(ex) 


1.7 x KT 10 


TRIUMF 


1988 


Ahmad, et al. (1988) 




+ Ti- 




' + Ca(gs) 


4.3 x 1(T 12 


PSI 


1993 


Dohmen, et al. (1993) 




+ Ti- 




' + Ca(ex) 


8.9 x KT 11 


PSI 


1993 


Dohmen, et al. (1993) 




+ Ti- 




- + Ca(gs) 


1.7 x 10~ 12 


PSI 


1998 


Kaulard, et al. (1998) 




+ Ti- 




' + Ca(ex) 


3.6 x KT 11 


PSI 


1998 


Kaulard, et al. (1998) 



The SINDRUM II Collaboration at PSI has reported on a search for the charge-changing 
yT-\-Ti — > e +J rCa in muonic atoms (Kaulard, et al., 1998). It was carried out simultaneously 
with a measurement of fi~ + Ti — > e~ +Ti. The e + momentum spectrum is shown in Fig.|32|. 
The results are given separately for the transition to the ground state and that to the giant 
dipole resonance. They are summarized in Table |X111|, together with the previous results. 



E. Muonium to anti-muonium conversion 



A muonium atom is a hydrogen-like bound state of fi + and e~ . The spontaneous conver- 
sion (or oscillation) of a muonium atom (/i + e _ or Mu) to its anti-atom, anti-muonium atom 
(/i~e + or Mu) is another interesting class of muon LFV process. In this Mu — Mu conver- 
sion, the lepton flavors change by two units (AL e /^ = ±2) in the ordinary law of separate 
additive muon and electron numbers, whereas it would be consistent with multiplicative 
muon or electron number conservation (Feinberg and Weinberg, 1961). The possibility was 
suggested by Pontecorvo in 1957 (Pontecorvo, 1957), even before the muonium atom was 
observed for the first time at the Nevis cyclotron of Columbia University (Hughes, et al., 
1960). 



1. Phenomenology o/Mu-Mu conversion 



Various interactions could induce the |AL$| = 2 processes, such as Mu — Mu conversion, 



as discussed in Section |III E[ . To discuss the phenomenology of the Mu — Mu conversion, as 
an example, the effective four fermion interaction of the (V — A)(V — A) type (Feinberg and 
Weinberg, 1961) is taken. It is given by 

Huum = (^p)/!7A(l - 75)e/l7 A (l - 7s)e + H.c. (152) 



61 



in which G Mu y^ is a coupling constant characterizing the strength of the interaction. 

In the absence of an external magnetic field, the muonium and the anti-muonium have 
the same ground-state energy levels. The possible new interaction in Eq. (|152 ) would cause 
a splitting of their energy levels of 

^ 2 <M|^|M>=^(%S), (153) 

where n is the principal quantum number of the muonium atom, and do is the Bohr radius 
of the muonium atom. For the ground state of the muonium atom (n = 1), 

5 = 1.5 x 1(T 12 • (eV). (154) 

^ Crp 

The Mu — Mu conversion is analogous to K° — K° mixing. If a muonium atom is formed 
at t = in a vacuum under no external electromagnetic field, it could oscillate into an anti- 
muonium atom with time. For a small 5 value, the probability (p Mu mh) is approximately 
given by (Willmann and Jungmann, 1998) 

PuuMlit) = sm 2 (|) ■ A M e- A ^ « (| )* • \,e~^\ (155) 

where A M = l/r fl (= 2.996 x 10~ 10 eV) is the muon decay width. The maximum probability of 
anti-muonium decay is t max = 2r M . Fig.|33] shows the oscillation pattern as a function of time. 

3>0 

MuMu 



The total conversion probability after integration over time (P° ^r-j-) in a zero magnetic field 



is 

P^m = jf Pu«m(t)dt = 2(| ^ |2 'f | ^ |2) = 2.56 x 1CT 5 ■ (^) 2 - (156) 

The experimental limit constrains the upper limit of magnitude of G Mu j^. The limit of 
G^MuMu is improved by the square root of the conversion probability. 

The presence of an external electromagnetic field would remove the degeneracy between 
the muonium and the anti-muonium atoms. It would reduce the probability of the muonium 
to anti-muonium conversion. The splitting of different muonium energy levels in the presence 
of a magnetic field is calculated by using the Breit-Rabi formula for the states of their total 
spin, F, and its component, mp. In a magnetic field, the (F,mp) = (1,±1) — > (1,±1) 
transitions become rapidly suppressed, even at a weak field, because of the Zeeman splitting 
of energy levels. The transitions between different F states are also highly suppressed, even 
in a zero magnetic field, owing to the muonium Is hyperfine splitting (of 1.846 x 10~ 5 eV). By 
taking into account the magnetic-field dependences of different energy levels of muonium and 
anti-muonium and their transition rates, Eq.( |156| ) can be modified for unpolarized muons 
by 

P Mum( B ) = 7 I! 9 /|A|2 Z I A 12 ZT\ \2\ = P MuVhl ' S b(B), (157) 



4^ F 2(|^ + |AP+|A MI , 

where A = E Mu (F, mp) — Ej^(F, mp), and 5 and A are functions of the magnitude of the 
magnetic field (B). The reduction factor, Sb(B), has been calculated for possible interac- 
tions of different types (Hou and Wong, 1995; Horikawa and Sasaki, 1996). Fig.|34| shows 
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the dependence of the Mu — Mu conversion probability on the external magnetic field and 
different coupling types. For example, for the traditional (V — A)(V — A) interaction, the 
conversion rate becomes one half at a magnetic field of about 10 mG and is further strongly 
suppressed for a magnetic field greater than 10 3 G. 



2. Event signature and backgrounds 

In experiments of the Mu — Mu conversion, an anti-muonium converted from a muonium 
initially produced is searched. The experimental signature of an anti-muonium decay is the 
emission of an energetic e~ from \x~ — > e~V [1 y e decay with a dissociated e + left behind with 
an average kinetic energy of 13.5 eV. This corresponds to the binding energy of the Is state 
of a muonium atom. 

The sensitivity to Mu — Mu conversion is known to be suppressed when the muonium 
atom is in matter. This occurs since a negative muon in anti-muonium is easily captured by 
surrounding atoms. Therefore, recent experiments have been performed by using muonium 
atoms in a vacuum. 

There are two major sources of potential backgrounds. One is accidental coincidences 
of energetic e~ produced by Bhabha scattering of e + from /i + decay in a muonium and the 
scattered e + . The second is the physics (prompt) background from the fi + — > e + b , e T> fJj e + e~ 
decay (whose branching ratio is 3.4 x 10~ 5 ), when the e~ becomes energetic and only one of 
the two e + s is detected. 



3. Experimental status of Mu — Mu conversion 



TABLE XIV. Historical progress and summary of Mu — Mu conversion. 



Place 


Year 




Reference 


TRIUMF 


1982 


< 42 


Marshall, et al. (1982) 


TRIUMF 


1986 


< 20 


Beer, et al. (1986) 


TRIUMF 


1990 


< 0.29 


Huber, et al. (1990) 


LANL 


1991 


< 0.16 


Matthias, et al. (1991) 


LANL 


1993 


< 6.9 


Ni, et al. (1993) 


PSI 


1996 


< 0.018 


Abela, et al. (1996) 


JINR 


1997 


< 0.14 


Gordeev, et al. (1997) 


PSI 


1999 


< 0.003 


Willmann, et al. (1999) 



The historical progress in the searches for Mu — Mu conversion is listed in Table |XIV| . A 
recent experiment was carried out at PSI (Willmann, et al., 1999). The experiment fully uti- 
lized the powerful techniques developed at the previous experiment at LANL (Matthias, et 
al, 1991), which requires the coincidence identification of both particles in the anti-muonium 
decay. Its experimental setup is shown in Fig.|3"5j. Muonium atoms were produced by stop- 
ping surface muons in a SiO"2 powder target, where some fraction diffused out through the 
target surface with thermal energy in a vacuum. To detect e~ from fi~ decay, a magnetic 
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spectrometer was used. It consisted of five concentric multiwire proportional chambers with 
64 segmented hodoscopes at a 0.1 T magnetic field. The e + with an average kinetic energy 
of 13.5 eV was detected by micro-channel plate detectors after electrostatic acceleration to 
8 keV. With the production of about 5.7 x 10 10 muonium atoms, their analysis yielded one 
event satisfying all of the required criteria with the expected background events of 1.7 ±0.2 



due to accidental coincidence. The Monte-Carlo data and real data are given in Fig.p6|. The 
90% CL. upper limit on the conversion probability at zero magnetic filed is 

P uum < 8 - 3 >< lO- n /S B (B), (158) 

where the factor Sb(B) describes the suppression of the Mu — Mu conversion in an external 
magnetic field, B. It could be translated into the upper limit on the effective coupling 
constant, G Mu ^, which is given by 

GunWl < 3-0 x 10~ 3 G F (159) 
at 90% CL. upper limit under a 0.1 T magnetic field. 
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VI. FUTURE PROSPECTS 



The field of muon decay physics is presently very productive, even after its long history of 
over 60 years. Currently, there are several new experiments which are being either prepared 
or planned. Some of them, which were mentioned in this article, are R77 at RIKEN-RAL and 
R-99-07 at PSI for the muon lifetime measurement fin Section IIV A 21), E614 at TRIUMF to 



measure the Michel spectrum and its asymmetry (in Section IVB 2 ), R-94-10 and R-97-07 



at PSI to measure the e + polarization in polarized /z + — > e + z/ e z/ M decay (in Section |1V C 2j ), 
R-99-05 at PSI for //+ -> e+7 decay (in Section |Vjg ), the new phase of SINDRUM-II at 
PSI, and E940 (MECO) at BNL for n~-e~ conversion (in Section [VTT^ ). Each of them 
is aiming at an improvement of about an order of magnitude or more over the previous 
experiments. The potential progress expected by each of such experiments is based not only 
on innovative ideas on detection methods, but also on muon beams of high intensity and 
good quality. In particular, the planned searches for muon LFV processes strongly rely on 
the beam, such as the PMC in SINDRUM-II and the superconducting solenoid capture and 
transport systems for the MECO experiment. More muon fluxes with less contamination 
are critical for further improvements. 

Currently, two out of the three meson factories are operational. One of the two opera- 
tional machines, the PSI cyclotron, has increased its proton current, achieving 1.5 mA, the 
highest proton current in the world. The muon beam intensities for various existing labora- 
tories are listed in Table fX"Vl In addition, the use of higher energy proton machines, such as 



the BNL AGS for negative muons, is being considered for the MECO experiment, where a 
pulsed-beam capability at the AGS and a larger cross section of negative pions at a few 10 
GeV proton energy are to be utilized. In the long-term future, there are several new projects 
to construct high-intensity proton accelerators: the JAERI/KEK Joint Project (previously 
JHF) (JAERI/KEK Joint Project, 1999), which consists of a 50-GeV proton synchrotron 
(50-GeV PS) with a 15 /iA beam intensity, and a 3-GeV proton synchrotron with a 330 //A 
beam intensity; the SNS (Spallation Neutron Source) at Oak Ridge; a possible European 
Spallation Neutron Source (ESS). Probably, a proton driver for a collider (Muon Col- 

lider Collaboration, 1996; Ankenbrandt et al, 1999) can be included in the long-term future. 
Note that among the above, only the 50-GeV PS is planned to have a continuous proton 
beam by slow beam extraction, whereas the others may have only fast beam extraction of a 
low repetition rate. 

Regarding LFV, besides the study of muon decays, a unique possibility of lepton flavor 
changing Rutherford scattering, where the conversion of incident electrons into muons of the 
same energy by scattering in the external electric field of a massive nucleus, is also discussed 
(Abraham and Lampe, 1996). However, the expected cross section is very small to compete 
with rare muon decay processes, and technical details have not yet been discussed. 



A. Towards new high-intensity muon sources 

Significant improvements in low-energy muon physics could be expected if a high- 
intensity muon source, having a beam intensity of 10 12 — 10 13 /i ± s/sec with a narrow energy- 
spread and less contamination, can be realized. The muon beam intensity envisaged would 
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TABLE XV. Intensities of existing muon beams available. 



facility 


Protons 


Time structure 








Muon Intensity 






PSI 


1.5 mA 


continuous 


3 x 10 8 




s/sec 


at 28 MeV/c 


(surface 


muons) 




590 MeV 


(50 MHz) 


1 X 10° 




s/sec 


at 100 MeV/c 




TRIUMF 


150 fiA 


continuous 


2 X 10 7 


/' H 


s/sec 


at 28 MeV/c 


(surface 


muons) 




500 MeV 


(50 MHz) 


3 x 10 6 


A 


"s/sec 


at 100 MeV/c 




RAL 


200 fiA 


pulsed 


1 x 10 6 




"s/sec 


^ at 28 MeV/c 


(surface 


muons) 




800 MeV 


(50 Hz) 


1 x 10 5 


fJ - 


"s/sec 


at 50 MeV/c 






MSL 2 ) 


6 fik 


pulsed 


1 x 10 5 


M H 


"s/sec 


at 28 MeV/c 


(surface 


muons) 




500 MeV 


(20 Hz) 


1 x 10 4 


/i" 


"s/sec 


at 55 MeV/c 






JINR 3 ) 


4 fiA 


continuous 


3 x 10 4 


M H 


"s/sec 


at 28 MeV/c 


(surface 


muons) 




660 MeV 




1 x 10 3 




"s/sec 


at 100 MeV/c 





^ The highest instantaneous intensity of 2 x 10 4 /i + s/200 nsec. 

2 ) Meson Science Laboratory at KEK, using the existing 500 MeV Booster ring. 

3 ) Phasotron, Joint Institute for Nuclear Research, Dubna, Russia. 



be four or five orders of magnitude higher than that available today. Ideas of such a high- 
intensity muon source are based on (i) solenoid pion capture, (ii) phase rotation, and (iii) 
muon cooling. A schematic view of the basic concept is shown in Fig.[37[ 

In solenoid pion capture, low-energy pions and muons are trapped in a high solenoidal 
magnetic field (such as 10 T or more). From Monte Carlo simulations with appropriate pion 
production cross sections, about 0.3 to 0.1 captured pions (of less than 0.5 GeV/c) per proton 
are estimated for proton beam energies of 50 GeV to 10 GeV, respectively. For a proton 
intensity of the existing and planned proton machines of about 10 13 — 10 14 protons/sec, a 
large number of captured pions sufficient for the aimed intensity are expected. 

The phase rotation is to accelerate slow muons and to decelerate fast muons by a strong 
radio-frequency (RF) electric field, yielding a narrow longitudinal momentum spread. To 
identify fast and slow muons by their time of flight from the production time, a very narrow 
pulsed proton beam must be used. An intensity enhancement of a factor of about ten in the 
longitudinal energy distribution is estimated from Monte Carlo simulations (Kuno, 1997b). 

The muon cooling, which is based on ionization cooling (Skrinskii and Parkhomuchuk, 
1981), is to reduce muon beam emittance. The longitudinal cooling is based on the repe- 
tition of energy loss by ionization and subsequent acceleration to restore the longitudinal 
momentum. The ionization cooling works only for muons. 

These ideas have emerged in studies of a yU + yU~ collider at the high-energy frontier (Muon 
Collider Collaboration, 1996; Ankenbrandt et ai, 1999). The physics potential with low 
energy muons available from the front-end of the /x 4 ' [i~ collider complex has been discussed. 
Although there are many common R&D items between a low-energy muon source and a 
collider, there are discussions on whether the front-end muon collider (FMC) could 
be directly used in experiments with muons. The FMC will run with a pulsed beam of slow 
repetition (at typically 15 Hz). However, most experiments with muons require a beam with 
a high duty factor, or a nearly continuous beam, because of the reduction of the instantaneous 



66 



rate (Molzon, 1997). The precise requirement on the beam time structure depends on the 
type of experiments. For instance, searches for fi + — > e + 7 and fi + — > e + e + e~ must use a 
continuous beam to reduce the instantaneous rate, whereas searches for e~ (or fi~—e + ) 
conversion, Mu — Mu conversion and a measurement of the muon lifetime need a pulsed 
beam with a pulse separation of an order of the muon lifetime (~ /isec). Thus, independent 
R&D items, in particular concerning phase rotation and the muon cooling system, exist in 
a low-energy muon source. These technical issues must be pursued separately. 

There are several dedicated R&D programs on low-energy muon sources with high in- 
tensity. One of those is the PRISM project at KEK in Japan (Kuno, 1998). The PRISM 
project, which is an acronym of Phase- Rot at ion Intense Secondary Meson beam, would 
combine solenoid capture, phase rotation and possibly modest muon cooling to produce 
a cooled muon beam. The requirements on muon cooling as a secondary beam is not as 
strict as that in the collider. Its R&D program starts from a relatively low repetition 

rate (~ kHz), and aims at a higher repetition in the future. The others are the MUONS 
project at TRIUMF (Blackmore, et al, 1997), and the Super-Super Muon Channel project 
at the RIKEN-RAL muon facility (Ishida and Nagamine, 1998). For the last project, a 
new scheme of the production of cooled /i + s by laser ionization of thermal muonium is also 
proposed (Nagamine, 1996). 

Once a highly intense muon source with a narrow energy spread and less contamination 
is available, physics programs with stopped muons, in particular searches for rare muon LFV 
processes, would be significantly improved. First of all, the potential sensitivity achievable in 
searches for rare processes is ultimately limited by the number of muons available. Therefore, 
a high-intensity beam is essential. Small beam contamination is necessary to further reduce 
any background associated with it. A narrow energy spread of the beam will allow a thin 
muon stopping target to improve the detection resolution. For instance, if about 10 19 — 10 20 
muons/year are available, a new experiment on pT — e~ conversion with a sensitivity of 1CT 18 
could be possible (Blackmore, et al., 1997). 

The high-intensity muon sources could be used not only for experiments with low-energy 
muons, but also for experiments with energetic muons if the muons thus produced are in- 
jected into additional accelerators for further acceleration. Potential programs might include 
the measurements of the muon anomalous magnetic moment and the muon electric dipole 
moment, and also a muon accumulator ring for neutrino sources (Geer, 1998; CERN Report, 
1999). In addition to particle-physics programs, a broad research field from materials science 
to biology would benefit from new highly intense muon sources, which would definitely open 
up a new era of muon science. 
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VII. CONCLUSIONS 



We have described the current theoretical and experimental status in the field of muon 
decay to search for physics beyond the Standard Model. Among many interesting topics 
of physics related to muons, we have discussed the precise measurements of normal muon 
decay, and the searches for muon LFV processes. In particular, we have highly stressed the 
importance of muon LFV processes. 

The physics motivation for LFV is extremely strong. LFV has recently attracted much 
attention from theorists and experimentalists, more than ever. This has happened because 
SUSY models predict large branching ratios for LFV processes. Their predictions are just 
as large as one or two orders of magnitude lower than the present experimental limits. They 
could be accessible and tested by future experiments. There are many scenarios of SUSY 
models which predict sizeable LFV effects. They are such as SUSY GUT, SUSY with right- 
handed (heavy) Majorana neutrinos, SUSY with i?-parity violation and others. They would 
provide an opportunity to give a hint on physics at very high energy scales, like either the 
GUT scale or the mass scale of a heavy right-handed Majorana neutrino for the see-saw 
mechanism. Of course, there are many other classes of theoretical models which predict a 
large LFV effect besides SUSY. Therefore, LFV searches have robust potential to uncover 
new physics beyond the SM. 

We have presented the phenomenology of muon LFV processes of |ALj| = 1, such as the 
/i + — > e + 7 and fi + — > e + e + e~ decays, the fi~—e~ conversion in a muonic atom, and those of 
|ALj| = 2 such as the muonium to anti-muonium conversion. We discussed all of the above 
processes of |ALj| = 1 within the same framework of the effective Lagrangian to illustrate 
how various contributions (such as photonic and non-photonic) can be disentangled with 
measurements of the three muon LFV processes. Thereby, searches for these three processes 
are equally important. If the muon is polarized, additional information could be obtained 
by measuring the angular distributions in the fi + — > e + 7 decay and fi~—e~ conversion, 
and T-odd and P-odd correlations in the fi + — > e + e + e~ decay. Furthermore, for the fi + — > 
e + 7 decay, the use of polarized muons would be useful to eliminate any background processes 
in the search. Experimentally, positive muons in a surface muon beam are known to be 100% 
polarized, and thereby the use of polarized muons will be feasible in the future. Then, we 
have briefly mentioned the most recent experimental results, such as MEGA for fi + — > 
e + 7 decay, SINDRUM-II for fi~—e~ conversion, the recent search for Mu — Mu conversion 
at PSI, and others, together with future experimental prospects. In addition, precision 
measurements of the normal muon decay, ji + — > e + z/ e 77 M , have attracted much interest. In 
the near future, new measurements of the muon lifetime, of the Michel spectrum and its 
asymmetry, and of the e + polarization in polarized fi + — > e + u e U lx decay will be carried out 
while aiming at an order of magnitude improvement. 

A new intense muon source with 10 12 — 10 13 /i ± s/sec would be strongly required to make 
substantial improvements in low-energy muon physics. The aimed intensity is four or five 
orders of magnitude higher than that available now. The ideas for such a muon source are 
based on (i) solenoid pion capture, (ii) phase rotation, and (iii) muon cooling. These ideas 
came from a fi + fi~ collider. However, the beam time structure must be of high duty factor 
for low-energy muon physics, leading to new technical challenges which do not exist for the 
collider R&D studies. To overcome these issues, several R&D programs dedicated to 
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low-energy muons are now being undertaken at KEK, RIKEN and TRIUMF. With increased 
muon fluxes, the searches for rare muon LFV processes, as well as precision measurements 
of muon decay, are expected to be significantly advanced. 

Muon physics becomes important with strong physics motivations. There are extraordi- 
nary opportunities which will allow us to explore discovery potentials of physics beyond the 
SM, with low-energy muons. 
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APPENDIX A: RADIATIVE MUON DECAY 



The differential branching ratio of the radiative muon decay, fi + — ► e + v e V^ , is given in 
Eq. (j52|) . The functions appearing in Eq . (|52|) , F(x, y, d), G(x, y, d) and H (x, y, d) in the SM, 
are given as follows: 

p = p(o) +r _pW + r 2 J F( 2 ) ; 
G = G® +rG« + r 2 G^, 

H = H^ + rH^+r 2 H^, (Al) 

where r = [vn e jm n ) 2 . m e and are the masses of an electron and a muon, respectively. 
Here, x and y are the normalized electron and photon energies, x = 2E e /m fJL and y = 
2E 1 /m ll ; d is given by d = 1 — ftp e -p-y- p e and p 1 are unit momentum vectors of the electron 
and the photon respectively, ft is defined as ft = |p^|/i£ e . 



p{o) 



pW 

F m 

G (o) 

G (l) 

G& 
H (o) 



x, y, d) = -{y 2 (3 - 2y) + 6xy(l - y) + 2x 2 (3 - 4y) - 4x 3 } 

+8{-.xy(3 - y - y 2 ) - x 2 (3 - y - Ay 2 ) + 2x 3 (l + 2y)} 

+2d{x 2 y(6 -by- 2y 2 ) - 2x 3 y(4 + 3y)} + 2d 2 x 3 y 2 (2 + y) (A2) 

x, y, d) = §{- y(3 ~ %) " (3 - 4y) + 2a;} + ^{y(6 - 5y) - 2x(4 + y) + 6x 2 } 

+8{x(4 - 3y + y 2 ) - 3x 2 (l + y)} + 6dx 2 y(2 + y) (A3) 

32 (4-3J/) o X 48y 

x,y,d) = — { 3 } + -T A4 

d z x d 
g 

x , y , d) = -{xy(l - 2y) + 2x 2 (l - 3y) - 4x 3 } + 4{-x 2 (2 - 3y - 4y 2 ) + 2x 3 (2 + 3y)} 

(jj 

-4dx 3 y(2 + y) (A5) 

32 8 
x, y, d) = — (-1 + 2y + 2x) + -(-xy + 6x 2 ) - 12x 2 (2 + y) (A6) 

96 

x,y,d) = -— (A7) 
g 

x, y, d) = -{y 2 (l - 2y) + xy(l - 4y) - 2x 2 y} + 4{2xy 2 (l + y) - x 2 y(l - 4y) + 2x 3 y} 
+2d{x 2 y 2 (l - 2y) - 4x 3 y 2 } + 2d 2 x 3 y 3 (A8) 

x, y, d) = + 2y} + ^{y(2 - 5y) - xy} + 4xy(2y - 3x) + 6rfx 2 y 2 (A9) 

d z x d 

96y 48y /A . 



APPENDIX B: MSSM LAGRANGIAN 

The Lagrangian for the minimal supersymmetric standard model (MSSM) is described. 
In SUSY theories, elementary fields are introduced as a pair of bosonic and fermionic fields. 
Such a pair is called a supermultiplet. There are two types of supermultiplets, a gauge 
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multiplet and a chiral multiplet. A gauge multiplet consists of a gauge field (Afy and its 
superpartner, a gauge fermion (or gaugino) field (A a ), which is a Majorana fermion field in 
the adjoint representation of the gauge group. In the MSSM, we have to introduce gaugino 
fields for SU(3)c, SU{2) L and U(l)y gauge groups. A chiral multiplet is a set of a complex 
scalar field (0) and a left-handed Weyl fermion field Its complex conjugate is called 

an anti-chiral multiplet which consists of 0* and i^r. In the MSSM, these chiral multiplets 
correspond to matter fields, namely quark, lepton and Higgs fields and its superpartners. 
The fields in the same chiral multiplet have the same quantum numbers for the gauge 
groups. Chiral multiplets necessary for the MSSM are listed with their quantum numbers 
in Table [XV J| . The right-handed squarks and sleptons are defined as follows: (u^ L )* = um, 
(d c iL )* = dm, and (e£ L )* = em- We sometimes use a notation and for superpartners of 
and ipL and <3> for a supermultiplet (0, -0l)- 



TABLE XVI. Chiral multiplets in MSSM. 





Qi(qiL,qiL) 






Li(j iL ,liL) E^UAl) 






SU(3) C 


3 


3 


3 


1 1 


1 


1 


SU(2) L 


2 


1 


1 


2 1 


2 


2 


U(1) Y 


l 

6 


2 
3 


i 

3 


-1 1 


i 

2 


i 

2 



In phenomenological applications of SUSY model, the SUSY Lagrangian consists of two 
parts, namely the SUSY invariant Lagrangian and the soft SUSY breaking terms. 

£ = C-SUSY inv + C-SUSY breaking (Bl) 

The SUSY invariant part of the MSSM Lagrangian is given as follows. Csusy inv can be 
decomposed into two parts, 

f-'SUSY inv C-gauge ~t~ £ super •'potential-, (B2) 

C-gauge depends on the gauge coupling constants, and is given by 

Cgauge= £ F<$ + £ (^7 ' ^iL + I^W) 

gauge multiplet chiral multiplet 

£ gaugino— matter C*J) term, (B3) 

where for each gauge group the C ga u g ino-matter and C D term terms are given by 

£gaugino-matter = ~£ ^Wi^^^iL + H.C., (B4) 



ID term = "EyE^W' ( B5 ) 



In addition to the normal gauge interactions defined in the covariant derivatives, these two 
types of interactions specified by the gauge coupling constants in Eqs. (p4|) and (|B5|) are 
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necessary to keep the SUSY invariance of the Lagrangian. C superpo tentiai is determined from 
the superpotential W(4>i), which is a function of scalar fields of the chiral multiplets, 



W((j)). 2 ld 2 W((p)- 



The superpotential therefore generates a set of bosonic and fermionic interactions. From 
the gauge invariance, the superpotential {W(<f>i)) for the MSSM is given by 

Wmssm = {y e ) l3 U x EtL, + {y d ) ij H l DiQ j + {y^H^fQj - »H X H 2 , (B7) 

where the contraction of the SU (2) indices is made by using the anti-symmetric tensor, 
e a p. Also, the i?-parity conservation is required (Section [111 C 3| ). By substituting these 
functions in Eq. (|B6|) , this superpotential would induce the ordinary Yukawa couplings and 



the higgsino mass terms as well as the other interactions which are necessary to ensure the 
SUSY invariance. 

The soft SUSY breaking mass terms are defined as terms which do not introduce 
quadratic divergence, and essentially serve as mass terms for superpartners. A general 
form of SUSY-breaking terms in the MSSM is given by 

£-SUSY breaking = ~{ m l)ij^*Ri^Rj ~ ( m 'l)ijlj J jLj ~ ( m d) ij^Ri^Rj 

-( m l)iAu Rj - (rrfyijciliqu - m 2 H ^HlH x - m 2 H2 H*H 2 
-[moiA^ijH^ej + mQ{A d ) ij H l DlQ j 
+m (A u ) ij H 2 U*Q j - iiBE x E % 



+-M X B R B L + -M 2 W R W L + -M 3 G R G L + h.c]. (B8) 

These terms are quadratic terms for scalar quarks, leptons and Higgs fields, scalar tri- 
linear terms (A terms) and gaugino Majorana mass terms. These terms are supposed to be 
generated from spontaneous symmetry breaking of SUSY, presumably at some high energy 
scale in some sector outside the MSSM dynamics, such as the supergravity or the dynamical 
SUSY breaking sector. For more details, see for example Nilles (1984) or Haber and Kane 
(1985). 



APPENDIX C: DIFFERENTIAL BRANCHING RATIO OF THE 

(i+ -> e+e+e" DECAY 



The differential branching ratio of the /i + — > e + e + e~ decay is given in Eq.( |115j ). The 
kinematical functions appearing in Eq.( |115[ ), ai(x x ,x 2 ), /?j(xi,x 2 ) and Ji(x x ,x 2 ), are given 
as follows (xi = 2Ei/m a (i = 1,2) are the normalized e + energies): 



a x (x x ,x 2 ) = 8(2 — x x — x 2 ){x x + x 2 - 1), 

a 2 (x x , x 2 ) = 2{x x (l - xx) +x 2 (l- x 2 )}, 

, . or 2xj-2x 2 + l 2x?- 2a;i + l, 

tts^i,^) = 8{ + }, 

1 — x x 1 — x 2 
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a 5 (xi,x 2 ) 
Pi(xi,x 2 ) 

P2(X 1 ,X 2 ) 



li(xi,x 2 ) 
7 2 (xi,x 2 ) 

73(^1,^2) 
74(^1, x 2 ) 



32(xi +x 2 - 1), 
8(2- Xl -x 2 ), 
Jxi + x 2 )(xj + x 2 2 ) 



3(xi + x 2 ) 2 + 6(xi + x 2 ) -4 



8 



(2 - X! - x 2 ) 

X 



(1 - xi)(l - x 2 )(2 - xi - x 2 ) 
{2(xi + x 2 )(a:? + z 2 ) - 4 ( a; i + x 2 )(2a% + x x x 2 + 2x\) 
+ (19x 2 + 30xia; 2 + 19x 2 ) - 12(2xi + 2x 2 - 1)}, 

^(1 -xi)(l -x 2 )(xi +x 2 - l)(xi - x 2 ) 



(2 - x 2 - xi) 



= 32, 



16, 



(xi + x 2 • 


-1) ( 




- x 2 ) 


(xi + x 2 ■ 






-x 2 ) ( 


(xi + x 2 - 





(xi +x 2 - l)(x 2 - xi), 



(2 - xi - x 2 )(x 2 - xi). 
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FIG. 1. Historical progress of LFV searches for various processes of muons and kaons. 
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FIG. 3. Constraints on the mass of Wr vs. its mixing angle (£) in the manifest left-right 
symmetric model. The experimental constraints of "Strovink" , "Peoples" , "MEGA" and "DO" are 
from Jodidio, et al. (1986), Derenzo (1969), the MEGA experiment (unpublished), and Abachi, et 
al. (1996), respectively. The aimed goal for E614 is also shown (provided by D.R.Gill). 
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FIG. 4. Decay probability distribution of fi — ► e 1/F7 decay as a function of the e energy 
(x = 2E e /m fl ) and photon energy (y = 2E 1 /m tJi ) for unpolarized muons. The decay probability is 
high at x ~ 1 and y ~ 0. 
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FIG. 5. Differential branching ratio of the fj^ — ► e ± uuj decay as a function of the photon 
energy (y = 2E 1 /m^). This branching ratio is obtained by integrating over the e + energy and the 
angle between an e + and a photon. 
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FIG. 6. Photonic penguin diagrams for [i — e transitions, such as fi + — > e + e + e~ or fi~—e~ con- 
version. The cases of (a) a heavy particle {^heavy) hi the loop, and (b) a light fermion (ipn g ht) m 
the loop are shown. $ is a scalar field. 




FIG. 7. Feynman diagram for /i + — > e + 7 decay induced by slepton flavor mixing (Amy. 
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FIG. 8. Feynman diagrams which contribute to the renormalization effect on the slepton masses 
from the Planck to the GUT energy scales. 




FIG. 9. Feynman diagrams for the / u+ -► e+7 decay in SU(5) SUSY GUT. The closed blobs 
represent the flavor transitions due to the off-diagonal terms of the slepton mass matrices. 
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FIG. 10. Predicted branching ratios for the /i + — ► e + 7 decay in the SU (5) SUSY GUT based 
on the minimal supergravity model as a function of the right-handed slepton mass for four different 
sets of the SUSY input parameters of M2 (the SU (2) gaugino mass) and tan j3 (the ratio of the two 
Higgs vacuum expectation values). For the other parameters, the trilinear scalar coupling constant 
Aq = and nit = 175 GeV. The following CKM matrix elements are used: ((Vc^m)^! = 0.04 and 
I {VcKM)td\ = 0.01. (a) and (b) correspond to a positive and negative sign of the higgsino mass 
parameter fj,, respectively. 



SS 




FIG. 11. Predicted /i + — > e + 7 branching ratios in the SU(5) SUSY GUT model with higher 
dimensional operators in the GUT superpotential. The branching ratios are shown as a function 
of the right-handed selectron mass for tan (5 = 6 (dashed lines) and 30 (solid lines). The thick 
lines are for the non-minimal case in which Ug and V] are the same as Vckm, and the thin lines 
are for the minimal case in which Vg = Vckm an d V\ = 1, where Vg is the mixing matrix for 
the right-handed sleptons, and V\ is that for the left-handed sleptons. The bino mass of M\ = 60 
GeV/c 2 , the trilinear scalar coupling constant of A$ = 0, the positive higgsino mass (// > 0), and 
the top quark mass of 175 GeV/c 2 are used. The experimental bound shown in the dashed line is 
B(n + — ► e + 7) < 4.9 x 10~ n (Bolton, et al. (1988)), and it is noted that the recent best limit is 
B(/j, + — > e + 7) < 1.2 x 10~ n (Brooks, et al. (1999)). For detail on the calculations, see Hisano, et 
al. (1998b) (after Hisano et al, (1998b)). 
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FIG. 12. Predicted branching ratios for the fi + e+e+e" decay in the SU(5) SUSY GUT 
based on the minimal supergravity model. The input parameters are the same as in Fig. 9. 
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FIG. 13. Predicted branching ratios for the fi~— e~ conversion in SU(5) SUSY GUT. The 
SUSY gaugino mass of M\ = 50 GeV, and the top Yukawa coupling of /t(M) = 2.4 are used. The 
left and right figures correspond to a positive and negative sign of the higgsino mass parameter 
H, respectively. The experimental bound shown is B(/i~Ti — > e~Ti) < 4.3 x 10 -12 (Dohmen, et 
al. (1993)), and it is noted that the recent best limit is B(fi~Ti — > e~Ti) < 6.1 x 1CU 13 (Wintz 
(1998)) (after Hisano et al. (1997)). 




FIG. 14. Feynman diagrams in SO (10) SUSY GUT which give dominant contributions to the 
fi + — > e + 7 process. (m~.)iiL and (rn~)LR are proportional to m T . 
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FIG. 15. Predicted branching ratios for fx + — > e + 7 decay in the S'O(IO) SUSY GUT based on 
the minimal supergravity model. Input parameters are the same as in Fig. 9. 
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u^ey in the MSSMRN with the MSW large angle solution u->ey in the MSSMRN with the MSW small angle solution 




FIG. 16. Predicted branching ratios of fi + — > e + 7 decay as a function of the Majorana mass of 
the second-generation right-handed neutrino (M U2 ) in the MSSM model with right-handed neutrino. 
They are given for the MSW large angle solution and the MSW small angle solution. The three 
curves correspond to tan/3 = 30, 10, and 3 from top to bottom for both figures. The other 
parameters are shown in the top of the figures (after Hisano and Nomura, (1999)). 
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FIG. 17. Tree diagrams for LFV processes in SUSY models with imparity violation. 
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FIG. 18. One-loop diagrams for LFV processes in SUSY models with i?-parity violation. 
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FIG. 19. Angular distribution of e + in polarized // + — * e + 7 decay. 
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FIG. 20. Examples of theoretical models to induce Mu — Mu conversion. They are mediated 
by (a) a doubly charged Higgs boson, (b) heavy Majorana neutrinos, (c) a neutral scalar particle 
like a tau sneutrino, and (d) a bilepton X (after Willmann, et al, (1999)). 
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FIG. 21. Schematic view of the E614 detector at TRIUMF (provided by D.R. Gill). 
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FIG. 23. Integrated rates of backgrounds from annihilation-in-flight (a dotted line) and radia- 
tive muon decay (a dashed line) as a function of the photon energy. The sum of the two is shown 
by the solid line. 
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FIG. 24. Angular distribution of the physics background from the // + — ► e + v e V^ decay from 
polarized muons (a solid line). /U + — ► e^7 (a dotted line) and fj, + — ► e^7 (a dashed line) decays 
are also shown (after Kuno and Okada, (1996)). 
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FIG. 25. Suppression factor of the accidental background in a /i + — > e + 7 search as a function 
of (half) the detector opening angle. The solid (dotted) line is for 100% (97%) muon polarization 
(after Kuno et al, (1997)). 
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FIG. 26. Schematic layout of the MEGA detector (provided by R. Mischke). 
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FIG. 27. Kinematics of the /i + — > e + e + e~ decay in the muon center-of-mass system, in which 
Pi, f>2 are the momentum vectors of the two e + s and ps is that of the e~ , respectively. The plane-I 
is the decay plane on which p±, p2, and p^ lie. The plane-II is the plane in which the muon 
polarization vectors, P and p$, are located (after Okada, et a/., (1999)). 
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FIG. 28. Fraction of muon decay in orbit normalized to the total nuclear muon capture rate 
as a function of the e~ energy for a titanium target. It represents an effective branching ratio of 
muon decay in orbit as a background to the fi~—e~ conversion. It was calculated by Shanker's 
formula in Eq.(147). The energy of the pT — e~ conversion signal in a titanium target is ^ e =104.3 
MeV. 
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FIG. 30. Electron momentum distributions for the [i~ +Ti — > e~ +Ti reaction, measured by the 
SINDRUM-II detector, after the consecutive analysis cuts. The expected signal at = 4 x 10~ 12 
is shown (provided by P. Wintz). 
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FIG. 31. Schematic layout of the MECO detector (provided by W.R. Molzon). 
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FIG. 32. Positron momentum spectra of the \i~ + Ti — ► e + + Ca reaction. /x~e + (gs) and 
/i _ e + (gr) are the expected signals for the transitions to the ground state and to the giant dipole 
resonance states, respectively (provided by P. Wintz). 
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FIG. 33. Time dependence of the probability of anti-muonium decay when a pure muonium 
atom is created initially. The solid line represents the exponential decay of muonium. The decay 
probability of anti-muonium is given for G m j^/Gf = 1000 (G m -^/Gf = 1.8 x 10~ 2 ) for the dotted 
line (for the dashed line). In the latter case, the maximum probability occurs at about twice the 
muon lifetime (after Willmann and Jungmann, (1998)). 



107 



effective Hamiltonian form: 

(V+A)x(V+A)or(S+P)x(S+P) 
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FIG. 34. Mu — Mu conversion rate for different interactions as a function of the external mag- 
netic field (after Willmann and Jungmann, (1998)). 
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FIG. 35. Schematic layout of the detector for muonium-antimuonium conversion at PSI (after 
Willmann et ai, (1999)). 
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FIG. 36. Distribution of the distance of closest approach between the e + track and e~ track 
versus their timing difference in the experiment to search for Mu — Mu conversion, (a) Monte Carlo 
data and (b) real data (after Willmann et al, (1999)). 
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FIG. 37. Schematic layout of a high-intensity muon source. 
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